Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,355,457	0	A	C	41.9%	110.3 / 14.5	62	T245P (ACC→CCC)	gltB	glutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base A (28/8);  new base C (0/26);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCG  >  NC_000913/3355314‑3355604                                                                                                                                                |                                                                                                                                                    gcTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGt                                                                                                                                                <  2:717461/149‑1 (MQ=255) gcTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTaa                                                                                                                                              <  2:996099/151‑1 (MQ=255) gcTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAACGCTTCTCCGCTATCACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGt                                                                                                                                                <  1:843727/149‑1 (MQ=255)             aTCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgc                                                                                                                                  >  1:743134/1‑151 (MQ=255)                     gACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc                                                                                                                          >  1:125602/1‑151 (MQ=255)                      aCCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg                                                                                                                         >  2:1058931/1‑151 (MQ=255)                       ccTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcgt                                                                                                                        >  2:760137/1‑151 (MQ=255)                            gTCTGGAATCGGCCTTTTCCCTGTTCCCCCGGCGCTTCCCCACAAACCCCGTACCGCGCTGGCCGCTGGCCCACCCGTTCCGCTTTCGGGCGCATAACGGTGAAATCACCCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACtt                                                                                                                    <  2:493466/150‑1 (MQ=255)                             tCTGGAATCGGCCTTTTGCCTGTTCCCCCAGCGCTTCTCCACTACCCCCGTCCCGCGCTGGCCGCGGGCGCACCCGTCCCGCTATCTGGCGCATAACGGTGAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTta                                                                                                                   <  2:125146/150‑1 (MQ=255)                              cTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTtat                                                                                                                  >  2:130313/1‑150 (MQ=255)                               tGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTtata                                                                                                                 >  2:536140/1‑150 (MQ=255)                                ggAATCGGCCTTTTGCCTTTTCCCCCAGCGCTTCTCCACAAACCCCGTCCCGCGCTGGCCGCTGGCCCACCCTTTCCGCTATCTGGCGCATAACGGTAAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATaaa                                                                                                               <  2:190515/151‑1 (MQ=255)                                ggAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATaa                                                                                                                >  2:918657/1‑150 (MQ=255)                                   aTCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTc                                                                                                            >  2:604914/1‑151 (MQ=255)                                   aTCGGCCATTTCCCTGTTCCACCAGCGTTTCTCCACTAACCCCGTCCCGCGCTGGCCGCTGGCGAACCCGTTCCGCTTTCGGGCGAATAGGGGGAAAAACAACACCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAAtt                                                                                                             <  2:283274/150‑1 (MQ=255)                                     cGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTcc                                                                                                           >  2:985396/1‑150 (MQ=255)                                      ggCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTcc                                                                                                           >  1:411871/1‑149 (MQ=255)                                                cTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                                                                                                >  2:463862/1‑150 (MQ=255)                                                  gTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATc                                                                                             >  1:607389/1‑151 (MQ=255)                                                     cccccaGCGCTTCCCCACAAACCCCGTCCCGCGCTGGCCGCTGGCGCACCCGTTCCGCTATCGGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCt                                                                                          <  2:146339/148‑1 (MQ=255)                                                      caccaGCGCTTCTCCACTAACACCGCCCCGCGCGGGCCGCTGGCGCACCCGTCCCGTTATCGGGCGCATAAGGGTAAAATCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTg                                                                                         <  2:601813/151‑1 (MQ=255)                                                       accaGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTg                                                                                         >  2:948103/1‑150 (MQ=255)                                                         caGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc                                                                                      >  1:346834/1‑151 (MQ=255)                                                          aGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCt                                                                                     >  2:346028/1‑151 (MQ=255)                                                               ttCCCCACTAACCCCGCCCCGCGCGGGCCGCTGGCGCACCCTTCCCGCTATCTGGCGCATAACGGGGAAATCAACACCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCt                                                                                     <  2:472036/146‑1 (MQ=255)                                                                 ctcCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacgac                                                                              >  1:674183/1‑151 (MQ=255)                                                                 ccccACTAACCCCGTCCCGCGCGGGCCGCTGGCGCAACCGTCCCGCTATCGGGCGAATAACGGTAAAATCAACACCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacga                                                                               <  2:991111/148‑1 (MQ=255)                                                                      ctacccccGTCCCGCGCTGGCCGCTGGCGCACCCGTCCCGCTATCTGGCGCATAAGGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacgacg                                                                             <  2:651209/141‑1 (MQ=255)                                                                       tAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGAcgccgc                                                                         >  1:829089/1‑150 (MQ=255)                                                                          ccccGTCCCGCCCGGGCCGCGGGCGCACCCGCACGCATATCTGCCAGATGCCGGTGAACACAACCCCATCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCAc                                                                       <  2:396593/147‑1 (MQ=255)                                                                                cccccgcggcccgctggcCCACCCTTCCCGTTTTCTGGCCCATAACGGTAAAACCACCCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTc                                                               <  2:393217/141‑1 (MQ=255)                                                                                 ccGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTAAAATAAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCa                                                              <  1:411810/151‑1 (MQ=255)                                                                                  cgcgCTGCCCGTTGGCCCTCCCGTCCCGCTACCTGGCGCAAAGCGGTAAAACCAACACCATCACGGTAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCaa                                                             <  1:273660/151‑1 (MQ=255)                                                                                        ggccgcggccgcACCCTTCCCGCTTTCTGGCGAATAACGGTAAAACCAACCCCATCACCGTAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAAc                                                       <  1:225286/151‑1 (MQ=255)                                                                                        ggccgcggccgcACCCGTCCCGATATCGGGCGCATAACGGTGAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACg                                                           <  1:82216/147‑1 (MQ=255)                                                                                          ccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCg                                                     >  2:161166/1‑151 (MQ=255)                                                                                          ccgcgggcGCACCCGTTCCGTTATCTGGCGAATACGGGAAAAACCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAAc                                                       <  2:641303/149‑1 (MQ=255)                                                                                            gctggcGCACCCGTCCCGCTATCTGGCGAATAAGGGTAAAATAAACCCCTCCCCCGTAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCgg                                                    <  2:130060/150‑1 (MQ=255)                                                                                                 cccacccttcccGCTATCGGGCGCATAACGGTAAAATCACCCCCTTCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGct                                                  <  1:286270/137‑1 (MQ=255)                                                                                                   cacccGTCCCGCTATCGGGCGCATAGGGGAAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGa                                              <  1:632173/146‑1 (MQ=255)                                                                                                    aaCCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTc                                           >  1:673895/1‑151 (MQ=255)                                                                                                     aCCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTcc                                          >  1:122075/1‑151 (MQ=255)                                                                                                     aCCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTcc                                          >  1:693374/1‑151 (MQ=255)                                                                                                      ccGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTcc                                          >  1:286605/1‑150 (MQ=255)                                                                                                      ccGTCCCGCTATCGGGCGCAAAGCGTTAAAATCACCCCCACCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCa                                         <  1:1058797/151‑1 (MQ=255)                                                                                                        gTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGt                                       >  1:496598/1‑151 (MQ=255)                                                                                                            cGCTATCTGGCGCATAACGGTAAAATCAACACCATCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGt                                       <  2:309450/147‑1 (MQ=255)                                                                                                            cGCTATCTGGCGCATAACGGTAAAACCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGa                                   <  1:94132/151‑1 (MQ=255)                                                                                                               tATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATgg                                >  2:623856/1‑151 (MQ=255)                                                                                                                aTCTGCCGAATACCGGAAAAACCAACCCCCTCCCGGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGtt                                      <  2:176261/144‑1 (MQ=255)                                                                                                                 tCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGAt                              >  2:649865/1‑151 (MQ=255)                                                                                                                  cTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGAt                              >  2:1100971/1‑150 (MQ=255)                                                                                                                      cgcATAACGGTAAAATCAACACCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACa                          <  2:1101650/150‑1 (MQ=255)                                                                                                                           acgggtaaaacaaaCACCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACa                          <  1:288360/133‑1 (MQ=255)                                                                                                                                taaaaCCAACCCCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAAct                <  1:648858/148‑1 (MQ=255)                                                                                                                                 aaaaTCAACCCCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAAc                 <  1:791848/147‑1 (MQ=255)                                                                                                                                 aaaaTCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActg               <  2:743239/149‑1 (MQ=255)                                                                                                                                  aaaTCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActgct             <  1:827903/151‑1 (MQ=255)                                                                                                                                       aaCCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActgctgc           <  1:1114673/148‑1 (MQ=255)                                                                                                                                         cacCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGc       >  1:1020936/1‑150 (MQ=255)                                                                                                                                           ccATCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAgg    <  2:690681/151‑1 (MQ=255)                                                                                                                                             aTCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCg  >  2:100467/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCG  >  NC_000913/3355314‑3355604

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