Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	47.5%	37.1 / 45.1	61	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (27/5);  new base A (0/29);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 2.02e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357740‑358033                                                                                                                                                             |                                                                                                                                                      cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                                       >  1:1343493/1‑151 (MQ=255)   gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaaca                                                                                                                                                      >  1:1655024/1‑150 (MQ=255)    aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                                     >  1:457800/1‑150 (MQ=255)    aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                                     >  1:1570263/1‑150 (MQ=255)     tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa                                                                                                                                                    >  2:727470/1‑150 (MQ=255)          gTAGCGCATCAGGCAATGCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAA‑CACCCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑ccgcc                                                                                                                                             <  1:1202641/151‑1 (MQ=255)            aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGt                                                                                                                                           >  1:1558092/1‑151 (MQ=255)              tgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTtc                                                                                                                                         >  2:657146/2‑151 (MQ=255)              ccaacagggCATTTCCAGCCCCAATCCGGTGTCACCGGCTCCGGTGATTCCTTTCCGCCAA‑CACCCGCGCTTTTAC‑CC‑AC‑CCCCTGGTGAAAAAAT‑GCTTTCGCGGGGCCTGTCTGCCAGTTATTTTG‑CCGCCGCAACAAAA‑ccgccg                                                                                                                                            <  1:1138059/141‑1 (MQ=255)                 aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTtctc                                                                                                                                       >  1:1967994/1‑150 (MQ=255)                   cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                                    >  2:599383/1‑151 (MQ=255)                       cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                                 >  1:1339555/1‑150 (MQ=255)                           tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGAgg                                                                                                                            >  1:1041513/1‑151 (MQ=255)                             cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTcgcg                                                                                                                          >  1:1292882/1‑151 (MQ=255)                              aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                         >  1:1028765/1‑151 (MQ=255)                                ccGCAGCCTTTTGTCGGGGGCTCCCGTTTTTCATTCCCCCAAA‑CACCCGCGCTTTTAC‑CC‑AC‑CGCCAGGTGAAAAAAT‑CCCTTGGCGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCagtagt                                                                                                                              <  2:1714937/144‑1 (MQ=255)                                ccGCAGACCTGTGTCGGCGGCTCCCGTTTTTCTTTCCCCCCAA‑CACCCGCGCATTTAC‑CC‑AC‑CCCCTTGTGAAAAAAT‑GCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                         <  2:1054142/149‑1 (MQ=255)                                  gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGTCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                         >  1:1382419/1‑147 (MQ=255)                                    aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                      >  2:386408/1‑147 (MQ=255)                                    aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                       >  1:132199/1‑147 (MQ=255)                                       ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                       >  2:1688461/1‑144 (MQ=255)                                        cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                       >  1:1377395/1‑143 (MQ=255)                                         tgtgGCCGCGTCTACCGTTTTTC‑TTTCCCCCAACCACCCCCGCATTTTC‑CC‑AC‑CCCCATGTGAAAAATT‑CCCTGCGGGGCGGCTGTCTCCATGTTTTTTCC‑CCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                  <  1:1594801/147‑1 (MQ=255)                                         gttgacAGCGGCCACCGTG‑TTTATTTCCCCCAACCACCCGCGTTTTTAC‑CA‑AA‑CGCTATGTACAAAAAG‑GCCTGCGGGGCGTCTCTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑aa                                                                                                             <  2:1254335/145‑1 (MQ=255)                                          gtgtCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                       >  1:1582594/1‑141 (MQ=255)                                              cAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggggcgacgacaaa                                                                                                           >  1:1542684/1‑137 (MQ=255)                                               ggcgtttcccGTTATTCTTTCCCCCCAA‑CACCCGCCTTTTTAC‑CC‑AC‑CGCCATGTGAAAAATT‑CCCTCCGGGGCGGCTTTCGCCAAGTTTTTTTC‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑c                                                                                                       <  2:1705003/141‑1 (MQ=255)                                               aGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑aaaaa                                                                                                          >  2:1905856/1‑145 (MQ=255)                                                 cGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGcgcgaggaaaaaattt                                                                                                       >  2:25812/1‑135 (MQ=25)                                                  ggCTCCCGTGTTACTTTCCCGCCAA‑CACCCGC‑CATTTCA‑CC‑AACCGCCATGTGAAAAAAT‑CCTTCCGGGGCGGCTGTCTCCCACTTTTTTTC‑CCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑Ctt                                                                                                     <  2:25050/148‑1 (MQ=255)                                                   gCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggggtgggggaaaa                                                                                                            >  1:842544/1‑132 (MQ=255)                                                    cccccttaattcattcCCCCCAA‑CACCCACAC‑TTTATCCC‑AC‑CGCCTTGTCAAAAAAT‑CCCTTGGGGGGGGCTGTCTCCCAGTTATTTCC‑CCCCCACAACAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTaa                                                                                                  <  2:795593/141‑1 (MQ=255)                                                       ccGTGTTTCATTCCCGCCAA‑CACCCGCGCTTTTAC‑CC‑AC‑CGCCATGTAAAAAATT‑CCCTCCGCGGGGGCTGTCTCCCATTTTTTTTC‑CCGCCGAAAAAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGcc                                                                                               <  2:1567470/149‑1 (MQ=255)                                                            aatCATTTCCCCCCA‑ACACCCCCCATTTAT‑CCAAA‑CGCTATGACAAAAAAT‑GCCTGGGGGGGGGTTCTCGGCCATTTTTTTTC‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGt                                                                                            <  1:1240835/146‑1 (MQ=255)                                                             ttCATTCCCCCAAA‑CACCCGCGCATTAAT‑CC‑AC‑CCCAATGTACAAAAAT‑GCCTCCCCGGCGGCTGTCTCCCATTTTTTTTC‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTc                                                                                           <  2:33222/147‑1 (MQ=255)                                                                   cccccaaa‑caccCGCGTTTTTAC‑CC‑AC‑CGCCATGTCAAAAAAT‑CCCTCGGGGGCGCCTTTCT‑CCAGGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCa                                                                                          <  2:728740/136‑1 (MQ=255)                                                                        caa‑caaCCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAGACTTTAAAccctcctttttt                                                                                     >  2:649615/1‑134 (MQ=25)                                                                         aa‑caccCGCGCTTTTTC‑CC‑AC‑CCCAATGTGAAAAAAT‑GCCTTCGCGGCGGTTGTCTGCCAGTTTTTTTC‑CCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCaa                                                                             <  2:584965/149‑1 (MQ=255)                                                                         aa‑caccCGCGCTTTTAC‑CC‑AC‑CCCCATGTGAAAAAAG‑GTCTCCGCGGCGGCTGTCTGCCATTTTTTTTC‑CCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                           <  2:250872/151‑1 (MQ=255)                                                                         aa‑caacCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGG‑AAAAAA‑CTTTAAccccccctgtt                                                                                       >  2:1412675/1‑133 (MQ=25)                                                                          a‑caacCGCGCTTTTAC‑CC‑AC‑CGCCAGGTGAAAAAAG‑GCCTTCGCGGCGGTGGTCTCCCACTTTTTTTC‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCaa                                                                             <  2:2038821/148‑1 (MQ=255)                                                                                      tttAC‑CC‑AC‑CGCCATGTAAAAAAAC‑CCCTTCGCGGCGGTTGTCTCCCATTTTTTTTC‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGt                                                                <  2:973504/148‑1 (MQ=255)                                                                                       ttAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑A‑AAAAAA‑ATTTTAACCCGGCCTTTTGTGTggggatcccctctggt                                                                  >  2:1169003/1‑133 (MQ=25)                                                                                        tAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑GAAAAA‑ATTTTTTCCCCTCCTGTTTTGTGGGGATCggcgc                                                                      >  2:964223/1‑138 (MQ=21)                                                                                              caaccGCCGTGTAAAAAAAT‑GCGCGCGCGGGGGCTGTCCCCCGCGTATTTTC‑CGCCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAAtg                                                         <  2:1288103/146‑1 (MQ=255)                                                                                                         taaaaaaaT‑CCCTTCGCGGGGGCTTTCTGCCATTTTTTTTC‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAgg                                                 <  2:1224093/144‑1 (MQ=255)                                                                                                         ggaaaaaaTGGCCTCCGGGGCGGTTTTCCGCCATTTTTTTCC‑CCGCCCAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGaa                                              <  1:1240995/148‑1 (MQ=255)                                                                                                             aaaaacccccgcgcgGCGGCTCTCGGACGGTTTTTTTC‑CCCCCTAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGaa                                              <  1:739128/136‑1 (MQ=255)                                                                                                                      tccgggccggtTGTCTCCCATTTTTTTTC‑CCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                                    <  2:977659/141‑1 (MQ=255)                                                                                                                         gcggcgGCTGTCTCCCAGTTTTTTCC‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                 <  2:478079/148‑1 (MQ=255)                                                                                                                           ggctgcTGCCGGCCGTCTTTTTTT‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                                <  1:1653421/147‑1 (MQ=255)                                                                                                                           ggcggcGGTCTCCCGTTTTTTTTC‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                                <  2:2002130/147‑1 (MQ=255)                                                                                                                                 tGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGG‑AAAAAA‑ATTTAA‑ACCGCCCCGTTGTGTGGGGGGCGCTG‑C‑TGGGGACAAAGTGGGAAGGGGAGGggggcggggaaaccccacaa                        >  1:1465921/1‑135 (MQ=255)                                                                                                                                     tGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAG‑AATTTAACCCCTCATTTTGTTTGGGAACCGCTGCC‑GGGTACAAATTGGGAAAGGGAAGGGGGGGGAA‑AAACTCCAGaaaccccc                  >  1:1007340/1‑145 (MQ=255)                                                                                                                                       ccattttttttc‑CCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                <  2:1558171/144‑1 (MQ=255)                                                                                                                                       ccAGTTTTTTTC‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa                         <  2:1343536/142‑1 (MQ=255)                                                                                                                                         agtttttttc‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                <  2:2062321/144‑1 (MQ=255)                                                                                                                                           cTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGACGA‑G‑AAAAAG‑ATTTAAAACCCCCCTGTTTTGTGGGAAACGCTG‑C‑TGGGTACAAAGTGGGAAGGGAAGGGGGGGGGAACAGCACCCCAACCGCGCTaact             >  1:1813753/1‑148 (MQ=255)                                                                                                                                            ttttttc‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt              <  2:457810/146‑1 (MQ=255)                                                                                                                                             tttttt‑CCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa               <  1:1724815/146‑1 (MQ=255)                                                                                                                                             tttttcccccccgaAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa               <  1:550695/137‑1 (MQ=255)                                                                                                                                               tttc‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa           <  2:309110/148‑1 (MQ=255)                                                                                                                                               tttc‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg          <  1:1560770/149‑1 (MQ=255)                                                                                                                                               tttc‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt         <  2:806388/150‑1 (MQ=255)                                                                                                                                                      gccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCCAAA‑GAAAAAAA‑CTTTAA‑CCCCCCCTGTTGTGTGGGGAACCCTG‑CTTGGTACAAATG‑GGGAAGGGGAGGGGGGCGGATAAACTCCAGAACAGCGCCaaattagattgcc    >  2:848338/1‑138 (MQ=255)                                                                                                                                                           aacaaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑G‑AAAAAA‑CTTTAAAACCCCCCCGTTTTGTGGGAATCGCTG‑C‑TTGTAAAAAATGGGAAAAGGGAGGGGGGGGGATCAACTCCCGAAACCCGCTCATTGAGTTCGGCTc  >  2:1208560/1‑146 (MQ=255)                                                                                                                                                             |                                                                                                                                                      CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTAT‑CC‑AA‑CGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑C‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357740‑358033

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