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

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035                                                                                                                                                             |                                                                                                                                                     cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGcaac                                                                                                                                                      >  4:43658/1‑151 (MQ=255)    aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAGGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGcaacaa                                                                                                                                                    <  4:468536/150‑1 (MQ=255)      aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACaaaa                                                                                                                                                  <  3:450645/150‑1 (MQ=255)      aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑c                                                                                                                                                >  1:1180187/1‑151 (MQ=255)       aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑cc                                                                                                                                               <  1:1316776/151‑1 (MQ=255)        gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑ccg                                                                                                                                              <  2:219326/151‑1 (MQ=255)          gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑ccgcc                                                                                                                                            >  1:717574/1‑151 (MQ=255)          gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑ccgcc                                                                                                                                            >  3:453719/1‑151 (MQ=255)           tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑ccgccg                                                                                                                                           >  2:111276/1‑151 (MQ=255)           tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑ccgcc                                                                                                                                            <  3:1385689/150‑1 (MQ=255)               gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTtct                                                                                                                                       >  1:1056865/1‑151 (MQ=255)                 aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTtctct                                                                                                                                     >  2:572437/1‑151 (MQ=255)                   cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                                   >  2:383205/1‑151 (MQ=255)                    aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                                   >  4:139719/1‑150 (MQ=255)                       cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                                >  3:440395/1‑150 (MQ=255)                           tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                           >  3:315960/1‑151 (MQ=255)                                ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                      >  4:480806/1‑151 (MQ=255)                                      cccTGTGGCGGCGTCTCCGGTAATTCATTCCCCCCAACCCCCCCGTTTTTACCCAC‑CGCAAGGT‑AAAAAATGCCT‑TCC‑GCGGCGTC‑TTTCTGCAAGCTTTTTT‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                <  3:44741/150‑1 (MQ=255)                                       ccTGTCTCAGCGCCTCCCGTATTTCTTTCCCCCCAACACCCGCGCATTTTCCCAA‑CCCTGTGT‑CAAAAAAGCCT‑TGC‑GGGGCGGC‑TTTCCACGAGCATTTTT‑GCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                                <  2:131663/150‑1 (MQ=255)                                        cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                      >  4:349995/1‑143 (MQ=255)                                         tgtgTCGGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAC‑CCCCATGTGAAAAAATCCCT‑T‑C‑GGGGCGGC‑TGTCTCCCAGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                                   <  3:693187/145‑1 (MQ=255)                                          gtgccggCGGCTCCCGTGATTCATTCCCCCAAACACCCGCCCTTTTACCCAC‑CGCCATGGAAAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCATTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                 <  3:430991/140‑1 (MQ=255)                                           tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                                    >  1:980659/1‑140 (MQ=255)                                            gtCAGCGGCTCCCGTAATTCATTCCCGCAAACACCCGCGCTTTTATCCAA‑CGCCATGTGAAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                            <  2:110276/149‑1 (MQ=255)                                                  ggCTACCGTTATTTTTTCCCGCAAACACCCGCTCATTTCTCCAC‑CCCAATGGGAAAAATTGCTG‑T‑G‑GGGGCGGC‑TGTCTGCGACTTTTTTC‑CCCACCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTaa                                                                                                   <  1:1181241/151‑1 (MQ=255)                                                   gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCCGTCGTGGCgggggagcggataaa                                                                                                           >  3:1159782/1‑132 (MQ=21)                                                     tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGgaaaaagttt                                                                                                      >  3:1115756/1‑139 (MQ=25)                                                      aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGaaaaaaaatt                                                                                                      >  3:498272/1‑138 (MQ=255)                                                       ccGTGTTTCTTTCCCCCCAACACCCGC‑CATTTTCCCAACCGCCATGTAAAAAAATCCCT‑G‑C‑GGGGCGGT‑TGTCTCCCAGTTTTTTT‑CCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑tt                                                                                                      <  1:521132/143‑1 (MQ=255)                                                         gttaTTCTTTCCCGCAAACACCCGCGCTTTTACCCAC‑CCCCATGTAAAAAAATGCCT‑G‑C‑GCGGCGGC‑TGTCTGCCAGTTTTTTT‑TCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGt                                                                                             <  3:452620/147‑1 (MQ=255)                                                              tCATTCCCGCCAACACCCCCGCTTTTCCCAAC‑CGCTTTGTGCAAAAGT‑CCC‑T‑G‑GGGGCGGCTTGTCCG‑CATGTTTTTT‑TCCCCCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTc                                                                                            <  3:1350203/146‑1 (MQ=255)                                                                 ttcccGCCAACACCCGCGCTTTTACCAAC‑CCCCATGT‑AAAAAATGGCT‑TCC‑GGGGCGGC‑TTCCTGCCGGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa                                                                                           <  2:135542/141‑1 (MQ=255)                                                                       caaacaccCGCGCATTTCTCAAG‑CCCCATGTAAAAAAATCCCT‑C‑C‑GCGGGGGC‑TGCCGGCCGGCTTTTTT‑CCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGc                                                                                <  4:335533/147‑1 (MQ=255)                                                                           cacccGCACATTCACCAAA‑CGCAATGTAAAAAAATGCCT‑G‑G‑GGGGCGGG‑TGCCTCCGAGTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                           <  1:734785/147‑1 (MQ=255)                                                                               cgcgCTTTTACCCAC‑CCCCATGGCAAAAAATGCCT‑C‑C‑GGGGCGGC‑TGTCTCCCAGTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                      <  3:730072/151‑1 (MQ=255)                                                                                     tttACCCAC‑CCCCATGTCAAAAAATCCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGTTATTTT‑CCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTa                                                                 <  3:1353339/149‑1 (MQ=255)                                                                                          ccAC‑CGCCATGTGAAAAAATGCCT‑C‑C‑GCGGGGGT‑TGTCTGCCAGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACaa                                                             <  2:93686/148‑1 (MQ=255)                                                                                           caaccaccATGTGCAAAAATCCCT‑T‑G‑GGGGCGTC‑GTCCGGCGAGTTATTTT‑CCCCCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAt                                                            <  4:177859/143‑1 (MQ=255)                                                                                             a‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAAAAAA‑TTTTATACCCGCCTTTTGTGTGGGGATCGCCGCT‑ggggaaaaagggg                                                        >  4:523705/1‑139 (MQ=19)                                                                                              accgctatgggaaaaaTGCCT‑T‑CGGGGGCGCC‑GTTCCGCCATTTTTTTG‑CCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtgtg                                                         <  1:1242790/139‑1 (MQ=255)                                                                                               ccctattagaaaaaaTGCCT‑G‑G‑GCGGCGGT‑TCCCTGCCCGTAATTTT‑CCCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaa                                                      <  2:43659/141‑1 (MQ=255)                                                                                                 ccgggtcaaaaaaTCCTT‑T‑C‑GGGGCGGC‑TTTCGGCCAGTTTTTTC‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAt                                                            <  1:1201826/135‑1 (MQ=255)                                                                                                 ccATGTAAAAAAATGCCT‑G‑C‑GGGGGGTC‑TGTCTGCCAGTTTTTTT‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg                                                    <  1:586258/151‑1 (MQ=255)                                                                                                   agtgacAAAAAACCCT‑T‑CGGGGGCGGC‑TTTCTGCCAGTTTTTTT‑CCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg                                                    <  1:1369467/145‑1 (MQ=255)                                                                                                     gtaaaaaaaTGCCT‑C‑C‑GGGGGGGT‑TGCCTGCAAGTTTTTTT‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAgg                                                   <  1:1235781/144‑1 (MQ=255)                                                                                                      ggCAAAAAAGCCTCC‑C‑GCGGGGCC‑TGTCGGCAGGTTATTTT‑CCCGCCAAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAt                                                            <  1:61135/138‑1 (MQ=255)                                                                                                         aaaaaGCCCT‑C‑C‑GGGGGGGT‑TGTCGGCCGTCTTTTTT‑CCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGaa                                                <  1:532283/147‑1 (MQ=255)                                                                                                          aaagTGCCT‑C‑C‑GGGGGGTC‑TGTCTGCCAGTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTgcg                                           <  1:1278485/151‑1 (MQ=255)                                                                                                          aaaaTGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAAA‑TTTAAAAACACCCTGTTGTGTGGGGATCAGTGCTGGGGGAAAATGTgg‑aagggaaggggg                                           >  4:198256/1‑147 (MQ=255)                                                                                                                 ccccgc‑gcgGCGGC‑TGTCT‑CCGGTTTTTTTCCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                                       <  1:1135738/145‑1 (MQ=255)                                                                                                                  t‑t‑C‑GCGGCGCC‑TGTCGGCAAGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                                       <  1:291809/147‑1 (MQ=255)                                                                                                                                 tgtcggcaggttttttt‑cCCGCCGCAACAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                            <  4:491137/134‑1 (MQ=255)                                                                                                                                       ccgtttttttt‑cCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                          <  1:368445/136‑1 (MQ=255)                                                                                                                                        cAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGAGAAAAAAA‑ATTATAACCGCCCTGTTGTTTGGGAATGGCCGCT‑TGTGACAAAGTGGGAAGGGGAGAGGGGGGGGACAGATCCAAAACAgcgcca                  >  1:1288505/1‑148 (MQ=255)                                                                                                                                            tttttt‑cCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg              <  3:61610/148‑1 (MQ=255)                                                                                                                                                     cgccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACATTTAAACCC‑TCCTGTTTTGTGGGGAACGCTGCT‑GGTGAAAAATTGGGAAGGGGAGGGGGGGGGAAAAGCTCCCGAAACCCCCCAATTGAGtttggc      >  1:1231584/1‑146 (MQ=255)                                                                                                                                                     cgccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAAAAAA‑TTTTAAACCATCATGGTTTTTGGGGATACCTTCT‑GGTGACAAAGTGGGAAGGGAAAGGGGGGGAATCAACCCCAGAACCCCGCCAATTTAGTTTGcc      >  3:78708/1‑150 (MQ=255)                                                                                                                                                          caacaaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGGGATAGAAAAAAA‑GTTAAAGCCCTCCTCCTGTGCGGGAAACGCTGGT‑GGGAAAAATATGGGAAGGGGAGGGCGGGGCGTCACCTCCAAAACCCCCCTAAATGaaatagccccca  >  2:60791/1‑138 (MQ=255)                                                                                                                                                          caacaaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAAA‑ATTAAACCCCTCCTTTTTTGTGGGCAACGCTGCT‑GGGTACAAAGTGGGAAGGGGAGGGGGGGGGAACAGCCCCACAACCCCCCCAATTTTGTTCGCCTc    >  1:564350/1‑147 (MQ=255)                                                                                                                                                             |                                                                                                                                                     CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035

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