Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,970,820 G→T 33.4% R193S (CGT→AGT)  motA ← proton conductor component of flagella motor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,970,8200GT33.4% 310.3 / 177.0 213R193S (CGT→AGT) motAproton conductor component of flagella motor
Reads supporting (aligned to +/- strand):  ref base G (84/58);  new base T (30/41);  total (114/99)
Fisher's exact test for biased strand distribution p-value = 2.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

Only 100 of 213 total aligned reads displayed.

CGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGC  >  CP009273/1970676‑1970964
                                                                                                                                                |                                                                                                                                                
cGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATc                                                                                                                                               >  2:2487640/1‑148 (MQ=255)
 gCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCg                                                                                                                                              <  1:83842/148‑1 (MQ=255)
    cTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCt                                                                                                                                           <  2:3358391/148‑1 (MQ=255)
    cTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCt                                                                                                                                           <  2:1596520/148‑1 (MQ=255)
    cTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCt                                                                                                                                           >  1:455812/1‑148 (MQ=255)
       catcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGaa                                                                                                                                        <  2:1821297/148‑1 (MQ=255)
       catcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGaa                                                                                                                                        <  2:17523/148‑1 (MQ=255)
        atcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAAc                                                                                                                                       >  1:1830972/1‑148 (MQ=255)
         tcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAAcc                                                                                                                                      <  1:1463422/148‑1 (MQ=255)
            ttttGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTGGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTaa                                                                                                                                   <  1:1326426/148‑1 (MQ=255)
             tttGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAAc                                                                                                                                  >  1:1167126/1‑148 (MQ=255)
                gCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAAcgcg                                                                                                                               <  1:3551416/148‑1 (MQ=255)
                gCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAAcgcg                                                                                                                               >  2:917642/1‑148 (MQ=255)
                gCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAAcgcg                                                                                                                               >  1:3144852/1‑148 (MQ=255)
                   ggTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGa                                                                                                                            <  1:134173/148‑1 (MQ=255)
                   ggTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCAGGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGa                                                                                                                            <  2:2477997/148‑1 (MQ=255)
                    gTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGaa                                                                                                                           <  1:1963679/148‑1 (MQ=255)
                    gTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGaa                                                                                                                           <  2:2558396/148‑1 (MQ=255)
                    gTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGaa                                                                                                                           >  1:3257756/1‑148 (MQ=255)
                     tttCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAAc                                                                                                                          >  2:3522838/1‑148 (MQ=255)
                        cGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGAc                                                                                                                       <  1:2844826/148‑1 (MQ=255)
                           cgcTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGAcccc                                                                                                                    >  2:768313/1‑148 (MQ=255)
                               tttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCAtt                                                                                                                >  2:398391/1‑148 (MQ=255)
                                   tGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAg                                                                                                            <  2:2184568/148‑1 (MQ=255)
                                   tGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAg                                                                                                            >  2:2270849/1‑148 (MQ=255)
                                   tGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAg                                                                                                            >  2:2270773/1‑148 (MQ=255)
                                   tGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAg                                                                                                            >  2:1063162/1‑148 (MQ=255)
                                        tAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGca                                                                                                       <  1:481347/148‑1 (MQ=255)
                                        tAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGca                                                                                                       <  1:658496/148‑1 (MQ=255)
                                        tAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGca                                                                                                       <  1:2553796/148‑1 (MQ=255)
                                        tAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTCACGCGTGAACGACCCCCATTACAGCCGca                                                                                                       >  2:2950851/1‑148 (MQ=255)
                                          aaaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaac                                                                                                     >  1:2446879/1‑148 (MQ=255)
                                           aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaaca                                                                                                    <  1:560598/148‑1 (MQ=255)
                                           aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaaca                                                                                                    >  2:775435/1‑148 (MQ=255)
                                                tCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATAcc                                                                                               >  1:1198973/1‑148 (MQ=255)
                                                  gCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaa                                                                                             <  1:2650559/148‑1 (MQ=255)
                                                    tAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAAc                                                                                           >  1:3249371/1‑148 (MQ=255)
                                                        ggggAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCg                                                                                       >  2:289451/1‑148 (MQ=255)
                                                           gAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGaa                                                                                    >  1:2164559/1‑148 (MQ=255)
                                                              ataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTg                                                                                 >  2:1687780/1‑148 (MQ=255)
                                                              ataaatCCATAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTg                                                                                 <  1:2838966/148‑1 (MQ=255)
                                                               taaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGa                                                                                >  1:2328675/1‑148 (MQ=255)
                                                                   tCCGTAAGCAAATAAAATGCCGAGTAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCCGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTcc                                                                            <  1:1624107/148‑1 (MQ=255)
                                                                    ccGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTccc                                                                           >  2:435102/1‑148 (MQ=255)
                                                                       tAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGa                                                                        >  2:2103880/1‑148 (MQ=255)
                                                                        aaGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGAc                                                                       >  2:684554/1‑148 (MQ=255)
                                                                         aGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCTTGTGCGATAAGCGCACCCAGCTCGGGGGCAGGACGATGGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGAcc                                                                      <  1:951183/148‑1 (MQ=255)
                                                                         aGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGAcc                                                                      <  2:1320377/148‑1 (MQ=255)
                                                                         aGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGAcc                                                                      >  2:331256/1‑148 (MQ=255)
                                                                            cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgc                                                                   >  2:3022681/1‑148 (MQ=255)
                                                                               tAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGcc                                                                >  1:2915205/1‑148 (MQ=255)
                                                                                 aaaTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAg                                                              >  2:3459947/1‑148 (MQ=255)
                                                                                     gCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCTGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTg                                                          >  2:3466817/1‑148 (MQ=255)
                                                                                     gCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTg                                                          >  1:900996/1‑148 (MQ=255)
                                                                                     gCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTg                                                          >  2:706916/1‑148 (MQ=255)
                                                                                      ccGTGGAAAGTCCCCACCATCGCTTGTGCTATATGCGCACCCAGCTCGGCGGCTGGACGATCGGCTGAACCTAACGCGTGAACGCCCCCCATTACAGCCGAAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGt                                                         <  1:3492628/148‑1 (MQ=255)
                                                                                           gAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGcc                                                    >  2:2675933/1‑148 (MQ=255)
                                                                                           gAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGcc                                                    >  1:226879/1‑148 (MQ=255)
                                                                                            aaaGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCg                                                   >  2:1726642/1‑148 (MQ=255)
                                                                                             aaGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCgg                                                  >  2:1838820/1‑148 (MQ=255)
                                                                                              aGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCggg                                                 <  2:1133753/148‑1 (MQ=255)
                                                                                               gTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGa                                                <  1:1749129/148‑1 (MQ=255)
                                                                                               gTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGa                                                <  1:2464364/148‑1 (MQ=255)
                                                                                               gTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGa                                                <  1:138878/148‑1 (MQ=255)
                                                                                               gTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGa                                                <  1:3239960/148‑1 (MQ=255)
                                                                                                  cccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACtt                                             <  1:3288590/148‑1 (MQ=255)
                                                                                                  cccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACtt                                             >  1:2274475/1‑148 (MQ=255)
                                                                                                  cccACCATCGCAGGGGCGATACGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACtt                                             <  1:2487747/148‑1 (MQ=255)
                                                                                                    caccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCt                                           >  1:1261490/1‑148 (MQ=255)
                                                                                                    caccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCt                                           >  2:2574329/1‑148 (MQ=255)
                                                                                                     accaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTg                                          >  2:2370563/1‑148 (MQ=255)
                                                                                                       caTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGcc                                        >  2:1492249/1‑148 (MQ=255)
                                                                                                         tCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTc                                      >  1:1906076/1‑148 (MQ=255)
                                                                                                           gCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGc                                    <  1:1687780/148‑1 (MQ=255)
                                                                                                             aTGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCtt                                  <  1:2036525/148‑1 (MQ=255)
                                                                                                              tgtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCttt                                 >  1:1316693/1‑148 (MQ=255)
                                                                                                              tgtgCGATAAGCGCACCCAGATCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCttg                                 >  2:1319699/1‑147 (MQ=255)
                                                                                                               gtgCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                >  2:3152337/1‑148 (MQ=255)
                                                                                                               gtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                <  2:232058/148‑1 (MQ=255)
                                                                                                               gtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTc                                >  1:1597517/1‑148 (MQ=255)
                                                                                                                     tAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGt                          >  1:746546/1‑148 (MQ=255)
                                                                                                                      aaGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtc                         >  2:1196201/1‑148 (MQ=255)
                                                                                                                       aGCGCACCCAGCTCGGCGGCATGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtct                        <  1:1039327/148‑1 (MQ=255)
                                                                                                                        gcgcACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCACACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtctc                       <  1:2397516/148‑1 (MQ=255)
                                                                                                                        gcgcACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtctc                       <  1:1554576/148‑1 (MQ=255)
                                                                                                                        gcgcACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtctc                       >  1:2119499/1‑148 (MQ=255)
                                                                                                                            aCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCACACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtc                   <  1:2980666/148‑1 (MQ=255)
                                                                                                                            aCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtc                   <  2:2810916/148‑1 (MQ=255)
                                                                                                                            aCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtc                   <  1:2270790/148‑1 (MQ=255)
                                                                                                                            aCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtc                   <  1:2980056/148‑1 (MQ=255)
                                                                                                                            aCCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAAtc                   >  1:299392/1‑148 (MQ=255)
                                                                                                                               cAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATctct                >  2:625182/1‑148 (MQ=255)
                                                                                                                                aGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCtt               >  2:44191/1‑148 (MQ=255)
                                                                                                                                 gCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTc              >  2:604203/1‑148 (MQ=255)
                                                                                                                                 gCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTc              <  1:3314541/148‑1 (MQ=255)
                                                                                                                                  cTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTca             >  2:1040499/1‑148 (MQ=255)
                                                                                                                                       cggcAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTcatccat        <  2:1333845/148‑1 (MQ=255)
                                                                                                                                          cAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCag     >  1:2208621/1‑148 (MQ=255)
                                                                                                                                            ggACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCagag   <  1:3370177/148‑1 (MQ=255)
                                                                                                                                             gACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGc  >  1:24216/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGC  >  CP009273/1970676‑1970964

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: