Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,233,680	T→C	intergenic (‑447/‑2227)	lysC ← / → yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,233,680	0	T	C	98.3%	190.8 / ‑5.2	60	intergenic (‑447/‑2227)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (34/25);  minor base G (0/1);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCG  >  NC_000913_3_bme_pgi/4233536‑4233823                                                                                                                                                 |                                                                                                                                                  gtaCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAGCAGTACGCTAAGGGCGCCaaaa                                                                                                                                               >  1:443229/1‑149 (MQ=255)   aCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAAcc                                                                                                                                             <  1:297447/149‑1 (MQ=255)    cATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCa                                                                                                                                            <  1:897660/149‑1 (MQ=255)     aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                           <  2:613736/149‑1 (MQ=255)     aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                           <  2:169444/149‑1 (MQ=255)     aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCACCGTGAATTATGATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                           <  1:534605/149‑1 (MQ=255)      tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATc                                                                                                                                          <  1:50402/149‑1 (MQ=255)        tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCCCCAAAACCATcac                                                                                                                                        >  1:54193/1‑149 (MQ=255)        tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATccc                                                                                                                                        >  1:17231/1‑147 (MQ=255)        tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATcac                                                                                                                                        <  1:525026/149‑1 (MQ=255)                  gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCtgt                                                                                                                              <  2:351856/149‑1 (MQ=255)                  gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAACAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCtgt                                                                                                                              <  2:283652/149‑1 (MQ=255)                          ggggCGGTTGTCAACGATGGGGTCATGCGGAGTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCACAACCATCACATTTGTCTTTGACgggtg                                                                                                                      >  2:318410/1‑147 (MQ=255)                                     cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                           <  2:865981/149‑1 (MQ=255)                                      aaCGATGGGGTCATGCGGATTTTTCATCCGCTCCTGGCGGTCAGAAGTTCCGCTGATTAATGCTTCACTGTGCTAAGGGTTTACACACAACATTACGCTAACGACACCAAAACCATCACATTTTTATGTGACTGGTGCTACAATGTTcc                                                                                                          >  1:845396/1‑149 (MQ=255)                                      aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGCCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCCACAATCTTcc                                                                                                          >  2:858098/1‑149 (MQ=255)                                         gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                       <  1:681660/149‑1 (MQ=255)                                          aTGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAACTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                      <  1:996369/149‑1 (MQ=255)                                           tGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGt                                                                                                     >  2:557109/1‑149 (MQ=255)                                                tCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACaa                                                                                                >  2:959535/1‑149 (MQ=255)                                                     cGGATTTTTCATCCACTCCTGGCGGACAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTCCAATCTTCCAGATTAACAATtcta                                                                                           >  2:302967/1‑148 (MQ=255)                                                        aTTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                                                        <  1:488669/149‑1 (MQ=255)                                                         tttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaaaa                                                                                       >  2:46183/1‑149 (MQ=255)                                                            ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                    >  2:428084/1‑149 (MQ=255)                                                               aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATc                                                                                                                                          <  2:539683/92‑1 (MQ=255)                                                               aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATc                                                                                                                                          >  1:539683/1‑92 (MQ=255)                                                                  cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACCTTTATTGAATACCTGTGTGACAATCTTCCAAGAACGCAATTCTCAAAATCAGAagag                                                                              >  1:976229/1‑149 (MQ=255)                                                                  cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCAACAATCTTCCAAAGTCACAAATCTCAAAAACAGAagag                                                                              >  1:309844/1‑149 (MQ=255)                                                                        tGGCGGTCGGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATAAGAAGAGTATTGc                                                                        >  2:27825/1‑149 (MQ=255)                                                                        tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                        >  2:658065/1‑149 (MQ=255)                                                                             gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATg                                                                 <  2:165802/149‑1 (MQ=255)                                                                             atCAGTAGGTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATg                                                                 <  1:221222/148‑1 (MQ=255)                                                                                  tagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAag                                                            >  1:630181/1‑149 (MQ=255)                                                                                  tagtTCAGCTAATAAATGCTGCACTGCGCTAAGGTTTTACACTCAACATTACGCTAACGGCACCACAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCATAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAag                                                            <  2:443229/149‑1 (MQ=255)                                                                                    gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAagag                                                          >  1:702211/1‑149 (MQ=255)                                                                                     ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAgaga                                                         >  1:848708/1‑149 (MQ=255)                                                                                         gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGAcgcg                                                     >  2:867417/1‑149 (MQ=255)                                                                                             ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGCCGccacaa                                                 >  2:996933/1‑147 (MQ=255)                                                                                                aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGc                                              <  1:959535/149‑1 (MQ=255)                                                                                                     ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAGATCAGAAGAGTATTGTCA‑‑ATGGCAAGAGATGCGACGAAGCTTGaa                                         <  2:181782/149‑1 (MQ=255)                                                                                                     ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGaa                                         >  2:630179/1‑149 (MQ=255)                                                                                                          tGCGCTAAGGGTTTACACTGAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                    >  2:319882/1‑149 (MQ=255)                                                                                                                  gggTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCTA‑‑ATGGCAAGATACGCGACGAAGCTTGAAGCAACGGTCGCaa                            >  2:491993/1‑149 (MQ=255)                                                                                                                   ggTTGACACTCAACATTACGTTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCaaa                           <  2:378441/149‑1 (MQ=255)                                                                                                                     tttACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAACGACACCAATCTCAACAT‑‑GAATGGCAGTACCAGCAGGGGATGAGAGGCGACGAAGCTTGACGCAACGGTCGCagatc                         >  1:444572/1‑145 (MQ=255)                                                                                                                         cacTCAACATTACGCTATCGGCTCGAAGACCATCACATTTTTCTGTGACTGGTGGTTCAATCTTCCAAAGTCGCAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCTAAGCTGaa                     <  2:845396/149‑1 (MQ=255)                                                                                                                          acTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                         <  2:198987/144‑1 (MQ=255)                                                                                                                          acTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                         >  1:198987/1‑144 (MQ=255)                                                                                                                           cttaaCATTACGCTAACGGCACCAAAACCATCACATGTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAACCCCAAAAAGCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                   <  2:418431/146‑1 (MQ=255)                                                                                                                            tCAACATTTCGCTAACGGCACCAAACCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                  <  1:179536/149‑1 (MQ=255)                                                                                                                            tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                  >  1:804889/1‑149 (MQ=255)                                                                                                                            tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGATTGAACCAACGGTCGCAAAGCTGaaaaa                  >  2:578824/1‑149 (MQ=255)                                                                                                                            tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTATCAAAATCAGAGGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGACGAAACGGTCGCAAAGATggcaga                  >  2:23733/1‑144 (MQ=255)                                                                                                                              aaCATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                >  2:448676/1‑149 (MQ=255)                                                                                                                                  ttACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAg                                               >  2:135550/1‑114 (MQ=255)                                                                                                                                  ttACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAg                                               <  1:135550/114‑1 (MQ=255)                                                                                                                                      gCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCGAAGCTGAAAAAGCATTGGGCg        >  1:764003/1‑149 (MQ=255)                                                                                                                                        tAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTACAAAGTCACAAGTCTAAAAATTAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGAAACGCTGAAAAAGCAttagcggga      >  2:166051/1‑142 (MQ=255)                                                                                                                                           cGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc   <  1:587199/149‑1 (MQ=255)                                                                                                                                            ggCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCg  >  1:617681/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCG  >  NC_000913_3_bme_pgi/4233536‑4233823

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑