Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,253,231 C→T 6.4% A199A (GCC→GCT leuC → 3‑isopropylmalate dehydratase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,253,2310CT 6.4% 631.1 / 21.3 204A199A (GCC→GCTleuC3‑isopropylmalate dehydratase large subunit
Reads supporting (aligned to +/- strand):  ref base C (106/85);  new base T (7/6);  total (113/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 206 total aligned reads displayed.

GCACTTCCGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTC  >  NC_002947/2253083‑2253378
                                                                                                                                                    |                                                                                                                                                   
   cTTCCGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTg                                                                                                                                                 <  2:5398601/150‑1 (MQ=255)
   cTTCCGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTg                                                                                                                                                 <  1:4536552/150‑1 (MQ=255)
    ttCCGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTgg                                                                                                                                                <  2:2878731/150‑1 (MQ=255)
    ttCCGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTgg                                                                                                                                                >  1:505832/1‑150 (MQ=255)
       cGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCaa                                                                                                                                             >  1:5164723/1‑150 (MQ=255)
       cGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCaa                                                                                                                                             >  1:904150/1‑150 (MQ=255)
        gAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAAc                                                                                                                                            >  2:4518048/1‑150 (MQ=255)
              gAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCAc                                                                                                                                      >  2:3325037/1‑150 (MQ=255)
              gAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCAc                                                                                                                                      >  2:3717107/1‑150 (MQ=255)
                    gTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATg                                                                                                                                <  1:3806149/150‑1 (MQ=255)
                    gTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATg                                                                                                                                >  2:1863059/1‑150 (MQ=255)
                     tGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATgg                                                                                                                               >  1:220581/1‑150 (MQ=255)
                      gCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGACGAACATGCTGGTGCGCGTGGAAGGCCAAGTGCCTGCCGGCGTCACCGCCAAGGATGTCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGa                                                                                                                              <  2:3142196/150‑1 (MQ=255)
                      gCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGa                                                                                                                              >  2:1608651/1‑150 (MQ=255)
                         cGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAAtt                                                                                                                           >  2:3936603/1‑150 (MQ=255)
                         cGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAAtt                                                                                                                           <  1:3253540/150‑1 (MQ=255)
                         cGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAAtt                                                                                                                           >  1:2924613/1‑150 (MQ=255)
                          gCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTc                                                                                                                          >  2:102707/1‑150 (MQ=255)
                          gCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTc                                                                                                                          >  1:1020721/1‑150 (MQ=255)
                          gCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTc                                                                                                                          >  2:4915589/1‑150 (MQ=255)
                          gCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTc                                                                                                                          >  2:245820/1‑150 (MQ=255)
                          gCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGCCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCACCGGTCACGCCATGGAATTc                                                                                                                          >  2:1671370/1‑150 (MQ=255)
                            cACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCgc                                                                                                                        >  1:1307697/1‑150 (MQ=255)
                             aCCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCgca                                                                                                                       >  2:3458641/1‑150 (MQ=255)
                             aCCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCgca                                                                                                                       <  2:172980/150‑1 (MQ=255)
                              cccAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCgcag                                                                                                                      <  2:4976091/150‑1 (MQ=255)
                                cAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCgcaggc                                                                                                                    >  1:212288/1‑150 (MQ=255)
                                 cgTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGTTATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCgcaggca                                                                                                                   <  1:1660646/149‑1 (MQ=255)
                                            gCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTcgc                                                                                                        <  1:2791138/150‑1 (MQ=255)
                                            gCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTcgc                                                                                                        >  2:4670855/1‑150 (MQ=255)
                                            gCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTcgc                                                                                                        <  1:1341370/150‑1 (MQ=255)
                                               aagaagATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGaa                                                                                                     <  1:3483581/150‑1 (MQ=255)
                                               aagaagATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTAACGCCATGGAATTCGCAGGCAGTGCCATTCGCGaa                                                                                                     <  2:1173646/150‑1 (MQ=255)
                                                    gaTGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTc                                                                                                <  1:272890/150‑1 (MQ=255)
                                                        aagaaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATg                                                                                            <  2:2704779/150‑1 (MQ=255)
                                                           aaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGaa                                                                                         <  2:739235/150‑1 (MQ=255)
                                                              aTGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGc                                                                                      <  1:3200932/150‑1 (MQ=255)
                                                                  tGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCa                                                                                  >  1:4794700/1‑150 (MQ=255)
                                                                   tgTGCGTGTTGATGGCGAACTGGCTCCAGGCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCAt                                                                                 <  1:2030851/149‑1 (MQ=255)
                                                                   ggTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCAt                                                                                 <  2:367000/150‑1 (MQ=255)
                                                                   ggTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCAt                                                                                 <  2:5451051/150‑1 (MQ=255)
                                                                       cgcgTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGAcc                                                                             >  2:2823369/1‑150 (MQ=255)
                                                                        gcgTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGTACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCa                                                                            <  1:2398399/150‑1 (MQ=255)
                                                                         cgTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCAt                                                                           <  2:5164723/150‑1 (MQ=255)
                                                                          gTGGAAGGCCAAGTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATc                                                                          <  1:3410174/150‑1 (MQ=255)
                                                                           tGGAAGGCCACTTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCt                                                                         <  2:789652/150‑1 (MQ=255)
                                                                           tGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCt                                                                         <  2:521733/150‑1 (MQ=255)
                                                                            ggAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTg                                                                        >  2:4220517/1‑150 (MQ=255)
                                                                             gAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGc                                                                       >  2:2415543/1‑150 (MQ=255)
                                                                               aGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCaa                                                                     >  1:533071/1‑150 (MQ=255)
                                                                                  ccAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACAt                                                                  <  2:2606042/150‑1 (MQ=255)
                                                                                    aaTTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGt                                                                <  1:4586151/150‑1 (MQ=255)
                                                                                     aTTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTc                                                               <  1:3690964/150‑1 (MQ=255)
                                                                                      ttGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTcc                                                              <  1:2847055/150‑1 (MQ=255)
                                                                                        gcctgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCAt                                                            <  1:4859381/150‑1 (MQ=255)
                                                                                         gctccaggCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGAtt                                                           >  1:2424002/7‑149 (MQ=255)
                                                                                         cctgccGGCGTCACCGCCAAGGATATCGTGCTGACCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATc                                                           <  2:4890160/150‑1 (MQ=255)
                                                                                          ctgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCg                                                          >  1:4178621/1‑150 (MQ=255)
                                                                                          ctgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAAtt                                                                                                   >  1:151245/1‑109 (MQ=255)
                                                                                          ctgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCg                                                          >  1:5426902/1‑150 (MQ=255)
                                                                                           tgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGa                                                         >  1:1091266/1‑150 (MQ=255)
                                                                                            gccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGaa                                                        <  1:3936603/150‑1 (MQ=255)
                                                                                            ccaggCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGaa                                                        <  1:3107287/147‑1 (MQ=255)
                                                                                              cGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGc                                                      >  2:2720067/1‑150 (MQ=255)
                                                                                              cGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGc                                                      >  1:4816019/1‑150 (MQ=255)
                                                                                               ggCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGcc                                                     >  1:2710635/1‑150 (MQ=255)
                                                                                                gCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGGCg                                                    <  2:2576532/150‑1 (MQ=255)
                                                                                                 cGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGGCgg                                                   <  2:132916/150‑1 (MQ=255)
                                                                                                  gTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGgc                                                  >  2:4411669/1‑150 (MQ=255)
                                                                                                     aCCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGcc                                               >  2:2532170/1‑150 (MQ=255)
                                                                                                          cAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGGCGGCGCAcgcgc                                          >  2:5039710/1‑149 (MQ=255)
                                                                                                              gATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGc                                      >  1:3151939/1‑150 (MQ=255)
                                                                                                               atatCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGcc                                     <  2:5006160/150‑1 (MQ=255)
                                                                                                               atatCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGcc                                     >  1:1396095/1‑150 (MQ=255)
                                                                                                                tatCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCt                                    >  1:2217676/1‑150 (MQ=255)
                                                                                                                  tCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtgg                                  <  1:2506448/150‑1 (MQ=255)
                                                                                                                    gTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtggtg                                <  1:194256/150‑1 (MQ=255)
                                                                                                                        tGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCaa                            <  2:1492310/150‑1 (MQ=255)
                                                                                                                        tGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATGCGCGAAGTGTCGCTGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCaa                            >  2:2386160/1‑150 (MQ=255)
                                                                                                                          gCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAAcc                          >  1:1598716/1‑150 (MQ=255)
                                                                                                                            cGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGa                        <  2:3692108/150‑1 (MQ=255)
                                                                                                                            cGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGa                        >  2:580451/1‑150 (MQ=255)
                                                                                                                            cGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGa                        <  1:4893069/150‑1 (MQ=255)
                                                                                                                            cGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGACCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGa                        >  2:4298612/1‑150 (MQ=255)
                                                                                                                              tGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATg                      <  1:2600310/150‑1 (MQ=255)
                                                                                                                               gATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGc                     >  1:1926251/1‑150 (MQ=255)
                                                                                                                                    gCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTa                >  1:622027/1‑150 (MQ=255)
                                                                                                                                     cAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTAc               >  2:3175045/1‑150 (MQ=255)
                                                                                                                                      aaGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTAcc              >  2:689031/1‑150 (MQ=255)
                                                                                                                                        gATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGt            >  1:4932364/1‑150 (MQ=255)
                                                                                                                                         aTCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGtt           <  2:2003621/150‑1 (MQ=255)
                                                                                                                                         aTCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGtt           >  1:1841286/1‑150 (MQ=255)
                                                                                                                                          tCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTg          >  2:3520724/1‑150 (MQ=255)
                                                                                                                                          tCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTg          <  2:296449/150‑1 (MQ=255)
                                                                                                                                           cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc         <  1:4010527/150‑1 (MQ=255)
                                                                                                                                           cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc         <  2:82142/150‑1 (MQ=255)
                                                                                                                                           cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc         >  1:4347729/1‑150 (MQ=255)
                                                                                                                                           cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc         <  2:2239236/150‑1 (MQ=255)
                                                                                                                                              cACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTa      >  1:4173622/1‑150 (MQ=255)
                                                                                                                                               aCCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTAc     >  1:68924/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCACTTCCGAGGTCGAGCACGTGCTCGCCACCCAGTGCCTGGTCGCCAAGAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTC  >  NC_002947/2253083‑2253378

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