Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	42.1%	82.9 / 22.2	58	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (22/11);  new base C (0/24);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 5.71e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357732‑358023                                                                                                                                                           |                                                                                                                                                  tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑Tcc                                                                                                                                                 >  1:1209959/1‑151 (MQ=255)  cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑Tccg                                                                                                                                                <  2:5390/151‑1 (MQ=255)   gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑Tccgc                                                                                                                                               >  3:1181377/1‑151 (MQ=255)     aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑Tccgccg                                                                                                                                             <  3:231892/151‑1 (MQ=255)       gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGca                                                                                                                                           <  4:113480/151‑1 (MQ=255)        ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGcaa                                                                                                                                          <  4:518925/151‑1 (MQ=255)         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                                                                                                                               >  3:453719/1‑151 (MQ=255)                  gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑ccgcc                                                                                                                               >  1:717574/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‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                       <  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‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                          <  1:521132/143‑1 (MQ=255)                                                                 gttaTTCTTTCCCGCAAACACCCGCGCTTTTACCCAC‑CCCCATGTAAAAAAATGCCT‑G‑C‑GCGGCGGC‑TGTCTGCCAGTTTTTTT‑TCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                 <  3:452620/147‑1 (MQ=255)                                                                      tCATTCCCGCCAACACCCCCGCTTTTCCCAAC‑CGCTTTGTGCAAAAGT‑CCC‑T‑G‑GGGGCGGCTTGTCCG‑CATGTTTTTT‑TCCCCCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                                <  3:1350203/146‑1 (MQ=255)                                                                         ttcccGCCAACACCCGCGCTTTTACCAAC‑CCCCATGT‑AAAAAATGGCT‑TCC‑GGGGCGGC‑TTCCTGCCGGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                               <  2:135542/141‑1 (MQ=255)                                                                               caaacaccCGCGCATTTCTCAAG‑CCCCATGTAAAAAAATCCCT‑C‑C‑GCGGGGGC‑TGCCGGCCGGCTTTTTT‑CCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                    <  4:335533/147‑1 (MQ=255)                                                                                   cacccGCACATTCACCAAA‑CGCAATGTAAAAAAATGCCT‑G‑G‑GGGGCGGG‑TGCCTCCGAGTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                               <  1:734785/147‑1 (MQ=255)                                                                                       cgcgCTTTTACCCAC‑CCCCATGGCAAAAAATGCCT‑C‑C‑GGGGCGGC‑TGTCTCCCAGTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                          <  3:730072/151‑1 (MQ=255)                                                                                             tttACCCAC‑CCCCATGTCAAAAAATCCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGTTATTTT‑CCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTa                                                     <  3:1353339/149‑1 (MQ=255)                                                                                                  ccAC‑CGCCATGTGAAAAAATGCCT‑C‑C‑GCGGGGGT‑TGTCTGCCAGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACaa                                                 <  2:93686/148‑1 (MQ=255)                                                                                                   caaccaccATGTGCAAAAATCCCT‑T‑G‑GGGGCGTC‑GTCCGGCGAGTTATTTT‑CCCCCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAt                                                <  4:177859/143‑1 (MQ=255)                                                                                                     a‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAAAAAATTTTATACCCGCCTTTTGTGTGGGGATCGCCGCT‑ggggaaaaagggg                                            >  4:523705/1‑139 (MQ=19)                                                                                                      accgctatgggaaaaaTGCCT‑T‑CGGGGGCGCC‑GTTCCGCCATTTTTTTG‑CCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtgtg                                             <  1:1242790/139‑1 (MQ=255)                                                                                                       ccctattagaaaaaaTGCCT‑G‑G‑GCGGCGGT‑TCCCTGCCCGTAATTTT‑CCCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaa                                          <  2:43659/141‑1 (MQ=255)                                                                                                         ccgggtcaaaaaaTCCTT‑T‑C‑GGGGCGGC‑TTTCGGCCAGTTTTTTC‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAt                                                <  1:1201826/135‑1 (MQ=255)                                                                                                         ccATGTAAAAAAATGCCT‑G‑C‑GGGGGGTC‑TGTCTGCCAGTTTTTTT‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg                                        <  1:586258/151‑1 (MQ=255)                                                                                                           agtgacAAAAAACCCT‑T‑CGGGGGCGGC‑TTTCTGCCAGTTTTTTT‑CCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg                                        <  1:1369467/145‑1 (MQ=255)                                                                                                             gtaaaaaaaTGCCT‑C‑C‑GGGGGGGT‑TGCCTGCAAGTTTTTTT‑CCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAgg                                       <  1:1235781/144‑1 (MQ=255)                                                                                                              ggCAAAAAAGCCTCC‑C‑GCGGGGCC‑TGTCGGCAGGTTATTTT‑CCCGCCAAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAt                                                <  1:61135/138‑1 (MQ=255)                                                                                                                 aaaaaGCCCT‑C‑C‑GGGGGGGT‑TGTCGGCCGTCTTTTTT‑CCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGaa                                    <  1:532283/147‑1 (MQ=255)                                                                                                                  aaagTGCCT‑C‑C‑GGGGGGTC‑TGTCTGCCAGTTTTTTT‑CCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTgcg                               <  1:1278485/151‑1 (MQ=255)                                                                                                                  aaaaTGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAAATTTAAAAACACCCTGTTGTGTGGGGATCAGTGCTGGGGGAAAATGTgg‑aagggaaggggg                               >  4:198256/1‑147 (MQ=255)                                                                                                                         ccccgc‑gcgGCGGC‑TGTCT‑CCGGTTTTTTTCCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                           <  1:1135738/145‑1 (MQ=255)                                                                                                                          t‑t‑C‑GCGGCGCC‑TGTCGGCAAGTTTTTTT‑CCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                           <  1:291809/147‑1 (MQ=255)                                                                                                                                         tgtcggcaggttttttt‑cCCGCCGCAACAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                <  4:491137/134‑1 (MQ=255)                                                                                                                                               ccgtttttttt‑cCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa              <  1:368445/136‑1 (MQ=255)                                                                                                                                                cAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGAGAAAAAAAATTATAACCGCCCTGTTGTTTGGGAATGGCCGCT‑TGTGACAAAGTGGGAAGGGGAGAGGGGGGGGACAGATCCAAAACAgcgcca      >  1:1288505/1‑148 (MQ=255)                                                                                                                                                    tttttt‑cCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg  <  3:61610/148‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                  TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCT‑T‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTT‑TCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357732‑358023

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