Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,253,280 T→C 6.5% L216L (TTG→CTG)  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,2800TC 6.5% 623.3 / 20.8 201L216L (TTG→CTG) leuC3‑isopropylmalate dehydratase large subunit
Reads supporting (aligned to +/- strand):  ref base T (86/101);  new base C (7/6);  total (94/107)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 206 total aligned reads displayed.

GAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGAAGCAAGCGGTCGAG  >  NC_002947/2253132‑2253429
                                                                                                                                                    |                                                                                                                                                     
   gaTGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTc                                                                                                                                                   <  1:272890/150‑1 (MQ=255)
   gaTGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTc                                                                                                                                                   >  2:715292/1‑150 (MQ=255)
   gaTGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTc                                                                                                                                                   >  1:2239236/1‑150 (MQ=255)
      gaagaaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGAt                                                                                                                                                >  2:4899391/1‑150 (MQ=255)
      gaagaaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCAGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGAt                                                                                                                                                >  2:3712306/1‑150 (MQ=255)
         gaaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGa                                                                                                                                             <  1:3209109/150‑1 (MQ=255)
         gaaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGa                                                                                                                                             >  1:3629747/1‑150 (MQ=255)
         gaaCATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGa                                                                                                                                             >  2:3690964/1‑150 (MQ=255)
          accATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGaa                                                                                                                                            <  2:5017310/148‑1 (MQ=255)
             aTGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGc                                                                                                                                         >  1:4827886/1‑150 (MQ=255)
                  tgTGCGTGTTGATGGCGAACTGGCTCCAGGCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCAt                                                                                                                                    <  1:2030851/149‑1 (MQ=255)
                         gTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATc                                                                                                                             >  1:1654717/1‑150 (MQ=255)
                         gTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATc                                                                                                                             >  2:4375262/1‑150 (MQ=255)
                            gAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGc                                                                                                                          >  2:2415543/1‑150 (MQ=255)
                               ggCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAAc                                                                                                                       <  1:666870/150‑1 (MQ=255)
                                 ccAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACAt                                                                                                                     <  2:2606042/150‑1 (MQ=255)
                                  cAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATg                                                                                                                    >  2:3413422/1‑150 (MQ=255)
                                   aaTTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGt                                                                                                                   >  2:1488167/1‑150 (MQ=255)
                                     ttGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTcc                                                                                                                 <  1:2847055/150‑1 (MQ=255)
                                        gctccaggCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGAtt                                                                                                              <  2:2693456/144‑2 (MQ=255)
                                        gctccaggCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGAtt                                                                                                              >  1:2424002/7‑149 (MQ=255)
                                         ctgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCg                                                                                                             >  1:4178621/1‑150 (MQ=255)
                                         ctgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAAtt                                                                                                                                                      <  2:151245/109‑1 (MQ=255)
                                          tgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGa                                                                                                            >  1:1091266/1‑150 (MQ=255)
                                          tgccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGa                                                                                                            >  1:4912820/1‑150 (MQ=255)
                                          tccaggCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGa                                                                                                            >  1:4865797/5‑150 (MQ=255)
                                           gccGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGaa                                                                                                           <  1:3936603/150‑1 (MQ=255)
                                           ccaggCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGaa                                                                                                           <  1:3107287/147‑1 (MQ=255)
                                             cGGCGTCACCGCCAAGGGTATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGc                                                                                                         >  1:4815848/1‑150 (MQ=255)
                                              ggCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGcc                                                                                                        >  1:2710635/1‑150 (MQ=255)
                                              ggCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGcc                                                                                                        >  2:3909052/1‑150 (MQ=255)
                                               gCGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGGCg                                                                                                       <  2:2576532/150‑1 (MQ=255)
                                                cGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGGCgg                                                                                                      <  2:132916/150‑1 (MQ=255)
                                                cGTCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGGCgg                                                                                                      <  2:4260560/150‑1 (MQ=255)
                                                 gTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGgc                                                                                                     >  2:4411669/1‑150 (MQ=255)
                                                  tCACCGCCAAGGACATCACGCTTGCAATCATTGGCGAGATTGGCACGGCTGGCGGCACAGGCTATGTGATCGAATATGCTGGCAGCGCCATTCGCTCGCTGTCGATGGAAGGCCGCATGACGGTCTGCAACATGTCGATTGAAGgcggcg                                                                                                    >  1:2670278/1‑150 (MQ=255)
                                                    aCCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGcc                                                                                                  >  2:4941209/1‑150 (MQ=255)
                                                    aCCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGcc                                                                                                  >  2:775167/1‑150 (MQ=255)
                                                           aGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTgg                                                                                           <  2:2710635/150‑1 (MQ=255)
                                                             gATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGc                                                                                         >  1:3151939/1‑150 (MQ=255)
                                                              atatCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGcc                                                                                        >  1:1396095/1‑150 (MQ=255)
                                                               tatCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCt                                                                                       >  1:2217309/1‑150 (MQ=255)
                                                                 tCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtgg                                                                                     >  2:5438251/1‑150 (MQ=255)
                                                                 tCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtgg                                                                                     <  1:2506448/150‑1 (MQ=255)
                                                                   gTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtggtg                                                                                   <  1:194256/150‑1 (MQ=255)
                                                                   gTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtggtg                                                                                   <  1:4652989/150‑1 (MQ=255)
                                                                    tGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtggtgg                                                                                  <  2:5304358/150‑1 (MQ=255)
                                                                    tGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCtggtgg                                                                                  <  1:3413422/150‑1 (MQ=255)
                                                                       tGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCCTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCac                                                                               >  1:4310338/1‑149 (MQ=255)
                                                                       tGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATGCGCGAAGTGTCGCTGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCaa                                                                               >  2:2386160/1‑150 (MQ=255)
                                                                         gCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAAcc                                                                             >  1:1598716/1‑150 (MQ=255)
                                                                           cGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGa                                                                           >  2:580451/1‑150 (MQ=255)
                                                                             tGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATg                                                                         <  1:2600310/150‑1 (MQ=255)
                                                                               aTCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGcc                                                                       >  2:5205465/1‑150 (MQ=255)
                                                                                   gCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTa                                                                   >  1:622008/1‑150 (MQ=255)
                                                                                   gCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTa                                                                   >  1:4574342/1‑150 (MQ=255)
                                                                                    cAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTAc                                                                  >  2:3175045/1‑150 (MQ=255)
                                                                                    cAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTAc                                                                  >  2:2746027/1‑150 (MQ=255)
                                                                                       gATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGt                                                               >  1:4932364/1‑150 (MQ=255)
                                                                                        aTCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGtt                                                              >  1:1841286/1‑150 (MQ=255)
                                                                                         tCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTg                                                             >  2:3520724/1‑150 (MQ=255)
                                                                                         tCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTg                                                             >  1:5231197/1‑150 (MQ=255)
                                                                                         tCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTg                                                             <  2:296449/150‑1 (MQ=255)
                                                                                          cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc                                                            >  1:4347729/1‑150 (MQ=255)
                                                                                          cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc                                                            <  1:5381208/150‑1 (MQ=255)
                                                                                          cGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGc                                                            <  1:4010527/150‑1 (MQ=255)
                                                                                             cACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTa                                                         >  1:4173622/1‑150 (MQ=255)
                                                                                              aCCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTAc                                                        >  1:68924/1‑150 (MQ=255)
                                                                                                 gccgGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTc                                                     <  2:1654717/150‑1 (MQ=255)
                                                                                                        gCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGgc                                              <  1:45605/150‑1 (MQ=255)
                                                                                                            cGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGccc                                          <  2:3542634/150‑1 (MQ=255)
                                                                                                               tCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTa                                       >  2:4696964/1‑150 (MQ=255)
                                                                                                                  cGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTAcgc                                    <  1:4518048/150‑1 (MQ=255)
                                                                                                                  cGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTAcgc                                    <  2:2769555/150‑1 (MQ=255)
                                                                                                                  cGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTAcgc                                    <  2:1020721/150‑1 (MQ=255)
                                                                                                                    ccATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACgcgc                                  <  1:1695871/150‑1 (MQ=255)
                                                                                                                     cATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGcc                                 >  2:1050475/1‑150 (MQ=255)
                                                                                                                       tGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGa                               >  1:4763392/1‑150 (MQ=255)
                                                                                                                            ttCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGc                          <  2:505832/150‑1 (MQ=255)
                                                                                                                            ttCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCAACGTCGAAGGGCGCCCGTACGCGCCGAAGGGc                          <  2:220581/150‑1 (MQ=255)
                                                                                                                              cGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGa                        <  1:1502645/150‑1 (MQ=255)
                                                                                                                               gcaggcagTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAg                       <  2:2924538/150‑1 (MQ=255)
                                                                                                                                 aggcagTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCa                     <  2:5107771/150‑1 (MQ=255)
                                                                                                                                   gcagTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGt                   >  1:3776445/1‑150 (MQ=255)
                                                                                                                                   gcagTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGt                   <  1:3717107/150‑1 (MQ=255)
                                                                                                                                   gcagTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGt                   <  1:2149865/150‑1 (MQ=255)
                                                                                                                                       tGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAATGGCGAGCAGTGGaa               <  2:2456381/150‑1 (MQ=255)
                                                                                                                                       tGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaa               <  1:102707/150‑1 (MQ=255)
                                                                                                                                       tGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaa               <  1:89457/150‑1 (MQ=255)
                                                                                                                                       tGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaa               <  1:1863059/150‑1 (MQ=255)
                                                                                                                                       tGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaa               <  1:2996109/150‑1 (MQ=255)
                                                                                                                                           aTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTAc                                                                  <  2:808614/95‑1 (MQ=255)
                                                                                                                                           aTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTAc                                                                  >  1:808614/1‑95 (MQ=255)
                                                                                                                                           aTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaagcaa           >  1:1931125/1‑150 (MQ=255)
                                                                                                                                           aTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaagcaa           >  1:3188591/1‑150 (MQ=255)
                                                                                                                                           aTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaagcaa           >  1:1528725/1‑150 (MQ=255)
                                                                                                                                           aTTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGaagcaa           <  2:3992420/150‑1 (MQ=255)
                                                                                                                                               gcgAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGAAGCAAGCgg       <  2:3629747/150‑1 (MQ=255)
                                                                                                                                                cgAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGAAGCAAGCGGt      <  2:1795593/150‑1 (MQ=255)
                                                                                                                                                 gAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGAAGCAAGCGGTc     <  2:212127/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAGATGAAGAACATGCTGGTGCGCGTGGAAGGCCAATTGCCTGCCGGCGTCACCGCCAAGGATATCGTGCTGGCCGTGATCGGCAAGATCGGCACCGCCGGTGGCAACGGTCACGCCATGGAATTCGCAGGCAGTGCCATTCGCGAATTGTCGATGGAAGGCCGCATGACCATCTGCAACATGTCCATCGAAGCCGGCGCCCGCGTGGGCCTGGTGGCAACCGATGCCACTACCGTTGCCTACGTCGAAGGGCGCCCGTACGCGCCGAAGGGCGAGCAGTGGAAGCAAGCGGTCGAG  >  NC_002947/2253132‑2253429

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: