Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	49.2%	59.0 / 26.9	63	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (23/9);  new base C (0/31);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 3.82e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357733‑358028                                                                                                                                                        |                                                                                                                                                         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                                                                                                                                       <  2:531534/151‑1 (MQ=255)                cGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                                       <  1:323628/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)                                                                                                                                             caagtttttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                      <  1:58961/133‑1 (MQ=255)                                                                                                                                                 tttttttc‑CCGCCGAAACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt        <  3:1372708/147‑1 (MQ=255)                                                                                                                                                 tttttttc‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt         <  2:1143671/146‑1 (MQ=255)                                                                                                                                                    ttttc‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt        <  3:424858/147‑1 (MQ=255)                                                                                                                                                      ttt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAACCCGTCATGTTGTGTGGCAATCGCTGCTGGGAACAATGTGGGAAGGGGAGAGCGGCGAATCAGCTCCAGAAACGCGCTAATTGAGttt  >  1:80278/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                         CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCT‑GTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CC‑TTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357733‑358028

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