Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0790CT41.1% 94.6 / 45.3 73L134L (CTG→CTA) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (4/39);  new base T (30/0);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.66e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGG  >  NC_000913/74937‑75226
                                                                                                                                               |                                                                                                                                                      
cTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGttt                                                                                                                                                   >  2:536301/1‑148 (MQ=255)
   gtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTggtttt                                                                                                                                             >  2:385505/1‑146 (MQ=255)
    tgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTtttt                                                                                                                                             >  1:317056/1‑150 (MQ=255)
     ggggTGTCATGGCCGTGGACTTTTCGCACCCATAGCAACAGACCCACCCCCGGTGTACTGGTGCGCGAATCCTGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGt                                                                                                                                           <  1:69523/151‑1 (MQ=255)
      gggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTc                                                                                                                                          <  1:838085/151‑1 (MQ=255)
          gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATaa                                                                                                                                      <  2:73309/151‑1 (MQ=255)
          gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTGtttt                                                                                                                                        >  2:139457/1‑146 (MQ=255)
            cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAAc                                                                                                                                    <  1:737450/151‑1 (MQ=255)
            cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAAc                                                                                                                                    <  1:1000676/151‑1 (MQ=255)
               cGCCATAGCCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                                 <  2:112384/151‑1 (MQ=255)
               cGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                                 <  2:275482/151‑1 (MQ=255)
               cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTCCTGGTGCGCGAATCCTGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACg                                                                                                                                   <  2:826180/149‑1 (MQ=255)
                 ccATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAgg                                                                                                                               <  1:590080/151‑1 (MQ=255)
                  cATAGACTTTTTGCATCCATAGCAACAGACACAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTCAAATTTTTCtttttgttaaaaaagg                                                                                                                               >  1:427683/1‑136 (MQ=255)
                                  ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑a                                                                                                             <  2:1038979/151‑1 (MQ=255)
                                    aTAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTATAAAAGAAGGCGAAGGAGCCATAATC‑aa                                                                                                            >  1:902584/1‑150 (MQ=255)
                                      aGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAAc                                                                                                          <  2:1023265/150‑1 (MQ=255)
                                        cAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTTCTCAACAAAAGAGAAGCAGCCATAAAC‑AAAAgg                                                                                                        >  2:125281/1‑150 (MQ=255)
                                          aCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTGTCATAAAAAAAAGCGAAGGAGCCATAATC‑aaa                                                                                                           >  1:1018061/1‑145 (MQ=255)
                                          aCAGACCCACCCCCGGTGTCCTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGt                                                                                                       <  2:730035/149‑1 (MQ=255)
                                           cAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTTTAAACAAAAGCGGAGGAGCCAAAAAC‑AAAAgg                                                                                                        >  2:571545/1‑147 (MQ=255)
                                            aGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTATTTTTTTTTTTCATAAAAAAAGGGGAAGGAACCACAATC‑aaa                                                                                                           >  2:1065760/1‑143 (MQ=255)
                                            aGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTTTTTTTGTTTTAAACAAAGGGGAAGGAGCCCTAAAC‑AAAAGGTACg                                                                                                    >  2:1124873/1‑150 (MQ=255)
                                               cccACCCCCGGTGTACTGTGGCGCGGAGCCTGATAAATTCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAg                                                                                                <  2:413275/151‑1 (MQ=255)
                                                  acccccGGTGTCCTGGTGCGCGGATCCTGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAgt                                                                                               <  1:300398/147‑1 (MQ=255)
                                                   gCCCCGGTTTTCTGGTTCCGGTATCTTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAgtgt                                                                                             <  2:856892/150‑1 (MQ=255)
                                                    cccgGGTGTAGTGGGGCGCGATTCTTCATAAAGCTCCCGCCAGTTTTAATCGCTCTCAGCCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCAGAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCa                                                                                           <  2:953720/151‑1 (MQ=255)
                                                       cGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTATTTTTTTTTTTTTTTAAAAAAAGGCGGAGGAGGCATAATT‑AAAAAGGACGaaaaagttctt                                                                                         >  2:495835/1‑142 (MQ=255)
                                                        ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTTTTTTTTTTATAAAAGAAAGCGAAGGAACCCTAATA‑AAACAGGAAGAAAGAGTGatt                                                                                         >  1:24030/1‑149 (MQ=255)
                                                        ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCatta‑t                                                                                      <  1:945630/151‑1 (MQ=255)
                                                        ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTCTTTTTTTTTTTCATAAAAAAAAGGGAAGAAGCCAAAATC‑CAAAGGTACGacaaagttatt                                                                                         >  2:1102781/1‑139 (MQ=255)
                                                          tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCatta‑tt                                                                                     <  1:625515/150‑1 (MQ=255)
                                                            tACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTTTTAAAAAAAGCGGAAGAGCCATAAAC‑AAAAAGGAAGAAAGGGTCatta‑tttc                                                                                   >  2:1127986/1‑148 (MQ=255)
                                                            tACTGGTGCGCGGACCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTcc                                                                                   <  1:46274/150‑1 (MQ=255)
                                                             aCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTTTTTTTTTCATAAACGAAGGGGGAGAAGCCCAAAAC‑AAAAGGGGCGaaaaaggtattatt                                                                                      >  1:1106227/1‑136 (MQ=255)
                                                             aCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTATTTTTTTCTTTTCTTAAAAGAAGGGGGAGTAGCCCAAATC‑AAAAGGTGAGaaaaaggtattatt                                                                                      >  1:377444/1‑136 (MQ=255)
                                                                ggggcgcGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAgc                                                                                <  2:845962/146‑1 (MQ=255)
                                                                ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAgccg                                                                              <  2:72708/151‑1 (MQ=255)
                                                                   gcgcGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAgccgcc                                                                            <  2:931874/150‑1 (MQ=255)
                                                                   gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAgccgcc                                                                            <  2:423990/150‑1 (MQ=255)
                                                                        gATCCTGAAAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAc                                                                      <  2:549933/151‑1 (MQ=255)
                                                                          tCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTTTTTTTTTTTTTCATAAAAGAAGGCGAAGGAGCCATAATT‑AAACGGGACGAAAGTGTCATTT‑TTTCAACCGCCcgg                                                                         >  1:373444/1‑146 (MQ=255)
                                                                           ccTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑Gtt                                                                  <  2:317014/151‑1 (MQ=255)
                                                                                 aaaTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTa                                                                 <  1:281729/146‑1 (MQ=255)
                                                                                 aaaTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTGTCTTAAACGAAAGCGAAAAAACCCTAAAT‑AAAAAGTAAGAAAAAGTCATTA‑TTTTAGCCCCCGGGGAC‑ggt                                                                  >  2:1046374/1‑143 (MQ=255)
                                                                                    tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTCATAAAAGAAAGGGAAGGAGCCATAATA‑AAAAGGGGCGAAAGAGTCATTA‑TTTCAGCCGCCGGGGAC‑GGTTACgg                                                             >  1:1124444/1‑147 (MQ=255)
                                                                                     cccccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATcc                                                        <  1:747635/149‑1 (MQ=255)
                                                                                     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATcc                                                        <  2:448468/151‑1 (MQ=255)
                                                                                     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATcc                                                        <  1:345280/151‑1 (MQ=255)
                                                                                      aCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTCGTTTTTTTTTTTCTTAAACGAAGGGGAAGGAACGCCAAAC‑AAACAGGAAGAAATTGTCAATATTTCTACCCCCCGGGAA‑‑GCTAACgg                                                             >  2:361013/1‑145 (MQ=255)
                                                                                       cccGCCATTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTGCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGc                                                      <  2:426413/151‑1 (MQ=255)
                                                                                        ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTCTAAAAAAAGGGGAAGGAACCATAAAC‑AAAAAGGACGAAAAAGTCATTA‑TTTCACCCGCCCGGGAA‑GCTAAAGGGATcc                                                        >  1:883951/1‑148 (MQ=255)
                                                                                          gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTTTTCTTTTTATTAAAGAAGGCGAAGGAGGCCTAATT‑AAACAGTACGAAAGAGTGATTA‑TTTCAGCCGCCCGGGAC‑GGTAACGGCAAcc                                                        >  2:547492/1‑146 (MQ=255)
                                                                                          gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTGTTTTAAACGAAAGCGAAGTAGCCATAATC‑AAACGGGACGAAAGTGTCATTA‑TTTCAGCCGCCGGGGAA‑GGTTACGGCATcccccgc                                                   >  1:62986/1‑146 (MQ=255)
                                                                                           ccAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGcc                                                  <  1:505736/151‑1 (MQ=255)
                                                                                            cAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCca                                                 <  1:926438/151‑1 (MQ=255)
                                                                                                 ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTAGTTT‑TTTTTTTTATAAAAAAACGGGGAGGAGACATAAAT‑AAAAAGGACGAAAGAGGCTTTA‑TTTCAGCCGGCGCGGAA‑GGTAAAGGCCTCCgccgcccc                                                >  1:618641/1‑145 (MQ=255)
                                                                                                  tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTCATAAAAAAAGGCGGAGTAGCCCTAATC‑AAACGGGACGAAAGAGTCATTT‑TTTCTTCCCCCCGGGAG‑ATTAAAGGCATCCGCCGCCCCAcccc                                           >  1:1121339/1‑149 (MQ=255)
                                                                                                   gATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGc                                           <  2:1132160/150‑1 (MQ=255)
                                                                                                    aTCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTTTTAAAAGAAGGCGGAGGAGCCATAAAA‑AAAAAGGGCGGAAAGGTTATTA‑TTTCCGCCGCCGGGGACGGTTAAAGG‑ATCCGCTGccccccccc                                           >  1:894073/1‑142 (MQ=255)
                                                                                                     tCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTTATAAAAAAAAGGGGAGGAGGCATAATT‑AAAACGGACGGAAGGGGCATTA‑TTTCAACCGCCCGGGAA‑GTTAAAGGCCTccccccccccccccc                                           >  2:130889/1‑135 (MQ=255)
                                                                                                       gCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTCTTTTTCTAAAAAAAGGAGAAGGAGGCATAATT‑AAAAGGTACGGAAGAGTTATTA‑TTTCTGCCGCCGGGGAA‑GGTAACCGCCTCCCCCGCCCCACCCCcttt                                       >  2:63128/1‑150 (MQ=255)
                                                                                                        ctctCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTTTTATAAAAAAAAGGGGAGAAGGCATAATT‑AAACGGGACGGAAGAGGCATTA‑TTTCAACCGCCGGGGAA‑GGTAAAGGCATCCGCTGCCACCCCCCcttt                                       >  2:818235/1‑149 (MQ=255)
                                                                                                         tctcAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTATAAACGAAAGCGGAGGAGCCCTAATCAAAACGGT‑CGGAAGGGTCCTTT‑ATTCAGCCCCCGGGGAA‑GGTAACGGCATCCGCCGCCACACCCCtttt                                       >  2:201495/1‑148 (MQ=255)
                                                                                                            cAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATGATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGc                                    <  1:537213/148‑1 (MQ=255)
                                                                                                                   ttCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc                           <  2:948300/150‑1 (MQ=255)
                                                                                                                         cAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGtttt                     <  2:993414/150‑1 (MQ=255)
                                                                                                                          aGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTAc                   <  2:502938/151‑1 (MQ=255)
                                                                                                                           ggCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACt                  <  2:645848/151‑1 (MQ=255)
                                                                                                                                       gggTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt        >  1:634949/1‑150 (MQ=255)
                                                                                                                                           ttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc     >  1:103269/1‑149 (MQ=255)
                                                                                                                                           ttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACa    <  1:489604/150‑1 (MQ=255)
                                                                                                                                            tttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCAGCGTCTAACAgg  >  1:50194/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATC‑AAACGGTACGAAAGTGTCATTA‑TTCCAGCCGCCGGGAAC‑GTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGG  >  NC_000913/74937‑75226

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