Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,448,259 T→G 14.4% L145R (CTG→CGG)  thiE → thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,448,2590TG14.4% 269.6 / 3.7 152L145R (CTG→CGG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (67/53);  new base G (20/6);  total (92/60)
Fisher's exact test for biased strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

CTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAC  >  NC_002947/5448111‑5448409
                                                                                                                                                    |                                                                                                                                                      
cTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCCGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATGCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGGCGGCCTg                                                                                                                                                       >  2:3311145/1‑150 (MQ=255)
 tCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGc                                                                                                                                                      >  1:1200447/1‑150 (MQ=255)
 tCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGc                                                                                                                                                      <  1:5293637/150‑1 (MQ=255)
   ggCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTc                                                                                                                                                    >  1:2526740/1‑150 (MQ=255)
    gCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCg                                                                                                                                                   <  1:3185199/150‑1 (MQ=255)
     ccGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGa                                                                                                                                                  <  1:4951614/150‑1 (MQ=255)
     ccGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGa                                                                                                                                                  <  1:1669000/150‑1 (MQ=255)
      cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGACGGCCTGCTCGAg                                                                                                                                                 <  2:2206737/150‑1 (MQ=255)
        atgggCAAGCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCa                                                                                                                                               >  2:2567934/4‑150 (MQ=255)
         cAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAg                                                                                                                                              <  2:2584562/150‑1 (MQ=255)
          aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGGGCCgg                                                                                                                                             >  1:1749719/1‑150 (MQ=255)
           ggCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGc                                                                                                                                            <  2:14048/150‑1 (MQ=255)
           cgCAATAATCGTCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGATTCTTCAATTCTGTCACCAACCCGTGCGCGCCCGCTGCCAAAGTCGGCCTTCTCGACCAGGc                                                                                                                                            <  2:598913/149‑1 (MQ=255)
              aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCcg                                                                                                                                         >  2:3885727/1‑150 (MQ=255)
               atcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGAGCAGGTCcgc                                                                                                                                        >  1:4244506/1‑150 (MQ=255)
                tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCTGCGCGAGCAGGCCcgcg                                                                                                                                       >  1:1383227/1‑150 (MQ=255)
                   tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGccc                                                                                                                                    >  2:873745/1‑150 (MQ=255)
                   tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGccc                                                                                                                                    <  2:4395926/150‑1 (MQ=255)
                    cGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCa                                                                                                                                   >  2:1015895/1‑150 (MQ=255)
                    cGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCa                                                                                                                                   <  2:4411676/150‑1 (MQ=255)
                     ggCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCGGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAg                                                                                                                                  <  2:5004466/150‑1 (MQ=255)
                        tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTg                                                                                                                               >  1:4907676/1‑150 (MQ=255)
                          cACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaa                                                                                                                             >  2:4859947/1‑150 (MQ=255)
                            ccTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCCAGCAGGCCCCCCCCCAGGTGAAAc                                                                                                                           >  2:5114531/1‑150 (MQ=255)
                             cggCCACGCCTGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACt                                                                                                                          <  2:4138786/148‑1 (MQ=255)
                             cTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCCGGCCCCCGCCCAGGTGGAAcg                                                                                                                          >  1:838607/1‑149 (MQ=255)
                                 cACGCCAGCCTCGAACGGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCg                                                                                                                      <  2:2155224/150‑1 (MQ=255)
                                    gCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATc                                                                                                                   <  2:5265859/150‑1 (MQ=255)
                                    gCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATc                                                                                                                   >  1:4856624/1‑150 (MQ=255)
                                     ccAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCCCCCAGGGGAAACAGCCGAGCg                                                                                                                  >  1:5234355/1‑150 (MQ=255)
                                       aGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGcc                                                                                                                >  1:4084520/1‑150 (MQ=255)
                                       aGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCCCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGcc                                                                                                                <  2:1202986/150‑1 (MQ=255)
                                       aGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGCCCACCGGAAACAGCCCAACGcc                                                                                                                >  1:1484922/1‑150 (MQ=255)
                                         ccTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCGGCGCGGGCAGGCCCGCGCCCCCGGGAAAATGCCGAACGccca                                                                                                              >  1:1148621/1‑148 (MQ=255)
                                          cTCGAACTGGCCGCCCGGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCCGGCCCGCGCCCCGGTGGAACAGCCGCTCGCCGTg                                                                                                             >  1:2145039/1‑150 (MQ=255)
                                          cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTg                                                                                                             >  2:1503799/1‑150 (MQ=255)
                                          cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGGGGAACTGCCGATCGCCggg                                                                                                             >  2:109795/1‑148 (MQ=255)
                                          cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGGGAAACTGCCGATCGCCggg                                                                                                             >  2:109768/1‑148 (MQ=255)
                                          cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGACACCAAGCCGGGCGCGCCCGCAGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCCAACGACGTg                                                                                                             >  2:2465022/1‑150 (MQ=255)
                                           tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGCGCCCGCCCGCGCCCCGGGGCAACTGCCGAGCGCCGGGa                                                                                                            >  1:390839/1‑150 (MQ=255)
                                               aCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATCGCCGTGATCgg                                                                                                        >  2:357121/1‑150 (MQ=255)
                                               aCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGCCACCAAGCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCgg                                                                                                        <  2:3427208/150‑1 (MQ=255)
                                                cTGGCCGCCCAGGCTGCCAGCGACGGCGCCAGCTACGTGGCCTTGGGCCGCTACTTCAATTCTGCCCCCACGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGCCCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                       <  1:2699488/150‑1 (MQ=255)
                                                   gccgccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGTCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCCCCCCCACGTGAAAATGCCGCTCGCCCTGGTCGGTGGc                                                                                                    >  1:1431828/1‑150 (MQ=255)
                                                       cccAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCCCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCa                                                                                                <  2:4524800/150‑1 (MQ=255)
                                                         cAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTGCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCCGGTGAACCTGCCGATCGCCGTGATCGGTGGccgcccc                                                                                              >  2:4397880/1‑144 (MQ=255)
                                                            gCTGCCCGCGATGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCATTTCTGTCACCACGCCGGGGGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTc                                                                                           <  2:390839/150‑1 (MQ=255)
                                                            gCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTc                                                                                           >  2:4513726/1‑150 (MQ=255)
                                                              tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCCCCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGa                                                                                         <  2:178492/150‑1 (MQ=255)
                                                                ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGCGAAACAGCCGATCGCCGGGATCGGTGGCATCACCCTCGACa                                                                                       >  2:70412/1‑150 (MQ=255)
                                                                   gCGACGGCGCCCGCTAGTTGGCCTTTGGCCGTTTCTTCAATTCTGCCACCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcg                                                                                    <  1:820502/150‑1 (MQ=255)
                                                                   gCGACGGCGCCCGCTAGGGGGCCTTGGGCCGCTTCTTCATTTCTGACACCCCGCGGGGCGCGCCCGCGGCCACCGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcg                                                                                    <  1:2157535/150‑1 (MQ=255)
                                                                   gCGAAGGCGCCAGCTCCGTGGCCTTTGGCCGCTTCTTCAATTCGGCCACCACGCCGGGCGCGCCCGCGGCCACCGGCGGCGTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcg                                                                                    <  2:1153200/150‑1 (MQ=255)
                                                                                    gTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCCATCGCCGTGATCGGTGGCATAACCCACGACAACGCCGCCCCACTGGTCGcc                                                                   >  2:73397/1‑150 (MQ=255)
                                                                                      ggCCTTTGGCCGCTTCTTCAATTGTGGCACCACGCGGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCa                                                                 <  1:3137073/150‑1 (MQ=255)
                                                                                               ccGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCTCCCCACAGGTCGCCCACGGCGCCGa                                                        >  2:3367066/1‑150 (MQ=255)
                                                                                                cGCTTCTTCAATACGGTCACCAAGCCGGGCGCGCCCGCTGCCAACGACGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                       <  2:3619060/150‑1 (MQ=255)
                                                                                                     cttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGGGAAACAGCCCATCGCCGTGATCCGTGGCATCACCCACGACAACGCCCCCCCACAGGTCGCCCACGGCGCCGACctgct                                                  >  2:5026977/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGTCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                 >  2:1424420/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCCCGGGAAACTGCCCATAGCCCTGATCGGTGGCCGCACCCGCGACACCCCCGCCCCCCAGGTCGCCCCCCGCGCCCCCcggctg                                                 >  1:562738/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCACGTGAAAATGCCGATCGCCGTGATCGGTGGCATCACCCACGACAAAGCCGCCCCACAGGTCGCCCACGGCGCCGCCctgctg                                                 >  1:1150732/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCACGCCCGCGCCCACGGGAAACTGCCGATCGCCCTGATCGCTGGCAACACCCTCGACAACCCCCCCCCACAGGTCGCCCACCGCGCCCAccaccag                                                 >  1:3358983/1‑145 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACCCCCCCCCACAGGTCGCCCACCGCGCCCACcggcag                                                 >  1:4413647/1‑148 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGCTCGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTCGCCCACGGCGCCGACctgcag                                                 >  2:4174197/1‑148 (MQ=255)
                                                                                                       tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTgg                                                >  1:4936542/1‑150 (MQ=255)
                                                                                                          ttccTGCCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCCCTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                             <  2:2320625/147‑1 (MQ=255)
                                                                                                                tCACCCCGCCGGGCGCGCCCGCGGCCCCCGCCGGCCGGCCCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATcc                                       <  1:3784303/150‑1 (MQ=255)
                                                                                                                  aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                     >  2:4978692/1‑150 (MQ=255)
                                                                                                                  aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATAGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                     >  2:1323140/1‑150 (MQ=255)
                                                                                                                  aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATAGGTGGCATCACCCACGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCAGCTGGCGGTGATCCAc                                     >  1:2852621/1‑150 (MQ=255)
                                                                                                                  aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGCCCACGTGAAAAAGCCGATCGCCGTGATAGGTGGCCTAAACCACGACAACCCCCCCCCACAGGTCGCCCACCGCGCCGCCCGCCTGGCGGGGTTCCAc                                     >  1:2852235/1‑150 (MQ=255)
                                                                                                                   cccccccGCGCGCCCCCGCCGCCACCGCCGCCCGGCTCGCGCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                    <  1:679128/145‑1 (MQ=255)
                                                                                                                   ccAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                    >  2:2699758/1‑150 (MQ=255)
                                                                                                                   ccAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                    >  2:291961/1‑150 (MQ=255)
                                                                                                                   ccAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                    >  2:4910624/1‑150 (MQ=255)
                                                                                                                    cacgCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                   <  2:4413647/147‑1 (MQ=255)
                                                                                                                    cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGTATCACCCGCGACAACGCCGCCCCACTGGTCGCCCACCGCGCCGACCTGCTGGCGGTGATCAACgg                                   >  1:249999/1‑150 (MQ=255)
                                                                                                                    cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGCCCCGGTGCAACTGCCGCTCGCCGGGATCGGTGGCATCACCCGCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCAGCTGGCGGTGATCCACgg                                   >  1:29858/1‑150 (MQ=255)
                                                                                                                     aaGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCGCGGc                                  >  1:378979/1‑150 (MQ=255)
                                                                                                                      aGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCACCACAACGCCGCCCCACAGGTCGCCCACCGCGCCGCCCGGCTGGCGGGGATCCACGGCt                                 >  1:3304088/1‑150 (MQ=255)
                                                                                                                       gCCGGGCGCGCCCGCGGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCtt                                <  2:1148621/150‑1 (MQ=255)
                                                                                                                       gCCGGGCGCGCCCGCCGCCTCCGTCGGCCCGCTCGAGCAGGCCCGCGCCAAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCTCCGACCTGCTGGCGGTGATCCACGGCtt                                <  1:3399273/150‑1 (MQ=255)
                                                                                                                            gcgcgcCCGCGGCCCACGTCGGCCTGCTCGCCCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCg                           <  1:2791793/150‑1 (MQ=255)
                                                                                                                                gcCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                       >  1:2914850/1‑150 (MQ=255)
                                                                                                                                 cccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGcc                      <  1:4994272/150‑1 (MQ=255)
                                                                                                                                  ccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAGCGCCGGGAGCGGTGGCATCACCCGCGACAACGCCGCCCCCCGGGTCGCCCACGGCGCCGACCGGCGGGCGGTGATCCACGGCTTGGTCGGTGCCg                     >  1:1863146/1‑150 (MQ=255)
                                                                                                                                     cgcccACCGCCGGCCTGCCCGCCCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACa                  <  1:630139/147‑1 (MQ=255)
                                                                                                                                         cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc              >  2:5395816/1‑150 (MQ=255)
                                                                                                                                         cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc              >  1:491736/1‑150 (MQ=255)
                                                                                                                                         cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc              >  2:3782488/1‑150 (MQ=255)
                                                                                                                                         cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc              >  2:3782336/1‑150 (MQ=255)
                                                                                                                                         cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCGTGTTCGGTGCCGACAgcgc              >  2:3809510/1‑150 (MQ=255)
                                                                                                                                         cAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCGCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc              >  2:2407164/1‑150 (MQ=255)
                                                                                                                                            cGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCa           >  1:1490841/1‑150 (MQ=255)
                                                                                                                                                ggCGTGCTCCCGCCGGGCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGACGACTGCGCTCAGGaa       <  2:1197675/150‑1 (MQ=255)
                                                                                                                                                  ccTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                         >  2:1077028/1‑130 (MQ=255)
                                                                                                                                                  ccTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                         <  1:1077028/130‑1 (MQ=255)
                                                                                                                                                  ccTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGt     <  2:1490841/150‑1 (MQ=255)
                                                                                                                                                    tGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa   >  2:2885265/1‑150 (MQ=255)
                                                                                                                                                    tGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACT‑GTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAc  >  1:393741/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAC  >  NC_002947/5448111‑5448409

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