Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 425,143 G→T 57.9% T61T (ACG→ACT secF → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273425,1430GT57.9% 28.4 / 311.2 200T61T (ACG→ACTsecFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (39/45);  new base T (56/59);  total (96/104)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 210 total aligned reads displayed.

CACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGG  >  CP009273/424994‑425292
                                                                                                                                                     |                                                                                                                                                     
cACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACg                                                                                                                                                       <  1:51920/150‑1 (MQ=255)
cACGGCCGGAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACg                                                                                                                                                       <  2:3669680/150‑1 (MQ=255)
   ggCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACtctc                                                                                                                                                    <  2:3620545/150‑1 (MQ=255)
   ggCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCCGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACtctc                                                                                                                                                    >  1:2456247/1‑150 (MQ=255)
     ccGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGa                                                                                                                                                  >  1:2717979/1‑150 (MQ=255)
     ccGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGa                                                                                                                                                  <  2:3569010/150‑1 (MQ=255)
       gTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaa                                                                                                                                                >  1:3699514/1‑150 (MQ=255)
       gTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaa                                                                                                                                                >  1:1640889/1‑150 (MQ=255)
            gTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCg                                                                                                                                           <  1:1312530/150‑1 (MQ=255)
             tCTATGACTTGATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCgg                                                                                                                                          <  2:4166870/150‑1 (MQ=255)
              cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGc                                                                                                                                         >  1:2322458/1‑150 (MQ=255)
               tATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCt                                                                                                                                        <  2:2704397/150‑1 (MQ=255)
                aTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTg                                                                                                                                       >  1:3931103/1‑150 (MQ=255)
                aTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTg                                                                                                                                       <  2:1878461/150‑1 (MQ=255)
                aTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTg                                                                                                                                       >  1:1635327/1‑150 (MQ=255)
                      ttATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTg                                                                                                                                 >  2:1083912/1‑150 (MQ=255)
                      ttATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTg                                                                                                                                 <  2:3010722/150‑1 (MQ=255)
                        aTGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGAc                                                                                                                               >  2:2005744/1‑150 (MQ=255)
                        aTGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGAc                                                                                                                               <  1:346441/150‑1 (MQ=255)
                        aTGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGAc                                                                                                                               <  1:2427735/150‑1 (MQ=255)
                         tGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACg                                                                                                                              >  2:883641/1‑150 (MQ=255)
                           cgcTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTCACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTa                                                                                                                            >  2:2620462/1‑150 (MQ=255)
                            gcTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTaa                                                                                                                           >  1:3302960/1‑150 (MQ=255)
                             cTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAAt                                                                                                                          >  1:4160215/1‑150 (MQ=255)
                               gggACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGc                                                                                                                        >  2:244619/1‑150 (MQ=255)
                                ggACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCg                                                                                                                       <  1:1181522/150‑1 (MQ=255)
                                  actactGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTg                                                                                                                     <  2:3516872/150‑1 (MQ=255)
                                   ctactGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGa                                                                                                                    <  2:3328848/150‑1 (MQ=255)
                                    tactGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGAt                                                                                                                   >  1:667580/1‑150 (MQ=255)
                                       tGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCa                                                                                                                >  1:2246941/1‑150 (MQ=255)
                                            tttCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCa                                                                                                           <  2:1067599/150‑1 (MQ=255)
                                            tttCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCa                                                                                                           <  1:3419293/150‑1 (MQ=255)
                                               cGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTCATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCaaaa                                                                                                        >  2:1835162/1‑150 (MQ=255)
                                                ggCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCaaaaa                                                                                                       >  2:3083456/1‑150 (MQ=255)
                                                   aTCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGcc                                                                                                    >  1:2914106/1‑150 (MQ=255)
                                                   aTCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGcc                                                                                                    >  2:2055482/1‑150 (MQ=255)
                                                        tGGTCTGCTGTTAATCGCTGATATCGTAATTATGGGCGTGCGCGGCATTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTTCAAAAAGAAGGttt                                                                                               >  2:1974500/1‑150 (MQ=255)
                                                         ggTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGtttt                                                                                              <  1:2490876/150‑1 (MQ=255)
                                                         ggTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGtttt                                                                                              >  1:3562244/1‑150 (MQ=255)
                                                          gTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTg                                                                                             <  2:3911153/150‑1 (MQ=255)
                                                          gTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGGACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTg                                                                                             <  2:1184182/150‑1 (MQ=255)
                                                           tCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGa                                                                                            <  1:309075/150‑1 (MQ=255)
                                                           tCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGa                                                                                            <  1:2391175/150‑1 (MQ=255)
                                                             tgctgTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCTCCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAag                                                                                          <  1:4127595/150‑1 (MQ=255)
                                                              gctgTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAaga                                                                                         <  1:1505233/150‑1 (MQ=255)
                                                               ctgTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAagag                                                                                        >  1:2862332/1‑150 (MQ=255)
                                                               ctgTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAagag                                                                                        >  1:2862316/1‑150 (MQ=255)
                                                                 gTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGcc                                                                                      <  1:2184183/150‑1 (MQ=255)
                                                                  ttAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCg                                                                                     <  2:2720677/150‑1 (MQ=255)
                                                                   tAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGa                                                                                    >  1:2828556/1‑150 (MQ=255)
                                                                        gctgctATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGAtgctg                                                                               >  1:3099227/1‑150 (MQ=255)
                                                                        gctgctATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGAtgctg                                                                               >  2:1572993/1‑150 (MQ=255)
                                                                         ctgctATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGAtgctgc                                                                              >  1:2649714/1‑150 (MQ=255)
                                                                          tgctATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCa                                                                             >  2:156789/1‑150 (MQ=255)
                                                                                cGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACtt                                                                       >  1:856631/1‑150 (MQ=255)
                                                                                 gTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACttt                                                                      >  1:116165/1‑150 (MQ=255)
                                                                                 gTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACttt                                                                      >  1:642864/1‑150 (MQ=255)
                                                                                  ttattaTGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTg                                                                     <  2:682315/150‑1 (MQ=255)
                                                                                  ttattaTGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTg                                                                     <  2:2779817/150‑1 (MQ=255)
                                                                                  ttattaTGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTg                                                                     <  2:378837/150‑1 (MQ=255)
                                                                                  ttattaTGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTg                                                                     <  1:2165213/150‑1 (MQ=255)
                                                                                    attatGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGt                                                                   >  2:126011/1‑150 (MQ=255)
                                                                                         gggCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTagcag                                                              <  1:3341194/150‑1 (MQ=255)
                                                                                          ggCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTagcagc                                                             <  1:2833849/150‑1 (MQ=255)
                                                                                           gCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGcc                                                            <  2:2392152/150‑1 (MQ=255)
                                                                                            cGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCa                                                           <  2:3726209/150‑1 (MQ=255)
                                                                                             gTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCAt                                                          >  1:404226/1‑150 (MQ=255)
                                                                                                cgcgGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAc                                                       >  1:3513645/1‑150 (MQ=255)
                                                                                                cgcgGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAc                                                       <  2:3645078/150‑1 (MQ=255)
                                                                                                 gcgGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGAca                                                      >  2:3122782/1‑150 (MQ=255)
                                                                                                      tttAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATg                                                 <  1:1083848/150‑1 (MQ=255)
                                                                                                         aaCTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTc                                              >  2:519610/1‑150 (MQ=255)
                                                                                                          aCTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTcc                                             >  2:4095699/1‑150 (MQ=255)
                                                                                                            tGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGt                                           <  1:3749052/150‑1 (MQ=255)
                                                                                                              ggggCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTAt                                         <  1:1924992/150‑1 (MQ=255)
                                                                                                               gggCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATg                                        <  1:3304609/150‑1 (MQ=255)
                                                                                                               gggCTGGATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATg                                        <  1:1212516/150‑1 (MQ=255)
                                                                                                                 gCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgcc                                      <  2:4083489/150‑1 (MQ=255)
                                                                                                                   tGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATgccgc                                    >  1:3593439/1‑150 (MQ=255)
                                                                                                                     gATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCct                                  >  2:484752/1‑150 (MQ=255)
                                                                                                                     gATTTCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCct                                  <  2:1427415/150‑1 (MQ=255)
                                                                                                                        ttCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCctgct                               <  2:146009/150‑1 (MQ=255)
                                                                                                                        ttCACCGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCctgct                               <  2:146124/150‑1 (MQ=255)
                                                                                                                         tCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCctgctg                              >  2:2350914/1‑150 (MQ=255)
                                                                                                                            ccGGTGGTACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAg                           <  2:3891595/150‑1 (MQ=255)
                                                                                                                                tggtACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGa                       <  2:771952/150‑1 (MQ=255)
                                                                                                                                tggtACGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGa                       <  2:4015198/150‑1 (MQ=255)
                                                                                                                                    aCGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAAcc                   <  2:2687837/150‑1 (MQ=255)
                                                                                                                                    aCGGTTATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAAcc                   <  2:2328740/150‑1 (MQ=255)
                                                                                                                                    aCGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAAcc                   <  1:1344235/150‑1 (MQ=255)
                                                                                                                                      ggTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACcgg                 <  1:1330317/150‑1 (MQ=255)
                                                                                                                                        ttATTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACcggcg               >  1:1948862/1‑150 (MQ=255)
                                                                                                                                        ttATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACcggcg               >  1:923339/1‑150 (MQ=255)
                                                                                                                                          aTTGAAATTACTCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGt             >  1:3259069/1‑150 (MQ=255)
                                                                                                                                          aTTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGt             <  1:690560/150‑1 (MQ=255)
                                                                                                                                           ttGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTc            >  2:3875427/1‑150 (MQ=255)
                                                                                                                                               aaTTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGgt        >  1:2841707/1‑150 (MQ=255)
                                                                                                                                               aaTTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGgt        >  2:434449/1‑150 (MQ=255)
                                                                                                                                                aTTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGgtg       >  2:2100363/1‑150 (MQ=255)
                                                                                                                                                     gCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTggg  <  2:2932851/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCTGAAATTGACGTAATGCGTGATGCATTGCAAAAAGCCGGTTTTGAAGAGCCGATGCTGCAAAACTTTGGTAGCAGCCATGACATCATGGTCCGTATGCCGCCTGCTGAAGGCGAAACCGGCGGTCAGGTGTTGGG  >  CP009273/424994‑425292

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: