Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,069	0	T	C	41.5%	88.9 / 32.1	65	G517G (GGA→GGG)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base T (29/9);  new base C (0/27);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 7.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGG  >  NC_000913/1575928‑1576218                                                                                                                                              |                                                                                                                                                      aCGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCa                                                                                                                                              <  2:211464/151‑1 (MQ=255) aCGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCa                                                                                                                                              <  1:420780/151‑1 (MQ=255)      cAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAat                                                                                                                                         >  1:121247/1‑151 (MQ=255)      cAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAat                                                                                                                                         <  1:77655/151‑1 (MQ=255)               tCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGc                                                                                                                                 <  2:670291/150‑1 (MQ=255)                  tACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCAt                                                                                                                              >  1:971241/1‑150 (MQ=255)                     cttttttCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCTTTTTAAAATCCTGATAACGCGTCAGAGTTTTATTCCCTGATACCGATTCCGTCCAGCTATGGCGGATATCACGTAACCCCATATAAAAAATATTCCCGCCAt                                                                                                                              <  2:967575/144‑1 (MQ=255)                     cGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGt                                                                                                                           >  2:820263/1‑150 (MQ=255)                        ttttCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACg                                                                                                                       >  1:797787/1‑151 (MQ=255)                           tCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGtt                                                                                                                    >  2:939916/1‑151 (MQ=255)                               ctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTtata                                                                                                                 >  1:583754/1‑150 (MQ=255)                                tgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAAc                                                                                                               >  1:836673/1‑151 (MQ=255)                                   aaccccATTGCGAGTTCATCGTTATAAGGCTTTTTAAAACCTTAATACCGTGTCAGTTTTTTATTCCCTAACCCCAATCCCGCCCACCTATTCCGGATATCATTTAACCCCAAATCAAAAATTTCCCCGCCATAGTAACGGTTATaa                                                                                                                <  2:422405/144‑1 (MQ=255)                                      cccattcccattTCATCTTAATAAGTCTTTTTAAAATCCTAATACCGCGTCAGATTTTTATTCCCTGATCCCTTTTCCGTCCAGCTATTGCGGATACCCCGTAACCCAATATCAAAAATATTCCCGCCATAGTAACGGTTATAACct                                                                                                             <  1:431245/147‑1 (MQ=255)                                      cccATTGCCAGTTCTTCTTTATAAGGCGTTTTCAATTCCTGATAACGCGTCAGATTTTTATTCCCTGATCCCGATTCCGCCCAGCTATTGCGGATATCACGAAACCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACctgc                                                                                                           <  1:434079/149‑1 (MQ=255)                                        cATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGt                                                                                                        >  1:633729/1‑150 (MQ=255)                                             cctgTTCATCTTTATAAGGCTTTTTAAAATCCTGATAACGCGTCAGATTTTTTTTCCCTGATACCGTTCCCGCCCAGTTTTTGCGAATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCa                                                                                                   <  1:806287/147‑1 (MQ=255)                                                tttcttcTTTAAAGGGCTTTTTAAAATCTTAAAAACGCTTAAGATTTTTTTTCCCTAATCCCATTTCCTCCCACCTTTTGGGGATATCATAAAACCCCAAATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTg                                                                                                <  2:669105/149‑1 (MQ=255)                                                 ttctttgtaataCGGTGTTTTCAACTCCTAATAAGGCTTCAGATTTTTATTACCTGATCCCGATCCCTCCCACCTATTGCGGATATCACGTACTCCCATATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTg                                                                                                <  2:721332/140‑1 (MQ=255)                                                     tttttaaagggCTTTTTAAAACCCTAAAAACGCGTAAGATTTTTTTTCCCTGATACCGATTCCGTCCAGTTTTTGCGGATATCACTTAACCCCATAAAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGtt                                                                                           <  1:922808/141‑1 (MQ=255)                                                     tttttATAGGGTTTTTTAAACTCTTAATACCGCGTAAGTTTTTTATTCCCTGATCCCGATCCCCCCCACTTTTTGGGAATATCTCTTAACCCAATATAAAAATTATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAgg                                                                                             <  2:132111/145‑1 (MQ=255)                                                        ttATAGGGCTTTTAAAATTCTTAATAACGCGTCAGTTTTTTATTCCCTGATACCGATCCCGCCCATCTTTGGGGGATATCACATAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTg                                                                                         <  2:655433/149‑1 (MQ=255)                                                         tataAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                      >  2:248833/1‑151 (MQ=255)                                                          ataGGGCTTTTTAAAATCGTGAAACCGCGTAAGAGTTTTATTCCCGAATCCCGATTCCGCCCAGCTATGGCGGATACCACTAAATCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGa                                                                                        <  1:242099/148‑1 (MQ=255)                                                          ataAGGCTTTTTAAAATCCTGATACCGCGTCAGATTTTTATTCCCTGATCCCGATTCCGCCCAGTTATTGCGGATATCACGTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTa                                                                                     <  1:22949/151‑1 (MQ=255)                                                           taAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                      >  2:30912/1‑149 (MQ=255)                                                           taAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                      >  2:243873/1‑149 (MQ=255)                                                             aGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                   >  2:723776/1‑150 (MQ=255)                                                               gctttttAAATTCTTAATAACGCGTCATTTTTTTTTTCCCTGATCCCGTTTCCGTCCAGCTATGGCGAATCTCACGTAACCCAATATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAg                                                                                    <  2:665068/144‑1 (MQ=255)                                                                cGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaa                                                                                >  2:600758/1‑150 (MQ=255)                                                                cGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaa                                                                                >  2:600882/1‑150 (MQ=255)                                                                 tttttAAATTCCTAATACCGCGTAAGTTTTTTATTCCCTAATCCCGTTCCCCCCCAGTTTTTGGGAATATCATTTAACCCAATATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                               <  2:470998/149‑1 (MQ=255)                                                                 gTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                               >  1:620220/1‑150 (MQ=255)                                                                  ttttCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAat                                                                             >  2:650614/1‑151 (MQ=255)                                                                    ttCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAatat                                                                           >  1:828292/1‑151 (MQ=255)                                                                              gATACCGCGTCAGATTTTTATTCCCTGATCCACTCCCCGCCCAGCTATTGCGAATATCACGTAATCCCAAATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCg                                                                  <  1:290725/150‑1 (MQ=255)                                                                               aaagcgctTCAGTTTTTTTTTCCCTAATCCCGATTCCGCCCAGTTTTTGCGAATATCACGTAACCCCATATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGt                                                                 <  1:49515/146‑1 (MQ=255)                                                                               aaaccgctTCAGATTTTTTTTCCCTGATACCGATCCCGCCCAGCTTTGGCGAATATCACATAACCCAATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCAtt                                                                     <  2:722905/142‑1 (MQ=255)                                                                               aTACCGCGTCATATTTTTTTTCCCTGATCCCGATCCCGTCCAGCTATTGCGGATATCACTTACCCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGt                                                                 <  2:598749/150‑1 (MQ=255)                                                                               aTAACGCGTCAGATTTTTTTTACCTGATCCCGATCCCGCCCAGCTTTGGCGGATACCACGTAACCCAAAAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCg                                                                  <  1:908674/149‑1 (MQ=255)                                                                               aTAACGCGTCAGATTTTTATTCCCTGATCCCGTTCCCGCCCACTTTTGGCGGATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTcc                                                                   <  1:305033/148‑1 (MQ=255)                                                                                             ttttttTCCCTAATACCGTTCCCGTCCAGCTATTGGGGATATCACTTAATCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCggg                                                    <  1:262876/149‑1 (MQ=255)                                                                                             ttttATTCCCTGATCCCGTTCCCGTCCAGCTATTGCGGATATCACTTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgg                                                     <  1:799975/148‑1 (MQ=255)                                                                                               ttattaCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgggg                                                   >  1:210451/1‑148 (MQ=255)                                                                                                tattcCCTAATCCCAATCCCGCCCAGTTATTGCGATTATAACGAAACCCAATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGgaga                                                <  2:357294/150‑1 (MQ=255)                                                                                                tattaCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAACCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGgaga                                                <  2:973239/150‑1 (MQ=255)                                                                                                     ccTGATCCCTATCCCGTCCAGCTATTGCGAATATCCCGTAATCCAATATAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATAtt                                            <  1:786012/149‑1 (MQ=255)                                                                                                     ccTGATCCCGATCCCGCCCAGCTATTGCGAATATCAATTAACCCAATAAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGatat                                             <  2:1023003/148‑1 (MQ=255)                                                                                                         aTACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGtt                                      >  2:1013968/1‑151 (MQ=255)                                                                                                          tACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGttt                                     >  1:738116/1‑151 (MQ=255)                                                                                                            ccgttcccGTCCAGCTTTTGCGGATACCCCTTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                    <  1:1001063/144‑1 (MQ=255)                                                                                                             cGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGa                                   >  2:511714/1‑150 (MQ=255)                                                                                                              gATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                                  >  1:206169/1‑150 (MQ=255)                                                                                                                  ccGTCCAGCTATTGGGGATACCACTTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATa                             <  2:443749/151‑1 (MQ=255)                                                                                                                   cGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGAt                              >  2:316189/1‑149 (MQ=255)                                                                                                                     tCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGtt                          >  1:310048/1‑151 (MQ=255)                                                                                                                     tCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGtt                          >  2:739998/1‑151 (MQ=255)                                                                                                                         gttATTGCGGATATCACGTAACCCAATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgt                        <  1:293199/147‑1 (MQ=255)                                                                                                                             tggggaatatCACATAACCCAAAAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatag                   <  2:298256/144‑1 (MQ=255)                                                                                                                             tggCGAATCTCACTAAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcata                    <  1:37378/147‑1 (MQ=255)                                                                                                                                 ggATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtcat               >  2:64043/1‑150 (MQ=255)                                                                                                                                       cagttAACCCCATAAAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCAc        <  1:939940/147‑1 (MQ=255)                                                                                                                                          gTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGcc     >  1:172525/1‑151 (MQ=255)                                                                                                                                            aaTCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCgg   >  1:250206/1‑151 (MQ=255)                                                                                                                                             aTCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCggg  >  1:294655/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                      ACGCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGG  >  NC_000913/1575928‑1576218

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