Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,492,834	0	T	G	43.9%	82.5 / 13.2	59	G89G (GGT→GGG)	tsgA	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base G (25/0);  total (38/21)
Fisher's exact test for biased strand distribution p-value = 1.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGAC  >  NC_000913/3492694‑3492981                                                                                                                                               |                                                                                                                                                    tCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑CTGCC‑CGgttgttgtt                                                                                                                                             >  1:920276/1‑149 (MQ=255)      tATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑GTTTC‑GGCTTGTTGTTTTGCcac                                                                                                                                      >  1:517646/1‑151 (MQ=255)       aTGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCCTGGTTTAAtttccg                                                                                                                                           >  2:40012/1‑141 (MQ=255)          aGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGTGGTTTAGCCACAAcc                                                                                                                                  >  2:544783/1‑151 (MQ=255)           gTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGGC‑GGGTTGGTGGTTAGCCACAGc                                                                                                                                   >  2:868610/1‑149 (MQ=255)              aTACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑GTGCC‑GGGTTGATGTTCAGCCACAGGCTgg                                                                                                                               >  1:337632/1‑150 (MQ=255)                accttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTTAGCCCCAGCCTGgcg                                                                                                                             >  2:887659/1‑150 (MQ=255)                accttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGgcg                                                                                                                             <  2:977082/150‑1 (MQ=255)                accttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTTAGCCCCAGCCCgg                                                                                                                               >  2:286812/1‑148 (MQ=255)                     caccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGtt                                                                                                                        <  2:974069/150‑1 (MQ=255)                        cttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTtctc                                                                                                                     >  2:110172/1‑150 (MQ=255)                           ccTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGGC‑GGTTTTATGTTCAGCCACAACCTGGCGCTGTTCTcggcg                                                                                                                 >  1:222066/1‑151 (MQ=255)                           ccTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTcggc                                                                                                                  <  1:38989/150‑1 (MQ=255)                             tCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTTAGCCACAGCCTGGCGCTGTTCTcggcgg                                                                                                                >  2:190784/1‑150 (MQ=255)                               aaCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGa                                                                                                             <  1:953268/151‑1 (MQ=255)                                    cGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTTATGTTCAGCCACAGCCTGGCGCTGTTTTCGGCGGCGGTGtt                                                                                                         >  2:488463/1‑150 (MQ=255)                                     ggCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCa                                                                                                       <  1:697936/151‑1 (MQ=255)                                      gCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTTGTGTTCAGCCACAGCCTGGCGCTGTTTTCGGCGGCGATGTTCAt                                                                                                      >  1:542881/1‑151 (MQ=255)                                          tttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTCGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATTTTTAGCCACAGCCCGGCGCTGTTCTCGGCGGGGGTGtt                                                                                                         >  1:237814/1‑144 (MQ=255)                                           ttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATtctc                                                                                                  <  2:604931/150‑1 (MQ=255)                                           ttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTTAGCCCCAGCCCGGCGCTGTTTTCGGCGGCGGTGTTTATtttc                                                                                                  >  1:47506/1‑150 (MQ=255)                                           ttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑CGGTTGATGTTCAGGCACAGCCTTGGGCCGTTCTCGGCGGGGGTGtt                                                                                                         >  2:317447/1‑143 (MQ=255)                                            tAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCCCAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATtct                                                                                                   <  1:115987/148‑1 (MQ=255)                                                   tATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGGTGTTCAGCCACAGCCTGGCGCTGTTCTTGGCGGGGCTGTTTCTTCTCGgggggg                                                                                           >  1:311840/1‑149 (MQ=255)                                                       ttCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTCATTTTCGGGGGGGGCAGCg                                                                                     >  1:248272/1‑151 (MQ=255)                                                         ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGc                                                                                   <  2:779124/151‑1 (MQ=255)                                                         ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGc                                                                                   <  1:611690/151‑1 (MQ=255)                                                         ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGc                                                                                   <  1:342253/151‑1 (MQ=255)                                                         ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTTAGCCCCAGCCTGGCGCTGTTCTCGGCGGCGGTGTTCATTCTCGGGGGGGTCAGCgg                                                                                    >  1:695051/1‑150 (MQ=255)                                                          cTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTTTCGGCGGCGGTGTTTATTCTCGGGGGGGGCAGCgg                                                                                    >  1:802241/1‑149 (MQ=255)                                                          cTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGTTGTTCAGCCACAGCCCGGCGCTGTTCTCGGCGGGGAGGTTTATTTTCGGGGGGGGCAGCgg                                                                                    >  2:117096/1‑149 (MQ=255)                                                          cTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑CTTCC‑GGGTTGATGTTTAGCCACAGCCCGGCGCTGTTCTCGGCGGCGGGGTTCATTCTCCGGGGGGGCCGCgg                                                                                    >  2:203515/1‑149 (MQ=255)                                                              aCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTTAGCCACAGCCTGGCGCTGTTTTCGGCGGCGGTGTTCATTCTCGGGGGGGTCAGCGGCATCAc                                                                              >  1:197955/1‑151 (MQ=255)                                                                      cTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCC‑GGGTTGATGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTCATTTTCGGGGGGGTCAGCGGCATCACCATGTCg                                                                       >  1:711261/1‑150 (MQ=255)                                                                         aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTg                                                                   <  2:216519/151‑1 (MQ=255)                                                                         aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTg                                                                   <  2:343187/151‑1 (MQ=255)                                                                         aTGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTTATTCTCGGGGGGGTCAGCGGCATCACCCTGTCGCtt                                                                    >  2:909415/1‑150 (MQ=255)                                                                                tCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTTATGTTCAGCCCCAGCCTGGCGCTGGTCTCGGCGGCGGTGTTCATTTTCGGGGGGGGCAGCGGCCTCACCCTGGCGGTTGGGGCAtt                                                            >  2:173946/1‑151 (MQ=255)                                                                                    cccGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTTATGTTCAGCCACCGCCTGGCGCTGGTCTCGGCGGCGGTGGTCCTTCCTGGGGGGGTCAGCGGCATCACCCTGTCGATTGGTACAtttc                                                          >  2:204121/1‑147 (MQ=255)                                                                                     ccGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATTTTCAGCCACAGCCTGGCGCTGTTCTCGGGGGCGATGTTCATTTTCGGGGTGGTCAGCGGCATCACCCTGTCGATTGGTTCATTCCTgg                                                       >  1:474280/1‑151 (MQ=255)                                                                                        ttGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCGCAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCGCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTa                                                     <  2:572855/150‑1 (MQ=255)                                                                                          gAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑CGGTTGATGTTTAGCCCCAGCCCGGCGGTGTTTTCGGCGGCGGTGGTTTTTCTCGGGGGGGGCAGCGGCCTCACCCTTTCGGTTGGGGCATTTCttg                                                       >  2:503750/1‑144 (MQ=255)                                                                                           aaaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTGCC‑GG‑TTGATGGTTCGCCCCAACCCCGCGCCGGTCTCTGGGGGGGTGTTTCTTTTCGGGGGGGGCAACGGCATCCCCATGTGGGGTGGTATCTTCCTggggaaa                                                  >  2:461653/1‑145 (MQ=255)                                                                                             aaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACaa                                               <  1:652049/151‑1 (MQ=255)                                                                                                gCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATg                                            <  1:626688/151‑1 (MQ=255)                                                                                                gCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTTATTTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTTTCGGGGGGGGCAGCGGCATCACCACGTCCGTTGGTACATTTCTGGGAAAACaaa                                              >  2:369286/1‑149 (MQ=255)                                                                                                    ttACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCC‑GGGTTGATGTTCAGCCACAACCTGGCGCTGTTTTCGGCGGCGGTGTTCATTTTCGGGGGGGGCAACGGCATCACCATGTCGGTTGGTACATTCCTGGTAACACAAATGTAt                                         >  1:696832/1‑150 (MQ=255)                                                                                                    ttACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGGGATGTTCATTCTCGGGGGGGTCAGCGGCATCACCCTGTCGATTGGTACATTCCTGGGAACACaaaaggat                                         >  1:930881/1‑145 (MQ=255)                                                                                                     tACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTTATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGGGGGCAGCGGCATCACCATGTCGGTTGGTAAATTCCTGGTAACACAAATGTa                                          >  1:175410/1‑148 (MQ=255)                                                                                                      aCGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGaa                                      <  2:964659/151‑1 (MQ=255)                                                                                                        gTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAgg                                    <  2:658251/151‑1 (MQ=255)                                                                                                        gTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAgg                                    <  2:688880/151‑1 (MQ=255)                                                                                                          tttGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCC‑GGTTTGATGTTTAGCCCCAGCCTGGCGCTGTTCTCGGCGGCGGTGTTCATTCTCGGGGGGGTCAGCGGCATCACCATGTCGATTGGGACATTTCTGGTAACACAAAAGTATGAAggg                                   >  1:267897/1‑150 (MQ=255)                                                                                                          tttGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAggg                                   <  2:738400/150‑1 (MQ=255)                                                                                                             ggCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTTATGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTCATTTTCGGGGGGGGCAAGGGCATCACCATGTCGATTGGTACATTTCTTGGAAAAAAAAAAGATGAAGGGCg                                 >  1:908292/1‑149 (MQ=255)                                                                                                               cTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGGTTGATGTTCAGCCACAGCCTGGCGCTGTTTTCGGCGGCGATGTTCCTTCTCGGGGGGGGCAGCGGCATCACCATGTCGATTGGGACATTCCTGGTAAAACAAATGTATGAAGGGCg                                 >  1:640661/1‑147 (MQ=255)                                                                                                                 ttCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAg                             <  2:156957/149‑1 (MQ=255)                                                                                                                                  gCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAt          <  1:892354/151‑1 (MQ=255)                                                                                                                                         gCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCg    >  2:810849/1‑151 (MQ=255)                                                                                                                                             ggTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGAc  >  2:893399/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                    TCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCC‑GGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGAC  >  NC_000913/3492694‑3492981

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