Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	50.8%	55.3 / 103.2	65	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (24/8);  new base C (0/33);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 2.65e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-11

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357731‑358027                                                                                                                                                            |                                                                                                                                                        cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                         >  2:23044/1‑150 (MQ=255) cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                         >  1:1078073/1‑150 (MQ=255) cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                                       <  2:351079/151‑1 (MQ=255) cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                                       <  2:916100/151‑1 (MQ=255) cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                                       <  1:741358/151‑1 (MQ=255)  tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                      <  2:1202813/151‑1 (MQ=255)  tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                      <  2:652731/151‑1 (MQ=255)       ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑ccgccg                                                                                                                                                  >  1:189711/1‑150 (MQ=255)       ggCCCGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑ccgccg                                                                                                                                                  >  2:491317/1‑150 (MQ=255)          cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaa                                                                                                                                               >  1:702184/1‑150 (MQ=255)          cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaa                                                                                                                                               >  1:535202/1‑150 (MQ=255)            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CCACCGCGCATTTATCCAACGCCA‑TGTGCCAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                            >  1:571782/1‑151 (MQ=255)            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                            >  2:1644016/1‑151 (MQ=255)            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                            >  1:439492/1‑151 (MQ=255)                 gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg                                                                                                                                       >  2:517849/1‑151 (MQ=255)                 gCGTAGCGCACCACGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGGGATT‑CATTTCCGCCAA‑CAACCCCCCATGTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCGGTCTGGCTGCTATTGTT‑ACTCCGCGACAAAAccg                                                                                                                                       >  1:1690266/1‑151 (MQ=17)                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgccg                                                                                                                                    >  2:1279743/1‑151 (MQ=255)                             aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                           >  2:1161824/1‑151 (MQ=255)                              ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                                                          >  1:1341210/1‑151 (MQ=255)                               gaatttCCCGCCGCACCCCTGGGGCGGCGGCTCCCGTGTTT‑ATTTTCCCCCAA‑CACCCGCGTTTTTATCCAACGCCT‑TGTGAAAAAA‑T‑GCCTCCGG‑GGCGGCTTTCTCCCAGCTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTcc                                                                                                                            <  2:1207122/144‑1 (MQ=255)                               cattccccACCCGCAGACTTTTGTCAGCGGCTCCGGTATTT‑CTTTCCCCCCAA‑CACCCGCGCTTTTACCCACCCCAT‑TGTCAAAAAA‑T‑CCCTCCGC‑GGCGCCTTTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCagt                                                                                                                         <  1:1192672/145‑1 (MQ=255)                                  aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                      >  1:1649749/1‑151 (MQ=255)                                    cccccccGCGGCCCGGTGCCAGCGCCCCCCGTTTTT‑CTTTCCCCCCAA‑CACCCGCGTTTTTACCCACCCCCA‑TGTGAAAAAA‑T‑GCCTCCGG‑GGCGGCTTTCCCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    <  2:1457282/146‑1 (MQ=255)                                      cacccGCAGCCCTGTGCCAGGGGCTCCGGTTATT‑CTTTCCCCCAAA‑CACCCGCCCTTTTACCCACCCCCA‑TGGGAAAAAA‑T‑CCTTCCGC‑GGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                 <  2:1475505/148‑1 (MQ=255)                                        gCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑g                                                                                                              >  1:1013436/1‑151 (MQ=255)                                              gACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCCACAAAACCGCCGTTCTCTCCAGTAGGGGC‑gg‑ggggaga                                                                                                        >  2:1085092/1‑146 (MQ=255)                                               ccgtttgtCGGCGCCTCCCTTGTTT‑CTTTCCCCCCAA‑CAACCGCACTTTAACCCACCGCCA‑TGTAAAAAAA‑T‑GCTTTCGG‑GGGGCCTGTCTGCAATCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑agag                                                                                                      <  2:94825/146‑1 (MQ=255)                                               aCCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGGC‑agag                                                                                                      >  2:1382828/1‑151 (MQ=255)                                                   gtgtCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑ggggcgggaaaa                                                                                                   >  1:238201/1‑141 (MQ=255)                                                   gtgtCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑ggggcgaggaaa                                                                                                   >  1:898773/1‑141 (MQ=255)                                                    tgtCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑ggggcgaagaaaa                                                                                                  >  2:1052313/1‑140 (MQ=255)                                                     tgtcggcgcctccctGATTCCATTCCCGCCAC‑CACCCGCGATTTTACCCACCCCCT‑TGTGAAAAAA‑G‑CCCTTCGC‑GGCGGTTGTCTCCCGGCTTTTTGC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAAt                                                                                                  <  2:916136/136‑1 (MQ=255)                                                      gcAGCGGCTCCGGTGTTT‑C‑TTTCCCCCAACCACCCGCGCTTTTACCCACCGCCA‑TGTAAAAAAT‑T‑GCCTCCGC‑GGGGGCGGTCGCCCATTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAg                                                                                                <  2:934236/149‑1 (MQ=255)                                                      ccggCGGCTACCGGTATT‑CTTTCCCCCCAA‑CACCCGCGCTTTTACCAACCGCAA‑GGTAAAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGc                                                                                               <  2:857561/148‑1 (MQ=255)                                                          cGGCTCCGGTAATT‑CTTTCCCCCAAA‑CACCCGCGCATTTACCCAACCCCA‑GGTGAAAAAA‑T‑GCCTCCGG‑GGGGGCGGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCtt                                                                                             <  2:189769/149‑1 (MQ=255)                                                             ctcccGTGTTT‑CTTTCCCCCAAA‑CACCCGCGCTTTTACCCACCGCCA‑TGTGAAAAAA‑T‑GCCTCCGG‑GGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGc                                                                                        <  2:1003489/148‑1 (MQ=255)                                                              tccggtgtTT‑CTTTCCCCCCAA‑C‑ACCCCCCATTTTACCCACCGCA‑AGTGGAAAAA‑TGGCCTCGGG‑GGGGGCTGTCGCCCAGTTTTTTTC‑CCCCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGcc                                                                                       <  2:595868/147‑1 (MQ=255)                                                              caCCGTGATT‑CATTCCCCCCAA‑CACCCGC‑CTTTTTCCCCACCCCAGTGTGAAAAAA‑T‑CCCTTCGC‑GGCGGTTTTCTCCCGGTTATTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGcc                                                                                       <  1:489991/150‑1 (MQ=255)                                                                ccGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGG‑GGCA‑AGGGAAAAAATTTTAAAcc                                                                                       >  2:17377/1‑150 (MQ=25)                                                                ccGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGGG‑AGAGgaaaaaagtttaaacc                                                                                      >  2:879901/1‑137 (MQ=25)                                                                    taTT‑CTTTCCCGCCAA‑CAACCGCCCTTTTACCCACCGCCA‑TTTGCAAAAA‑T‑CCCTCCGC‑GGCGCTTGTCGGCCGGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAt                                                                                 <  1:241940/150‑1 (MQ=255)                                                                      tt‑CTTTCCCCCCAA‑CACCCGCCCTTTTACCCACCGCAA‑TGTGAAAAAA‑G‑CCCTTCGG‑GGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgt                                                                               <  2:1482349/151‑1 (MQ=255)                                                                         ctttccccccAA‑CACCCGCCATTTTACCCAACGCAA‑TGTGAAAAAA‑T‑CCCTCCGC‑GGCGGTTGTCTGCCAGTTATTTTC‑CCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1487730/143‑1 (MQ=255)                                                                         cATTCCCCCCAA‑CACCCGCGCATTTACCCAACGCCA‑TGTGAAAAAA‑T‑GCCTCCGC‑GGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1197807/151‑1 (MQ=255)                                                                         cATTCCCCCCAA‑CACCCCCCCTTTTACCCAACCCCA‑TGTGAAAAAT‑T‑GCCTTCGG‑GGGGGCTGTCGGCCAGTTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1425841/151‑1 (MQ=255)                                                                              cccccac‑caCCCGCGTATTAACCCAACGCCA‑‑GGGCAAAAATT‑CCCTCCGC‑GGGGGTTGTCGGCCAGTTTTTTCC‑CCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  1:24692/148‑1 (MQ=255)                                                                              ccccaaa‑caccCGCGCTTTTCCCCACCCCCG‑TGTACAAAAA‑T‑CCCTCCGC‑GGGGGCTGTCGGCCAGTTATTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  2:1401319/146‑1 (MQ=255)                                                                                        acCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGCGGCAGAGAGAAA‑AGTCTCAAGCCCCTCAGTTTTTGTGGCacccccg                                                                >  1:71896/1‑145 (MQ=21)                                                                                              cATTTTTCCACCCCCT‑TGTCCAAAAT‑T‑CCCTCCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                         <  2:1217167/151‑1 (MQ=255)                                                                                              cATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGG‑GGCA‑AAGGAAAAAACTTTAACCCCCCCCCGTTGT‑TGGGAAGCgctgcgcggg                                                         >  1:1380470/1‑146 (MQ=25)                                                                                                tttACCCACCGCAA‑TGTGAAAAAA‑T‑GCTTCCGG‑GGGGCCTGTCGGCCGTCTATTTTC‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  1:505403/151‑1 (MQ=255)                                                                                                tttACCCAACCCCA‑TGTAAAAAAA‑T‑CCCTTCGG‑GGGGGCTGTC‑GCCAGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  2:1205507/151‑1 (MQ=255)                                                                                                ttaacccACCGCCG‑GGTGAAAAAA‑T‑CCTTTCGC‑GGGGGCTGCCCCACAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  2:639640/146‑1 (MQ=255)                                                                                                    ctcaaccgctgggtgAAAAAA‑T‑GCCTGCGG‑GGCGGCTTTCGGCCAGTTATTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                    <  1:626670/139‑1 (MQ=255)                                                                                                    cccACCGCCA‑TGTGAAAAAA‑T‑CCTTTCGG‑GGCGGCTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                   <  2:1063106/150‑1 (MQ=255)                                                                                                        aCGCCA‑GGGGAAAAAA‑G‑CCCTTCGG‑GGGGGCTTTCGCCAAGTTTTTTTC‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                               <  2:1299161/151‑1 (MQ=255)                                                                                                         cccccttgtaaaaaaa‑T‑CCCTTCGG‑GGGGGCTGTCTGCCATTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                               <  2:1654513/140‑1 (MQ=255)                                                                                                         cGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCC‑AGAGAAAAAATTTTTAACCCCCCATGTTTTTTGGGCAACGCTGCTGGtgaaaaagaggga                                              >  2:887553/1‑139 (MQ=25)                                                                                                         acccccatttaaaaaa‑a‑TCCTTCGGGGGCGGTTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                <  2:1564697/140‑1 (MQ=255)                                                                                                               tgtgAAAAAA‑T‑CCCTCCGG‑GGGGGCTGCCGCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  2:365101/151‑1 (MQ=255)                                                                                                               tgaaaaaaaa‑T‑CCCTCCGC‑GGCGCCTTTCTGCCAGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  1:1161254/146‑1 (MQ=255)                                                                                                                                          ttGTCTGCC‑GCTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:684132/150‑1 (MQ=255)                                                                                                                                          cttccTGCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:637447/147‑1 (MQ=255)                                                                                                                                          cggccggCCAGTTTTTTTC‑CCGCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:517635/145‑1 (MQ=255)                                                                                                                                          cTGTCTGCCGGTTTTTTTC‑CCGCCACAAAAAGACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:1555412/151‑1 (MQ=255)                                                                                                                                          cTGTCGGCCAGCTTTTTTC‑CCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  2:1352039/151‑1 (MQ=255)                                                                                                                                          cTGTCAGCCAGTTGTTTCC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:246352/151‑1 (MQ=255)                                                                                                                                                       tttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg    <  2:393153/150‑1 (MQ=255)                                                                                                                                                         tttc‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:439645/151‑1 (MQ=255)                                                                                                                                                            |                                                                                                                                                        CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357731‑358027

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