Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA49.2% 27.4 / 42.5 63L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base A (0/31);  total (18/45)
Fisher's exact test for biased strand distribution p-value = 2.62e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357740‑358038
                                                                                                                                                          |                                                                                                                                                       
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑caac                                                                                                                                                         >  2:894672/1‑151 (MQ=255)
   aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACaaa                                                                                                                                                      <  2:1080937/151‑1 (MQ=255)
      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTTCCAGCTATTTTTCCGCCG‑CAACAAAAcc                                                                                                                                                   <  2:1093957/151‑1 (MQ=255)
      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAcc                                                                                                                                                   <  2:725644/151‑1 (MQ=255)
        cGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAccgc                                                                                                                                                 <  2:497664/151‑1 (MQ=255)
            gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGtt                                                                                                                                             >  2:831014/1‑151 (MQ=255)
               cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctc                                                                                                                                          >  2:1059713/1‑151 (MQ=255)
                aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctct                                                                                                                                         >  2:1248067/1‑151 (MQ=255)
                 tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctctc                                                                                                                                        >  1:726310/1‑151 (MQ=255)
                  cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGCGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCCATTTTTCCGCCG‑CAACCAAACCGCCGTTCTCTcc                                                                                                                                       >  2:747198/1‑151 (MQ=255)
                   aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTcc                                                                                                                                       >  1:579912/1‑150 (MQ=255)
                       aaTCCCAGCCGCAGCCTTGTGTCAGCGCCTCCCTTGTTTCATTCCCGCCAACACCCCCGCTTTTATCAAACGCCA‑TGGGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTATTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                   <  2:694596/150‑1 (MQ=255)
                           ccAGCCGCAGCCCTGTGTCGGCGGCTCCCGTTTTTCTTTCCCGCAAACACCCGCGCTTTTACCCACCGCCA‑TGTAAAAAAATCCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                                <  1:1303261/149‑1 (MQ=255)
                           ccACCCCCAGCCTTGTGTCAGCGCCTCCCTTGATCCATTCCCCCCAACCCCCGCCCTTTTACCCACCCCCA‑TGTAAAAAAATGCTTC‑C‑GCGGCGCTTTTCTGCAAGTTTTTTTCCCCCCG‑AAACAAAACCGCCGTTCTCTCCagta                                                                                                                                   <  2:534347/146‑1 (MQ=255)
                            cccccGCAGCCTTGTTTCCGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGTTTTTACCCACCCCCA‑TGCACAAAAATCCTGC‑C‑GCGGCGGTTGTCTGCAAGTTTTTTTTCCGCCG‑CAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                   <  2:479296/142‑1 (MQ=255)
                            cAGCTGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑GGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                             >  2:782552/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                             >  1:892113/1‑151 (MQ=255)
                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                            >  1:283682/1‑151 (MQ=255)
                               ccGCAGCCTTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCATTTACCCACCGCCA‑TGTGCAAAAATGCCTC‑C‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                          <  1:433027/151‑1 (MQ=255)
                                cGCAGCCTTGTGTCGGCGGCCCCCGTAATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCCA‑TGTCAAAAAGTCCCTT‑C‑GCGCGGGCTGTCTCCCAGTTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                                <  2:772590/144‑1 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                         >  2:691697/1‑150 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                           >  1:415606/1‑149 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                            >  2:1226613/1‑148 (MQ=255)
                                  cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                            >  2:1063172/1‑146 (MQ=255)
                                    ggccTGTGTCGGCGCTTGCCGTGATTCTTTCCGGCAAACACCCGCGCATTTCTCCAACGCAA‑TGT‑CAAAAAAGCCTTCC‑GGGGCGGTTGTTGGCCGGTTTTTTTCCCGCCC‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                     <  1:966288/149‑1 (MQ=255)
                                      cttgtgTCAGCGGCTCCCGTGATTCTTTCCCCCAAACACCCGCGCTTTTTCCCAACTCCA‑TTCAAAAAATTGCTTC‑C‑GCGGCGCCTCCCTCCCCGTTTTTTCCCCCCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g                                                                                                                   <  1:895282/148‑1 (MQ=255)
                                        tgtgTCAGCGCCTCCCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCA‑TGTGAAAAAATCCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g                                                                                                                   <  1:916213/148‑1 (MQ=255)
                                                cgcctcccgtttttCATGCCCCCCAACACCCGCCCTTTTCCCCACCGCCA‑TGAGAAAAAATGCCTC‑C‑GCGGGGGCTGTCTGCCGGTTATTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATa                                                                                                             <  2:356490/133‑1 (MQ=255)
                                                 ggCTCCCGTGTTTCTTTCCCCCCAAAACCCGCGCTTTTATCCAACGCCT‑TGTGAAAAAATGCTTT‑C‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCttt                                                                                                        <  1:1197288/149‑1 (MQ=255)
                                                            ttCATTCCCCCCAACACCCCCGCATTTCCCCACCGCCA‑TGTGCAAAAATGCTTT‑C‑GCGCCGCCTTTCTCCCAGCTTTTTTCCCCCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCa                                                                                              <  1:417763/148‑1 (MQ=255)
                                                              cttgcccGCCAACACCCGCGCATTCCCCCACCGCCG‑GGGGCAAAAAGGCCGC‑G‑GGGGGGGCGGGCGGCGGGTTATTTTCCCGCCG‑AAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCg                                                                                                 <  2:78909/138‑1 (MQ=255)
                                                              cttccccGCCAACACCCGCGCTTTAATCCACCGCCA‑TGTG‑AAAAATGCCCT‑CGGCGGCGCTTGTCTCCCAGCTATTTTCCCGCCA‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCa                                                                                              <  2:1103724/141‑1 (MQ=255)
                                                                ttccccccaccaCCCGCGCTTTTA‑CCAACCCCATTGTAAAAAAATGCCTT‑C‑GCGCGGCCTGTCTGCCGTTTATTTGCCCGCCA‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgt                                                                                           <  2:824138/141‑1 (MQ=255)
                                                                ttccccccAACACCCCCGCTTTTACCCACCCCAT‑TGTGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                                        <  1:109085/144‑1 (MQ=255)
                                                                 ttCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑AAAAAACTTTTAGCCCCTCCtgtt                                                                                          >  1:37287/1‑147 (MQ=25)
                                                                  tCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGAGG‑AAAAAACTTTaaaaccctaccttttgt                                                                                       >  2:601505/1‑135 (MQ=21)
                                                                  ccccccAACACCCGCTCTTTCACCAAACGCCA‑TGTCCAAAAATCCCTC‑C‑GGGGCGCTTTCCGGCCAGTTTTTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttg                                                                                         <  2:39567/143‑1 (MQ=255)
                                                                      ccaccaccccctcATTTTCCCACCCCCG‑TGTCCAAAAAGCCCTG‑C‑GCGGCGGCTGTCGGCCAGCTATTTTCCCGCCA‑AAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                                 <  2:164420/139‑1 (MQ=255)
                                                                      ccAACACCCGCGTATTTACCCACCGCCA‑TGTGAAAAATTGCCTT‑C‑GGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                                   <  1:862678/149‑1 (MQ=255)
                                                                        aacaccCGCGCTTTTATCCACCGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGTTATTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                               <  1:54560/151‑1 (MQ=38)
                                                                         acaacCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAACTTTAA‑CCCCTCATGTTTTGTGGGaa                                                                                 >  1:1272931/1‑148 (MQ=25)
                                                                          cacccccgattttccccAACGCCA‑GGTGAAAAAATGCTTT‑C‑GGGGGGGTGGTCTGCCAGTTATTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                              <  2:891505/135‑1 (MQ=255)
                                                                                   aTTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACACAACCGCCGTTCTCTCCAGTAGTGGCGGCGCACGC‑G‑AAATATCATTTAAACCGTCTGGTGGTGTGGGCGTCACTGCTgg                                                                       >  2:770576/1‑148 (MQ=21)
                                                                                    ttcccccAACGCCA‑TGTGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCGGCTTTTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                                    <  2:2457/145‑1 (MQ=255)
                                                                                            acccccttgacaaaaaaTGCCTC‑G‑GCGGCGCCTGTCGGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTgg                                                           <  2:443190/139‑1 (MQ=255)
                                                                                               ccA‑TGTAAAAAAACCCCTT‑C‑GCGGCGGTTGTCTGCCATCTTTTTCCCCGCCG‑AAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑aa                                                                <  2:1022275/142‑1 (MQ=255)
                                                                                                ccttgtgAAAAAATGCCTT‑G‑GCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGaaa                                                        <  1:1144117/147‑1 (MQ=255)
                                                                                                     acaaaaagTGCTTT‑C‑GCGGTGTCTGCTTCCCAGTTTTTTTCCCGCCA‑AAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGaaa                                                        <  2:550665/142‑1 (MQ=255)
                                                                                                        aaaaaTGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTg                                                <  1:1215201/150‑1 (MQ=255)
                                                                                                               ccTG‑C‑GGGGCGGCTGTCTCCCAGTTTTTTCCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAAt                                        <  1:225901/151‑1 (MQ=255)
                                                                                                                  tcc‑gcgGGGCCTGCCGGCCAGCTTTTTTCCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGaa                                         <  2:176129/146‑1 (MQ=255)
                                                                                                                        ggcggtTTTCTCCCGGTTTTTTTCCCGCCA‑CAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAAt                                        <  2:156785/144‑1 (MQ=255)
                                                                                                                        ggcggtTGTCTGCCAGCTTTTTTCCCCCCC‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGaa                                         <  1:73981/143‑1 (MQ=255)
                                                                                                                                    ccactttttttcCCGCCC‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                     <  1:266082/144‑1 (MQ=255)
                                                                                                                                    ccATCTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                           <  1:1165358/145‑1 (MQ=255)
                                                                                                                                     ccgtcttttttCCCCCGAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                            <  2:532071/138‑1 (MQ=255)
                                                                                                                                       gttATTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                    <  1:579935/147‑1 (MQ=255)
                                                                                                                                       gCTTTTTCTCCCCGG‑CAACAAAAGCGGCCTTCTCTCCAGGAGTGGCGGGGCCAGA‑G‑GAATGGGTTGAAGGCCGTCATGTTGGGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                         <  2:889798/144‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                 <  1:987768/148‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCG‑AAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                 <  2:360540/148‑1 (MQ=255)
                                                                                                                                           ttttcCCGCCA‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa                <  1:1171027/149‑1 (MQ=255)
                                                                                                                                            ttttCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAACTTTTAACCCCCCCTGTTGTGTGGGGATCGCTGCTGGGGACAAATGGGGAAGGGGAGGGGGGGGGAA‑CAACTCCCGAAACCCGCCAAtt                  >  1:300546/1‑146 (MQ=255)
                                                                                                                                            tttcCCCCCG‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt              <  1:701640/150‑1 (MQ=255)
                                                                                                                                               cccgccg‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg           <  2:1310865/149‑1 (MQ=255)
                                                                                                                                                        aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa  <  2:1238899/151‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                       
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA  >  NC_000913/357740‑358038

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