Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,603,923	0	T	G	52.7%	36.2 / 66.2	55	V302G (GTG→GGG)	lsrC	autoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base G (29/0);  total (40/15)
Fisher's exact test for biased strand distribution p-value = 6.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-08

CATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAAC  >  NC_000913/1603774‑1604071                                                                                                                                                      |                                                                                                                                                         cATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGgtg                                                                                                                                                         <  1:1142180/151‑1 (MQ=255) cATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGgtg                                                                                                                                                         <  1:245389/151‑1 (MQ=255)       tttGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGAt                                                                                                                                                   <  2:58590/151‑1 (MQ=255)             cGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAACGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGc                                                                                                                                             <  1:759355/151‑1 (MQ=255)               gtggtTCCGGGGCGATCATTGGTGCGGTACTCGGCGCATGGGTCCTGTCGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGCCTGGTTCTGCGGGCGGTGCTGGGGTTTGATGGAGGcca                                                                                                                                           >  2:415711/1‑150 (MQ=255)                  gtTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCg                                                                                                                                        >  1:984172/1‑151 (MQ=255)                   ttCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGt                                                                                                                                       <  1:1239893/151‑1 (MQ=255)                     ccGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTtg                                                                                                                                     <  1:1056017/151‑1 (MQ=255)                            gATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGGCGCCTTCGTTTTgcg                                                                                                                                 >  1:630041/1‑148 (MQ=255)                               cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGaa                                                                                                                           <  2:457354/151‑1 (MQ=255)                               cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGaa                                                                                                                           <  2:300117/151‑1 (MQ=255)                                    gTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTTATGGACGCCTGCGTTGTGCGCTGGGACGTaa                                                                                                                      >  1:305763/1‑151 (MQ=255)                                    gTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTaa                                                                                                                      >  2:768304/1‑151 (MQ=255)                                    gTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGGGCTGGTGGTTGATGGACGCCTGCGTTGTGCGCTGGAAAGtga                                                                                                                      >  2:552275/1‑149 (MQ=255)                                       cGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATc                                                                                                                    <  1:347932/150‑1 (MQ=255)                                       cGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCt                                                                                                                   <  2:1117565/151‑1 (MQ=255)                                         gTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGGACGGAATCTa                                                                                                                  >  1:123295/1‑150 (MQ=255)                                         gTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTa                                                                                                                  >  2:1136121/1‑150 (MQ=255)                                         gTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTa                                                                                                                  >  1:124953/1‑150 (MQ=255)                                           aCTCGGCGCATGGTTCCTGACGCATATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCAAGGGGGAAGGATTTTACCGCGGGGCTGGTTCTGCTGGCGGTGCCGGGGTTTGATCGTCGGCTGCGGTGTGCGCTGGCAAGGAAACTAcgg                                                                                                               >  2:27286/1‑151 (MQ=255)                                           aCTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTAcgg                                                                                                               <  2:988481/151‑1 (MQ=255)                                              cGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGGCGCCTGCGGTGTGCGCTGGAACGCAATCTAcggc‑gg                                                                                                           >  2:1077907/1‑151 (MQ=255)                                                gcgcATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCa                                                                                                         >  2:244321/1‑151 (MQ=255)                                                gcgcATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGGTTGATGTAAGCCTGCGGTTTGCGCTCGAACGTCATCCAccgc‑ggga                                                                                                         >  2:762482/1‑149 (MQ=255)                                                       gTTCCTGACGCAGCTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCCTGGTGGAATGCTTTTATCGCGGGTCTGGTTCTGCTGGCGGAGCTTGGGTTTGATGGGAGGCTGGGTTGGGTGCTGGGACGGTATCGTCGGC‑GGTCAACatat                                                                                                  >  2:1155279/1‑151 (MQ=255)                                                       gTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTTATGGGCGCCTGCGGTGTGCGCTGGGACGTAATCTACGGCGGGCaaaaaat                                                                                                   >  2:1163888/1‑149 (MQ=255)                                                         tCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAAAATg                                                                                                 >  1:936411/1‑150 (MQ=255)                                                         tCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGGGCTGGGGTTTGATGGGCGCCCGCGTTGTGCGCTTGGACGTAATCTACGGC‑GGCCAAAAAATg                                                                                                 >  2:695158/1‑150 (MQ=255)                                                                 cAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGGGCCTGGGGTTGATTGACGCCCGCGGTGTGCGCTGGAACGTAATTTAAGGC‑GGCCAAAAAAAGGCCC‑Ctttt                                                                                       >  2:694898/1‑150 (MQ=255)                                                                   gATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTTTGCGCTGGAACGTAATCTACGGC‑GGCCAAAAAAAGCCCG‑Ctttttt                                                                                     >  2:1183090/1‑148 (MQ=255)                                                                   gATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGGTGGACGCCTGCGGTTTGCGCTGGGACGTAATCTACGGC‑GGGAAAAAAATGCCCG‑CCTTTT‑g                                                                                    >  1:376872/1‑151 (MQ=255)                                                                   gATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGGGCTGGGGTTTGATGGACGCCCGCGTTGTGCGCTGGAACGTAATCTCCAGC‑GGCAAAAATATGCCCC‑CTTTTT‑g                                                                                    >  2:982611/1‑151 (MQ=255)                                                                   gATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGGGGAATGCTTTTATCGCGGGTCTGGTTCTGCTGGCGGTGGTGGGGTTTGATGGACGCCTGCGTTTTGGGCTGGAACGTAATCCACGGC‑GGGCAAAATATTGCCG‑CTTTTT‑g                                                                                    >  1:567528/1‑151 (MQ=255)                                                                     tCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGGGCTGGGACGGAATTTACGGC‑GGCAAAAAGATGCCCC‑CTTTTT‑gg                                                                                   >  1:583563/1‑149 (MQ=255)                                                                        aTAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGTTGGACGCCTGGGTTGGGCGCTGGGAAGGCATCTACGGG‑GGCCAAAATATTCCCG‑GTTTTT‑gcgcc                                                                                >  1:995534/1‑146 (MQ=255)                                                                             gTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCCGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCg                                                                          >  2:349965/1‑150 (MQ=255)                                                                              tACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGGACGTAATTTACGGC‑GGCCAAAATATGCCCG‑CTTTTT‑TACGCC‑ACCGc                                                                         >  2:240379/1‑150 (MQ=255)                                                                              tACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGGGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATTTACGGC‑GGGAAAAATATGCCCG‑GTTTTT‑GTCGCC‑CCCGc                                                                         >  1:920278/1‑150 (MQ=255)                                                                                 tGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATc                                                                     <  1:1226392/151‑1 (MQ=255)                                                                                  ggTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGGGCTGGAACGTAATTTACGGC‑GGCAAAAATATGCCCG‑CTTTTT‑GACGCC‑CCCGCCCTc                                                                     >  1:311243/1‑150 (MQ=255)                                                                                   gTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTTGGACGGAATCTACGGC‑GGCAAAAATATGCCCC‑CTTTTTGGACGCA‑CCCGCCTTc                                                                     >  2:113447/1‑150 (MQ=255)                                                                                                ttCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGGTGGGCGCTGGAACGTAATCTACGGC‑GGCCAAAATAAGCCCGCCTTTTG‑GAGGCC‑A‑CGCCCTTCGTTTAAACCGcct                                                       >  2:864779/1‑148 (MQ=255)                                                                                                ttCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCCGCGTTGTGGGCTGGAACGGAATCTACGGC‑GGCCAAAAAATGCCCG‑CCTTAT‑GACGCCAACCCCCATCCTTAAAAccccat                                                        >  2:19355/1‑147 (MQ=255)                                                                                                       aTGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAAccc                                                           <  1:78412/139‑1 (MQ=255)                                                                                                       aTGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGGTTTGGGCTGGGACGTAATCTACGGC‑GGGCAAAAAATTCCCC‑CGTTTT‑TGCGCC‑ACCGCCCACCCTTAAAcccccgtttcttcta                                               >  2:78394/1‑139 (MQ=255)                                                                                                         ggtggAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGt                                                                  <  2:825088/130‑1 (MQ=255)                                                                                                         ggtggAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGGTTGATGGACGCCTGCGTCGTGCGCTGGGACGTAATCTAAGGC‑GGCAAAAAAATGCCCG‑CCTTTT‑GACGCC‑AACGCCCTCCGTTAAAACCCCCTCGTCAGGGaa                                             >  2:1086162/1‑151 (MQ=255)                                                                                                          gtggAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTTT‑GACGCC‑CCCCCCCTCCGTTAAACCCGCTTCGGTAGGGaaa                                            >  1:906467/1‑151 (MQ=255)                                                                                                                tGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCCGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCCAAAATATGCCCG‑CTTTAT‑GGCGCC‑CCCGCCATCCGGTAAACCCGCTTCGTCAGGGaaaaaaaag                                      >  2:345094/1‑148 (MQ=255)                                                                                                                  aTTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTTT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGGAAAAAACgggg                                    >  1:278055/1‑150 (MQ=255)                                                                                                                  aTTTTATCGCGGGTCTGGTTCTGCTGGCGGGGCTGGTGTTTGATGGACGCCCGCGGTGTGCGCTGGGACGTCATCTACGGG‑CGGAAAAATATGCCCG‑CCTTTT‑TATGCC‑CCCCCCCACCCGTAAACCCCCTTCGTCAGGGAAAAAAAgggg                                    >  2:689857/1‑150 (MQ=255)                                                                                                                             ggTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAAGGAAATTCACGGG‑GGAAAAAAATTCCCCC‑CTTTTT‑ACGGCC‑CCCCCCATCCGTTAAACCCCCTTTCTCGGGTAAAAAACGGGGGGCCCCATaag                         >  1:1096576/1‑150 (MQ=255)                                                                                                                             ggTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATaa                          <  2:376778/150‑1 (MQ=255)                                                                                                                               tCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATg                        <  2:1036400/150‑1 (MQ=255)                                                                                                                                     tCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTc                 <  2:765686/151‑1 (MQ=255)                                                                                                                                      ctgctgGCGGTGCTGGGGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGGAATCTACTGC‑GGGAAAAAAATGGCCG‑CTTTAT‑GAAGCC‑ACCGCCATCCGGTTAAACCGCTTCGTCAGGGAAAAAAAGGGAGGCCGCATAATGCGCAttt                 >  2:64767/1‑149 (MQ=255)                                                                                                                                      ctgctgGCGGTGCTGGGGTTTGATGGACGCCCGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGGCAAAAAATGCCCC‑CCTTTT‑TAAGCC‑ACCCCCACCCGCTAAAACCCCCTCGTCAGGGAAAAAAAAGGGGGGCCCCTAATTCGTATTc                 >  2:358103/1‑150 (MQ=255)                                                                                                                                                     gtgtTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGaa   >  1:604567/1‑149 (MQ=255)                                                                                                                                                     gtgtTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAAGCCGCATAATGCGTATTCGCTACGGTTGGGAAc  <  2:1259942/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         CATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGCCTGCGTTGTGCGCTGGAACGTAATCTACGGC‑GGCAAAAATATGCCCG‑CTTTAT‑GACGCC‑ACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAAC  >  NC_000913/1603774‑1604071

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