Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,957,975 C→A 6.3% A445S (GCT→TCT)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,957,9750CA 6.3% 1064.3 / 28.7 282A445S (GCT→TCT) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base C (157/107);  new base A (10/8);  total (167/115)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 286 total aligned reads displayed.

TGGTTAAGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGG  >  CP009273/1957827‑1958124
                                                                                                                                                    |                                                                                                                                                     
   ttAAGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGtt                                                                                                                                                   <  1:3767121/150‑1 (MQ=255)
    tAAGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGAGAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTc                                                                                                                                                  <  2:2742120/150‑1 (MQ=255)
     aaGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTcc                                                                                                                                                 >  2:3339392/1‑150 (MQ=255)
     aaGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTcc                                                                                                                                                 >  1:4389856/1‑150 (MQ=255)
      aGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCa                                                                                                                                                <  2:3625692/150‑1 (MQ=255)
      aGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCa                                                                                                                                                <  2:3625529/150‑1 (MQ=255)
      aGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCa                                                                                                                                                <  2:1773885/150‑1 (MQ=255)
       gATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAg                                                                                                                                               <  2:2925528/150‑1 (MQ=255)
        aTGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAgg                                                                                                                                              >  2:1200068/1‑150 (MQ=255)
        aTGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAgg                                                                                                                                              >  1:2858827/1‑150 (MQ=255)
        aTGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGAGGTTCCAgg                                                                                                                                              >  2:3480825/1‑150 (MQ=255)
            gTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTa                                                                                                                                          >  2:977205/1‑150 (MQ=255)
            gTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTa                                                                                                                                          <  2:3275152/150‑1 (MQ=255)
            gTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTa                                                                                                                                          <  2:3275125/150‑1 (MQ=255)
            gTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGAGGTTCCAGGGTTa                                                                                                                                          <  1:3562531/150‑1 (MQ=255)
               gCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATc                                                                                                                                       <  2:448698/150‑1 (MQ=255)
               gCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATc                                                                                                                                       <  1:1681001/150‑1 (MQ=255)
                cTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCg                                                                                                                                      <  2:2437467/150‑1 (MQ=255)
                  accaccGTGCTGGACGCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCgcc                                                                                                                                    >  1:3685903/1‑150 (MQ=255)
                    caccGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCgccag                                                                                                                                  <  2:2518223/150‑1 (MQ=255)
                     accGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCgccagc                                                                                                                                 >  2:1335277/1‑150 (MQ=255)
                     accGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCgccagc                                                                                                                                 >  1:4248539/1‑150 (MQ=255)
                          gCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAgcg                                                                                                                            >  1:3069268/1‑150 (MQ=255)
                           cTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAgcgc                                                                                                                           >  2:3974733/1‑150 (MQ=255)
                            tGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGcc                                                                                                                          >  1:778184/1‑150 (MQ=255)
                              gCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGcccc                                                                                                                        <  1:1260132/150‑1 (MQ=255)
                              gCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGcccc                                                                                                                        <  2:4459510/150‑1 (MQ=255)
                               ccTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCg                                                                                                                       <  2:2538474/150‑1 (MQ=255)
                                 tCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCACACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGa                                                                                                                     >  2:2357480/1‑150 (MQ=255)
                                      cACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGc                                                                                                                >  1:3725511/1‑150 (MQ=255)
                                      cACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGc                                                                                                                <  2:132847/150‑1 (MQ=255)
                                       aCTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCa                                                                                                               >  1:2323632/1‑150 (MQ=255)
                                        cTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCAt                                                                                                              <  1:1471520/150‑1 (MQ=255)
                                         tgtgTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATc                                                                                                             <  2:2965664/150‑1 (MQ=255)
                                         tgtgTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATc                                                                                                             <  1:2571186/150‑1 (MQ=255)
                                             tACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATcacca                                                                                                         >  2:906583/1‑150 (MQ=255)
                                               ccccTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACt                                                                                                       <  2:3574868/150‑1 (MQ=255)
                                               ccccTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGAGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACt                                                                                                       >  1:825611/1‑150 (MQ=255)
                                               cccaTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACt                                                                                                       <  1:4001888/150‑1 (MQ=255)
                                                cccTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTg                                                                                                      >  2:1872317/1‑150 (MQ=255)
                                                 ccTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGc                                                                                                     >  1:2785303/1‑150 (MQ=255)
                                                    gAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCaa                                                                                                  >  2:3448872/1‑150 (MQ=255)
                                                       tCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAAtct                                                                                               <  2:4163952/150‑1 (MQ=255)
                                                       tCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAAtct                                                                                               >  1:1428972/1‑150 (MQ=255)
                                                        cTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAAtctc                                                                                              >  1:2203595/1‑150 (MQ=255)
                                                         tGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTcc                                                                                             >  1:3305965/1‑150 (MQ=255)
                                                              cTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGcc                                                                                        >  2:1183960/1‑150 (MQ=255)
                                                               tGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCt                                                                                       >  2:3881007/1‑150 (MQ=255)
                                                                   cTTTTAGCGCACTTTCAATCCAGATCGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtca                                                                                   <  2:4036454/150‑1 (MQ=255)
                                                                   cTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTtca                                                                                   >  2:1029491/1‑150 (MQ=255)
                                                                       tAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCaa                                                                               <  1:4043206/150‑1 (MQ=255)
                                                                          cgcACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAAt                                                                            >  1:3116735/1‑150 (MQ=255)
                                                                           gcACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATg                                                                           >  1:1046331/1‑150 (MQ=255)
                                                                            cACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGc                                                                          <  2:3193407/150‑1 (MQ=255)
                                                                                ttCAATCCAGATAGCATTCAGGCCAACTGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATa                                                                      <  2:4528141/150‑1 (MQ=255)
                                                                                 tCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAg                                                                     >  1:1761740/1‑150 (MQ=255)
                                                                                  cAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGc                                                                    <  2:2053398/150‑1 (MQ=255)
                                                                                     tCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCggg                                                                 <  1:4318395/150‑1 (MQ=255)
                                                                                     tCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCggg                                                                 <  2:3848923/150‑1 (MQ=255)
                                                                                      ccAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGc                                                                <  2:3760050/150‑1 (MQ=255)
                                                                                      ccAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGAGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGc                                                                >  2:1068829/1‑150 (MQ=255)
                                                                                        aGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCGCCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAg                                                              >  2:485535/1‑150 (MQ=255)
                                                                                           gaGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGtt                                                           <  2:1371572/149‑1 (MQ=255)
                                                                                            aGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTg                                                          >  2:1958695/1‑150 (MQ=255)
                                                                                               aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAg                                                       >  2:2763941/1‑150 (MQ=255)
                                                                                               aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAg                                                       >  2:4089622/1‑150 (MQ=255)
                                                                                               aTTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAg                                                       >  1:1106692/1‑150 (MQ=255)
                                                                                                  cAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATc                                                    >  2:1683680/1‑150 (MQ=255)
                                                                                                  cAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATc                                                    >  1:1363444/1‑150 (MQ=255)
                                                                                                  aaGGCCAAATGCCGGATCGACGGTCGCCTAACCAGGTAACGTCCCGGCAGCGGTACCAGGTTTAATCGCCAGCCAGCGCCCCGGATACGCATCACCACACCCAATCTCCACGCATATCATCAAAATCCGATAGCGGGCAGGTAGCAGaac                                                    >  2:4128920/2‑148 (MQ=255)
                                                                                                      ccAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCata                                                >  1:298150/1‑150 (MQ=255)
                                                                                                      ccAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCata                                                <  2:3644192/150‑1 (MQ=255)
                                                                                                       cAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGAGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCatat                                               <  2:825611/150‑1 (MQ=255)
                                                                                                        aaaTGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATAtt                                              >  1:678309/1‑150 (MQ=255)
                                                                                                            gCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTcg                                          >  2:620344/1‑150 (MQ=255)
                                                                                                              cGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTcgcg                                        <  2:3842247/150‑1 (MQ=255)
                                                                                                                gATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGaa                                      <  2:1980027/150‑1 (MQ=255)
                                                                                                                 aTCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAAt                                     >  1:961453/1‑150 (MQ=255)
                                                                                                                  tCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAAtg                                    <  2:680668/150‑1 (MQ=255)
                                                                                                                  tCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAAtg                                    >  2:1673250/1‑150 (MQ=255)
                                                                                                                    gACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAAtgtg                                  <  1:2248127/150‑1 (MQ=255)
                                                                                                                    gACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAAtgtg                                  <  1:1758790/150‑1 (MQ=255)
                                                                                                                     aCGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGACTGTGc                                 >  2:4314916/1‑150 (MQ=255)
                                                                                                                      cGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCACTCTCCACGCCTTTCATCAAAATGCGATAGCTGGCAGGTTGCAGATCCATATTGTCGCGAATGTGca                                >  2:1812220/1‑150 (MQ=255)
                                                                                                                      cGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGca                                <  2:4324870/150‑1 (MQ=255)
                                                                                                                       ggTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGcac                               >  1:3788778/1‑150 (MQ=255)
                                                                                                                       ggTCGCCTCACCAGGTAACGTCCCGGCAGAGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAAGGTGcac                               >  2:2702415/1‑150 (MQ=255)
                                                                                                                         tCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGcacca                             >  1:1744223/1‑150 (MQ=255)
                                                                                                                                aCCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCTAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTggcgg                      >  1:761426/1‑150 (MQ=255)
                                                                                                                                 ccAGGTAACGTCCCGGCAGAGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTggcggc                     <  2:1631227/150‑1 (MQ=255)
                                                                                                                                      tAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGaaa                >  2:3790432/1‑150 (MQ=255)
                                                                                                                                        aCGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATc              >  1:309330/1‑150 (MQ=255)
                                                                                                                                          gTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATccc            >  1:2254548/1‑150 (MQ=255)
                                                                                                                                             ccGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCAtc         <  1:2008063/150‑1 (MQ=255)
                                                                                                                                              cGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCAtct        <  1:561980/150‑1 (MQ=255)
                                                                                                                                              cGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCAtct        <  2:2682516/150‑1 (MQ=255)
                                                                                                                                                 cAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCt     <  2:4345041/150‑1 (MQ=255)
                                                                                                                                                  aGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTg    <  1:3360869/150‑1 (MQ=255)
                                                                                                                                                    cGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTggg  <  2:4286194/150‑1 (MQ=255)
                                                                                                                                                    cGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTggg  <  2:3739563/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGTTAAGATGCGTTGCTACCACCGTGCTGGCCTCAACCACTGTGTACCCCTGAATCTGCGCCTGTTCTTTTAGCGCACTTTCAATCCAGATAGCATTCAGGCCAAATGCCGGATCGACGGTCGCCTCACCAGGTAACGTCCCGGCAGCGGTTCCAGGGTTAATCGCCAGCCAGCGCCCCGGATAAGCATCACCACTGCCAATCTCCACGCCTTTCATCAAAATGCGATAGCGGGCAGGTTGCAGATCCATATTGTCGCGAATGTGCACCACTGGCGGCAGAAATCCCATCTCCTGGG  >  CP009273/1957827‑1958124

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: