Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0790CT43.8% 74.0 / 33.9 64L134L (CTG→CTA) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (4/32);  new base T (28/0);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCT  >  NC_000913/74931‑75220
                                                                                                                                                      |                                                                                                                                               
gCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTc                                                                                                                                                 <  2:782232/149‑1 (MQ=255)
    ggCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTg                                                                                                                                           <  1:75060/151‑1 (MQ=255)
       ttGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTAGttgt                                                                                                                                           >  1:533654/1‑147 (MQ=255)
        tgtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTc                                                                                                                                        <  2:1024232/150‑1 (MQ=255)
           ggggCGTCATCGCCATAGGCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTg                                                                                                                                     <  1:763248/150‑1 (MQ=255)
           ggggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCTCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGt                                                                                                                                    <  1:846682/151‑1 (MQ=255)
             ggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCa                                                                                                                                  <  2:121183/151‑1 (MQ=255)
                gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTGt                                                                                                                                          >  1:515736/1‑140 (MQ=255)
                  cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTGTtttttttct                                                                                                                                 >  2:421117/1‑142 (MQ=255)
                   aTCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAtttttttttttttctaaa                                                                                                                              >  1:954496/1‑141 (MQ=255)
                   aTCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTtttttttttttctaaa                                                                                                                              >  2:606469/1‑141 (MQ=255)
                     cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                          <  2:612331/151‑1 (MQ=255)
                     cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTTTTTGTTTTTTTTATaaa                                                                                                                              >  2:799302/1‑147 (MQ=255)
                      gCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGGTTTCTTTTTTCTTTAACGaa                                                                                                                          >  1:822137/1‑150 (MQ=255)
                       ccATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                          <  2:31790/149‑1 (MQ=255)
                       ccATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTCATAAAAGaa                                                                                                                          >  2:3783/1‑149 (MQ=255)
                        cATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTTTTTTTTCTTTTCATTAAAGaa                                                                                                                          >  1:686883/1‑148 (MQ=255)
                        cATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTtttttttttttcttaaagaaag                                                                                                                        >  2:580973/1‑136 (MQ=255)
                           aGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGCTTTTTTATTTTTTCTTTTCTTACAAGaaaggta                                                                                                                     >  2:552826/1‑145 (MQ=255)
                            gatattgtGCATCCAAAGCAACGGACCCAACCCAGGTGACGT‑GGTGTGCAGACCCTGAAAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCGTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGaa                                                                                                                    <  1:945735/143‑1 (MQ=255)
                              cTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTTTAAACCAAggggga                                                                                                                    >  2:971700/1‑144 (MQ=255)
                               tttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCATAAAAAAAGGGGAAg                                                                                                                   >  2:137389/1‑148 (MQ=255)
                                     cATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GATGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGCTTTTTCCCTTGTTCTTGTCATAAACCCAGGcccagccggc                                                                                                              >  1:108938/1‑138 (MQ=255)
                                       tCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTATAAAAAAAAGCGAAGGAGCCCTa                                                                                                           >  2:638834/1‑148 (MQ=255)
                                             gCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTATAAACAAAAGCGCAGTAGGAATAATaaaa                                                                                                     >  1:294433/1‑148 (MQ=255)
                                               aaCAGACCCAGCCCCGGTGTCCT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACg                                                                                                   <  1:804705/148‑1 (MQ=255)
                                                  aGACCCAGCCCCGGTG‑ACTGGGTGCGCGGATCCTGAGAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGTGAAGTAGCCATAATCAAACGGTa                                                                                                <  1:821031/148‑1 (MQ=255)
                                                      ccAGCCCCGGTTTACT‑GGTGCGCGGATCCTTATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTTTTGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAAc                                                                                                    <  2:563686/140‑1 (MQ=255)
                                                      ccAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgt                                                                                         <  2:827209/151‑1 (MQ=255)
                                                         gCCCGGGTGTACT‑GGTGCGCGGTTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAATAGCCATAATCAAACGGTACGAAAgtgt                                                                                       <  1:474603/150‑1 (MQ=255)
                                                         gCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTGTCTTAAACAAAGGGGGAGTAGGCCCAATCAAAAGGGACggcagag                                                                                        >  1:515420/1‑143 (MQ=255)
                                                         gCCCCGGTGTACG‑GGTGCGCGGATCCGGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgtgt                                                                                       <  1:226925/150‑1 (MQ=255)
                                                          ccccGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCa                                                                                     <  2:130300/151‑1 (MQ=255)
                                                           cccGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTCATAAACGAAGGCGAAGTAGCCAAAATTAAACAGGATGAAAGGGTCat                                                                                    >  1:691142/1‑151 (MQ=255)
                                                                tgtACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTTTTTTTGTTATAAAAGAAAGGGAAGGAGCCATAATCAAACAGGACGaaaagggcattttt                                                                                >  2:873586/1‑139 (MQ=255)
                                                                  tACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTcc                                                                             <  2:722221/150‑1 (MQ=255)
                                                                   tccgggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTc                                                                              <  1:605522/145‑1 (MQ=255)
                                                                       ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAgccg                                                                        <  2:496465/151‑1 (MQ=255)
                                                                         tgcggggccctggATAGATCGCCCGCCAGTTTTGATCGCTGTCAACGATGTCTTTCAGGCTTTGTGGC‑GGGTTTTGCAGTTTGAGCGGGTCCTAAGGGAAGACGAAGTAGCCATATTCAAGCGGGACGAAAGTGTGATCA‑TTCCAtcc                                                                         <  1:562951/136‑1 (MQ=255)
                                                                          gcgcGGACCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAgccg‑c                                                                      <  2:161164/149‑1 (MQ=255)
                                                                            gcgGATCTTGATAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCggg                                                                  <  1:391096/151‑1 (MQ=255)
                                                                              ggAGCCTAATAAATCTCCCCCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGa                                                                 <  1:964197/150‑1 (MQ=255)
                                                                                    tGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCCTAAACCAAGGGGCAGGAGGCCTAATTAAACAGGACGGAAGTGTCCTTA‑TTTTAGCCC‑CCggg                                                                  >  1:780059/1‑143 (MQ=255)
                                                                                        aaaTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTGTTCTTTTCATAAAAGAAGGCGAAGGAGCCATAATTAAACGGGACGAAAGTGTCATTATTTCCACCCC‑CCGGGAA‑GTTaa                                                          >  2:494560/1‑147 (MQ=255)
                                                                                           tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTTTAAACCAAGGCGAAGGGAGCATAATTCAACGAGGTGAAAAGGTCATTA‑ATCCAGCCG‑CCGGGCACGTTAACGGCaac                                                   >  1:212704/1‑149 (MQ=255)
                                                                                           tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTCAGTTTTTTTTTTTCATAAAACAAAGCGAAGAAGCCCAAATCAAACGGGACGAAAGGGTCATTA‑TTCCAGCCG‑CCGGGGACCTTAAAGGCATc                                                   >  1:115408/1‑151 (MQ=255)
                                                                                                cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTCGGTTTTTTTTTTCATAAAAGCAGGGGGAGGGGGCATTATCAAAAGGGACGGAAATGTCATTA‑TTTCAGCCG‑CCGGGGAAATTTACGGCATcccc                                                >  1:758205/1‑147 (MQ=255)
                                                                                                 gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTCGTTTGTTTTTTTTATAAAAGAAAGCGGAGGAGCCATAATCAAACGGGACGAAAAAGTGAATA‑TTCCCACCG‑CCGGGAACGGTAACGGCAAcc                                                  >  2:482359/1‑146 (MQ=255)
                                                                                                   cAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCca                                           <  1:918560/151‑1 (MQ=255)
                                                                                                    aGTTTGGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTGCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCca                                           <  2:873655/150‑1 (MQ=255)
                                                                                                        ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACAcc                                       <  2:581329/150‑1 (MQ=255)
                                                                                                        ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGCCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACAcc                                       <  2:629440/150‑1 (MQ=255)
                                                                                                        ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCTTAAACGAAGGCGAAGGAGCCATAATCAAAAGGGACGAAAGAGGCATTA‑TTTCAACCG‑CCGGGGACGGTAAAGGGAACCGCTGCCACAcc                                       >  2:700021/1‑150 (MQ=255)
                                                                                                             cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGGTTTTTCTTTTTATTAACTAAGGCGAAGGAGTCATAATCAAACAGTAGGAAAGAGTGATTT‑TTTCAGCCC‑CCGGGGACAGTAACAGGATCCCCAGcccaaccaccttt                                 >  1:850034/1‑140 (MQ=255)
                                                                                                               ctctCAACCATTTCTTTCAGGCTTTGTCGC‑TGTTTTTTCAGTTTATTCTTGTCATAAACAAAGGCGAAGTAGCCATAATCAAACCCTACGAAACTCTCATTA‑TTCCAGCCC‑CCGAAAACGTTAACCTCATCCCCTCCCACACCCCtat                                  <  1:86097/148‑3 (MQ=255)
                                                                                                               ctctCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTCAGTTTTTTTTTTTCATAAAAGAAGGCGAAGTAGCCATAATTAAAAGGTACGAAAGGGTCATTA‑TTTCAGCCG‑CCCGGGACGGTTACGGCATCCCCTGCCCCACCGCctttt                                >  2:656544/1‑149 (MQ=255)
                                                                                                                     aCCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACa                         <  2:996369/151‑1 (MQ=255)
                                                                                                                        aGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTGTCATAAACCAAGGCGGAAGAGGCACAATCAAAAAGGACGAAAGGGTTATTA‑TTTCAGCCGCCCGGGAAGGTTAA‑GGGATCCGCAGCCCCACCCCCTTTTGCaaaaaat                       >  2:701223/1‑146 (MQ=255)
                                                                                                                         gTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTAGTTTGTTTTTTTCATAAACGAAGGGGGAGGAGGCCTAATCAAAAAGGACGAAAAAGTCATTA‑TTTCAGCCG‑CCGGGAAAGTTAACGGCATCCGCCGCCACACCCCCTTTTGCAAACAGAcc                     >  2:1397/1‑151 (MQ=255)
                                                                                                                              ttCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGt                  <  2:60978/149‑1 (MQ=255)
                                                                                                                                      ttGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggc         <  2:790869/150‑1 (MQ=255)
                                                                                                                                        gtgGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcggc      <  2:324247/151‑1 (MQ=255)
                                                                                                                                            c‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt  <  2:1004805/151‑1 (MQ=255)
                                                                                                                                            c‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt  >  1:116146/1‑151 (MQ=255)
                                                                                                                                            c‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGCCAAACAGTCCGGTTTTACTGgcggcg     >  1:121574/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACT‑GGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTA‑TTCCAGCCG‑CCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCT  >  NC_000913/74931‑75220

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