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

AGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTT  >  NC_000913/74934‑75228
                                                                                                                                                  |                                                                                                                                                           
aGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGtt                                                                                                                                                         <  1:145010/150‑1 (MQ=255)
    ttGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTtgtttgtt                                                                                                                                                     >  2:561958/1‑150 (MQ=255)
     tgtgGGGCGTCACCGCCATAGTCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGtt                                                                                                                                                     <  2:159878/149‑1 (MQ=255)
          ggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCa                                                                                                                                              <  2:295551/151‑1 (MQ=255)
           gCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTtttttttttttat                                                                                                                                             >  2:910416/1‑146 (MQ=255)
               cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTTATaaa                                                                                                                                          >  1:401907/1‑150 (MQ=255)
                 tCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTGTTTTTTTTTTTTATaaa                                                                                                                                          >  1:964781/1‑148 (MQ=255)
                   gCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTttttttttttttttaa                                                                                                                                          >  1:317745/1‑138 (MQ=255)
                             ttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTtttttttttttat                                                                                                                                             >  1:160778/1‑128 (MQ=255)
                              tttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTGTTTTTTTTATAAAAGAAgggggaggagg                                                                                                                           >  1:1077197/1‑141 (MQ=255)
                                    tCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTATTTTGTTTTTTTTCTTAAAGAAAGCGaaagaaccatca                                                                                                                      >  1:732513/1‑139 (MQ=255)
                                     ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCa                                                                                                                 <  2:486176/151‑1 (MQ=255)
                                     ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCa                                                                                                                 <  2:756836/151‑1 (MQ=255)
                                        tAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTTTTTTTTTTTTTTATTAAAGAAGGCGCA‑GGA‑GCCCTAATTaaa                                                                                                               >  1:258866/1‑150 (MQ=255)
                                           cAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑g                                                                                                            <  2:984598/149‑1 (MQ=255)
                                            aaCAGCCCCAGCCCCGGTGTCCTGGTGCGCGGACCTTGATAACTCCCCCGCCAGTTTTGATCGCTCTCATCCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTa                                                                                                         <  1:991239/151‑1 (MQ=255)
                                             aCAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAc                                                                                                        <  1:5743/151‑1 (MQ=255)
                                             aCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTATTTTGTTTTTTTCTTTAAAGCAAGCGCA‑GTA‑GTCACAATCaaaac                                                                                                             >  1:387107/1‑145 (MQ=255)
                                                gACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑g                                                                                                      <  2:28103/149‑1 (MQ=255)
                                                 aCCCAGCCCCGGTGTACTGGTGCGCGGATCCTAATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑g                                                                                                      <  1:721638/148‑1 (MQ=255)
                                                  cccAGCCGCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑Ga                                                                                                     <  2:413440/148‑1 (MQ=255)
                                                       ccccGGTGTACTGGTGCGCGGCTCCTGATAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCa                                                                                            <  2:153843/151‑1 (MQ=255)
                                                          cGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTCTTTGTTTTTGTTATAAAACAAAGCGGA‑AGA‑GGCAAAATTAAACGGGTAC‑‑AAA‑GGGTCatt                                                                                          >  2:473885/1‑150 (MQ=255)
                                                           ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCatta                                                                                         <  2:840005/150‑1 (MQ=255)
                                                             tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCattatt                                                                                       <  2:1117658/150‑1 (MQ=255)
                                                                 cTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTc                                                                                      <  1:607903/147‑1 (MQ=255)
                                                                  tGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTACTTTGTTGTTTTCCTAAAAAAAGTGGAA‑AAA‑AGCACAAAAAAAA‑GGGAC‑GAAA‑GTGTTTTTATTcc                                                                                     >  2:347363/1‑147 (MQ=255)
                                                                    gTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTTTAAAAGAAGGGGGA‑GTA‑GCCCTTAACCAAA‑GGGAC‑GAAA‑AGGGCATTTTTTCAtccgc                                                                               >  2:600286/1‑151 (MQ=255)
                                                                       cgcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAATTTGTTTTTGTCTTAAATGAAAGAGAAGGAA‑GCCA‑AATAAAAC‑AGATC‑GAAA‑GTGGCATTATTTCAagcgcccg                                                                            >  2:103228/1‑149 (MQ=255)
                                                                        gcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGGTTTTTTTTTTTC‑TAAAAAAAAGGGAA‑GGA‑ACCCAAATTAAAA‑CGGAC‑GAAA‑AGGTCATTTTTTCAACCGCC‑cgg                                                                          >  2:1020552/1‑151 (MQ=255)
                                                                          ggtcctggATAAACCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGaa                                                                        <  1:231093/144‑1 (MQ=255)
                                                                            aTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTCATAAAACAAAGGGAA‑GGA‑GGCCTAATCAAAA‑GGGAC‑GAAA‑AGGTCATTATTTCCACCCCC‑gggg                                                                         >  1:582065/1‑147 (MQ=255)
                                                                              catGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTGTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGt                                                                     <  1:159565/148‑1 (MQ=255)
                                                                               cTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTTTTTTTTTAAAAGAAAGCGAA‑GTA‑GGCATAATCAAAA‑GGGAC‑GAAA‑GAGTCATTATTTCAGCCCCC‑GGGGACGt                                                                     >  1:951164/1‑149 (MQ=255)
                                                                                   tAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAAc                                                                 <  2:1110143/149‑1 (MQ=255)
                                                                                        cccccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATc                                                           <  1:1111594/148‑1 (MQ=255)
                                                                                           ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCt                                                       <  1:672609/151‑1 (MQ=255)
                                                                                            cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTg                                                      <  2:872419/151‑1 (MQ=255)
                                                                                            cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTGTTTTTTTGATAAAAGCAAGCGAA‑GGA‑GGCATAATTAAAC‑AGGAC‑GGAA‑GTGTCATTATTTCCGCCGCCGGGGACCGTTAAAggattcgc                                                         >  1:856764/1‑143 (MQ=255)
                                                                                            cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGCTTTTTGTTTTTTTTTTTTTCTAAACGAAGGCGAA‑GGA‑GCCATAATTAAAA‑GGGAC‑GAAA‑ATGGCCTTATTTTCACCGCC‑GGGGACGTTAAAGGCCTcccc                                                        >  1:175867/1‑147 (MQ=255)
                                                                                             gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTTTCTTTTTTCCTTTCTTTAACAAAGGCGCA‑GAA‑GACCAAAACAAAC‑AGTGC‑GAAAGGGGCCGTTTTTCCGCCCCCC‑GGGACGGTTAACggatcc                                                           >  2:134477/1‑142 (MQ=255)
                                                                                                 gTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTAGTTTTTTTTTTTCATAAAAGAAAGGGAA‑GGA‑GCCATAATCAAAA‑GGGAC‑GAAA‑GGGTCATTATTTCAGCCCCC‑GGGGACGTTAACGGCCTCCGCTGcc                                                    >  2:306704/1‑148 (MQ=255)
                                                                                                         cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCttt                                          <  2:976357/150‑1 (MQ=255)
                                                                                                             ctcAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTgg                                       <  1:910292/149‑1 (MQ=255)
                                                                                                             ctcAACCAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTTGTTTTTTCTTTTTATAAAAAAAGGAGAA‑GGAGGCAAAAATAAAAG‑GGAAG‑GAAA‑CGGT‑ATTTTTTCAGCCGCC‑GGGGAAGGTAAAGGGATCCCCCGCCACACCGCcttt                                         >  1:665580/1‑147 (MQ=255)
                                                                                                               cAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCaaa                                   <  2:530605/151‑1 (MQ=255)
                                                                                                               cAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCaa                                    <  1:1098632/150‑1 (MQ=255)
                                                                                                                  ccAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg                                <  2:1060786/151‑1 (MQ=255)
                                                                                                                  ccAGTTCTTTCAGGCTTTGTGGC‑GGGGTTTTTATTTTTTTTTTTTCTTAAAACAAGGGGAA‑GGA‑GCCCTAATCAAAA‑AGTACGGAAA‑GGGT‑ATTTTTTCAGCCCCC‑CGGGACGTTAAAGGCATCCCCTGCCCCACCCCCTTTTGCcaa                                   >  1:369531/1‑148 (MQ=255)
                                                                                                                           tCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGGCCGGttt                        <  1:906052/150‑1 (MQ=255)
                                                                                                                                  ttGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggc                 <  2:240308/150‑1 (MQ=255)
                                                                                                                                      ggC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGgcggcg             >  1:992827/1‑150 (MQ=255)
                                                                                                                                          gggTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt          >  1:158537/1‑150 (MQ=255)
                                                                                                                                            gTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc       >  2:415461/1‑151 (MQ=255)
                                                                                                                                               tttCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAgg    >  2:601867/1‑151 (MQ=255)
                                                                                                                                                  cAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGgtt  >  1:649555/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                           
AGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGC‑GGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAA‑GTA‑GCCATAATCAAAC‑GGTAC‑GAAA‑GTGTCATTATTCCAGCCGCC‑GGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTT  >  NC_000913/74934‑75228

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