Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,875	0	A	T	49.2%	54.1 / 50.7	67	N272K (AAT→AAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (25/9);  new base T (0/33);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 3.96e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357725‑358019                                                                                                                                                           |                                                                                                                                                     aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                                      <  4:156452/151‑1 (MQ=255)    gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                                   <  2:978248/151‑1 (MQ=255)    gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                                   >  2:1230817/1‑151 (MQ=255)     aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                                                  <  3:655467/151‑1 (MQ=255)     aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                                                  >  3:1143726/1‑151 (MQ=255)         cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGAGTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                               <  3:1123951/149‑1 (MQ=255)           tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccgcc                                                                                                                                           >  3:708152/1‑151 (MQ=255)               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGca                                                                                                                                        <  4:240954/150‑1 (MQ=255)                 tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                      >  4:76071/1‑150 (MQ=255)                 tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                      >  2:1372639/1‑150 (MQ=255)                  gATACGCGTTGCGCATCTGCCCATCCCCGTCCCGGTCCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTCTCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                              <  3:302769/141‑1 (MQ=14)                     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc                                                                                                                                 <  3:241710/151‑1 (MQ=255)                        cGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                              <  1:323628/151‑1 (MQ=255)                        cGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                              <  2:531534/151‑1 (MQ=255)                          tAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgccg                                                                                                                            >  2:567950/1‑151 (MQ=255)                            gcgcATCAGCCATTCCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCTTTCCCGCC‑AACAACCGCGCTTTTATCCAACGCAATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGTTTTTTTT‑CCGCCGCAAAAAAAccgccg                                                                                                                            <  4:219912/149‑1 (MQ=17)                                aTCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                                      >  2:1184139/1‑151 (MQ=255)                                       attcccACCCGCAGACCT‑GTGTCAGCGCCTACCGTGATTCTTCCCCCCA‑AACACCCGCGTTTTTACCCACCCCCATGTAAAAAAATG‑CT‑TCCGCGGCGGCTGTCGCCCATCTTTTTTC‑CCGCCACAACAAAACCGCCGTTCTCTCCagta                                                                                                                <  4:351636/146‑1 (MQ=255)                                       aaTTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  1:310924/1‑149 (MQ=255)                                       aaTTCCAGCCGAAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  4:239817/1‑149 (MQ=255)                                       aaTCCCACCCGCAGCCCT‑GTGCCGGCGGCTCCCGTGTTTCATTCCCGCC‑AACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATC‑CC‑TTGGCGGCGCCTGTCTCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                <  2:1192369/150‑1 (MQ=255)                                         ttCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                             >  2:1201444/1‑151 (MQ=255)                                           cccagcGCCGACCTGGTGTCAGCGGCTCCGGTGTTTCATTCCCGCC‑AACACCCGCGCTTTTATCCACCGCCTTGTGAAAAAATC‑CC‑TCCGCGGCGCCTGCCTGCCGGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                            <  4:60074/146‑1 (MQ=255)                                            cAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                           >  3:472123/1‑150 (MQ=255)                                               ccGCAGCCCT‑GTGTCAGCGGCTACCGTGTTTCTTTCCCGCC‑AACACCCGCGCTTTTGCCCAACGCCATGTGAAAAAATG‑CC‑TCCGCGGCGGCTGTCTGCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                         <  1:208563/148‑1 (MQ=18)                                               ccGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                       >  3:132699/1‑150 (MQ=255)                                                cGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                      >  2:89874/1‑150 (MQ=255)                                                   accccT‑GTGTCAGGGCCGCCGGTATTTCTTTTCCCCC‑ACCACCCGCGATTTTACCCACCGCCATGTG‑AAAAATGCCC‑TTCGCGGCGGTGGTCTGCAGGTTTTTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCaa                                                                                                   <  3:1255113/147‑1 (MQ=255)                                                   aGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggg                                                                                                     >  4:398861/1‑147 (MQ=255)                                                   aGACCT‑GTGTCAGCGGAGACCGGGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCACAAATG‑CC‑TTCGCAGCGGCTGTCTGCAAGCTATTTTT‑CCGCCGCAACAAAACCGCAGGTCTCTCCAGGAGTgg                                                                                                            >  4:75442/1‑142 (MQ=18)                                                     aCCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggg                                                                                                     >  2:512379/1‑145 (MQ=255)                                                      ccT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                      >  2:1032102/1‑144 (MQ=255)                                                      ccT‑GTGCCGGCGGCCCCCGTAATTATTTCCCGCA‑AACACCCGCTAATTCCCCCAACCCCATGTGAAAAAATG‑CC‑TCCGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                       <  2:544430/143‑1 (MQ=255)                                                        t‑gttTCGGCGGCTCCCGTATTTCATTCCCCCA‑AACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGG‑CC‑TTCGGGGCGGCTGTCTGCAGGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAgag‑a                                                                                              <  4:39771/149‑1 (MQ=255)                                                          gtgccggccgCTCCCGTTATTCTTTCCCCCA‑AACACCCGCCTTTTTCCCCACCCCCATGTAAAAAAATG‑CC‑TGGGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga                                                                                                 <  3:1277781/137‑1 (MQ=255)                                                            gtCGGCGCCTCCCGTTATTCTTTCCCCCA‑AACACCCGCGATTTTACCCACCGCCTTGTAAAAAAATC‑CT‑TCCGGGGCGTCTGTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAagag                                                                                                <  2:18861/145‑1 (MQ=255)                                                                 cGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGGGGC‑ggggggga                                                                                                  >  3:106887/1‑134 (MQ=255)                                                                  gtttcccGTGTTTCCTCCCCCCC‑AACACCCGCTCTTCTCCCCAACCCCATGAGAAAAAATG‑CC‑TGGGGGGGGGCTTCCTGCAAGTTTTTTTC‑CGGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑ttt                                                                                    <  3:766991/144‑1 (MQ=255)                                                                  ggCTACCGTGTTTCTTTCCCGCA‑AACACCCGCCCTTTTACCCACCCCCATGTGAAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAAt                                                                                           <  2:450386/143‑1 (MQ=255)                                                                     tACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCCAGGG‑aaaaaaattt                                                                                     >  2:1123874/1‑139 (MQ=255)                                                                       ccgtttttCTTTCCCCCC‑ACCACCCCCGCTTTTCCCCACCCCCATGTAAAAAAAGC‑CC‑TGGGCGGCGTTGGTCTGCAGGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑tt                                                                                     <  3:1250431/137‑1 (MQ=255)                                                                       ccGTGTTTCTTTCCCCCC‑AACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATG‑CC‑TTCGGGGCGGCTGTCTGCCAGCTTTTTTC‑CCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑aaa                                                                                            <  4:310921/137‑1 (MQ=255)                                                                             ttCTTTCCCGCC‑AAAACCCGCGCTTTTACCCACCGCCATGTGAAAAAATC‑CC‑TCCGGGGGGGCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtg                                                                        <  4:142717/150‑1 (MQ=255)                                                                             ttCATTCCCCCC‑AACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATG‑CC‑TTCGGGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAt                                                                         <  1:109989/149‑1 (MQ=255)                                                                                ctttCCCCCAAACATCGGCGCATTAATCCACCGCTATGGCCAAAAACG‑CT‑TGGGGGGGGACTGTATCCCCGCTTTTTTC‑CCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtg                                                                        <  1:131047/147‑1 (MQ=255)                                                                                 tttCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGG‑AAAAAGA‑TTTAAACCC‑TCAtgtt                                                                      >  4:202596/1‑148 (MQ=25)                                                                                    ccGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAGT‑TTTAATCCCG‑CCtgttgtt                                                                   >  2:1241322/1‑147 (MQ=25)                                                                                          aacaccCGCGTTTTTACCCAACCCCTTGTCAAAAAATC‑CC‑TCCGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATc                                                           <  3:89887/151‑1 (MQ=255)                                                                                             acccccccATTTACCCACCGCCGTG‑GCAAAAAAGCCC‑TTCGGGGGGGCGGTCGGCCTGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                         <  3:799873/143‑1 (MQ=255)                                                                                                     aTTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAG‑AAAAAAA‑TTTTAACCCCCCCTTTTTTGTGGGCATCgctggcgg                                                   >  2:556734/1‑144 (MQ=25)                                                                                                     aTTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCCAGAG‑AAAAAAA‑CTTAAACCCGGCCTGTTTTGTGGGGATCgctggt                                                     >  1:369503/1‑144 (MQ=37)                                                                                                      tttgcccAACGCTGTGGGAAAAATCG‑CT‑TCCGCGGCGGCGTTCCGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                 <  3:1076775/144‑1 (MQ=255)                                                                                                      tttaccaacccccATGTAAAAAAATG‑CC‑TCCGGGGCGGCTGTCGCCAAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                <  2:1307144/139‑1 (MQ=255)                                                                                                      tttACCCACCGCCAGGGAAAAAAATG‑CT‑TCCGGGGGGGTTGTCCGCCGGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                   <  3:853810/147‑1 (MQ=255)                                                                                                           ccACCGCCATGTGAAAAAATG‑CC‑TTCGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                           <  1:73279/150‑1 (MQ=255)                                                                                                                gCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGCGAG‑AAAAAAA‑ATTAAAAACCGCATGTTGTGTGGGGATCGCGGCTgggtaaaa                                             >  2:606128/1‑137 (MQ=21)                                                                                                                    tgtgAAAAAGTG‑CC‑TGCGCGGCGGCTTTCGGCCAGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                  <  4:376332/150‑1 (MQ=255)                                                                                                                       tgaaaaaTG‑CCTTTCGGGGCGGTTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                     <  1:398865/143‑1 (MQ=255)                                                                                                                        aaaaaaTC‑CC‑TCCGGGGGGCCTCTC‑GCCAGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                     <  3:440943/142‑1 (MQ=255)                                                                                                                          aaaaTG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAG‑AAAAAACATTTAAACCCCCCA‑CTTGTGTGGGGGTCGCTTCTGGTGGAAAAGTGGGAAGGGGAggggg                           >  3:412172/1‑147 (MQ=255)                                                                                                                              agcgcctcgggggcggcGGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                         <  2:226777/140‑1 (MQ=255)                                                                                                                                 cc‑TGCGCGGCGGCTGTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                        <  3:63941/148‑1 (MQ=255)                                                                                                                                       ggggcggcGGTCTGCAAGTTTTTTTC‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                   <  3:1241381/146‑1 (MQ=255)                                                                                                                                            ggttGTCTCCCAGTTTTTTTT‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg            <  1:280430/147‑1 (MQ=255)                                                                                                                                            ggtggTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa             <  3:1319981/145‑1 (MQ=255)                                                                                                                                            ggcTGTCTGCCAGTTATTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa             <  1:202598/149‑1 (MQ=255)                                                                                                                                                     ccAGCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc    <  2:1043078/149‑1 (MQ=255)                                                                                                                                                     ccAGCTTTTTTC‑CCCCCAAAAAAATCCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  2:472088/151‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                     ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357725‑358019

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: ‑