Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,970,820 G→T 41.1% 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,8200GT41.1% 205.7 / 206.5 199R193S (CGT→AGT) motAproton conductor component of flagella motor
Reads supporting (aligned to +/- strand):  ref base G (66/51);  new base T (52/29);  total (119/80)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 200 total aligned reads displayed.

TGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTT  >  CP009273/1970673‑1970966
                                                                                                                                                   |                                                                                                                                                  
 gACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGa                                                                                                                                                   <  2:1864814/148‑1 (MQ=255)
  aCGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTAt                                                                                                                                                  <  2:1408870/148‑1 (MQ=255)
     cACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCgg                                                                                                                                               <  2:2433817/148‑1 (MQ=255)
      aCTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGc                                                                                                                                              <  1:3359730/148‑1 (MQ=255)
          catcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGaa                                                                                                                                          <  1:3558981/148‑1 (MQ=255)
            tcatTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCAGCGCATGTGCGATAAGCGCACCCCCCTCGGCGGCAGGACTATCCGCTGAAcc                                                                                                                                        >  2:3050702/1‑148 (MQ=255)
               ttttGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTaa                                                                                                                                     <  2:1480273/148‑1 (MQ=255)
                    cTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAATGCGt                                                                                                                                <  2:2498328/148‑1 (MQ=255)
                          tCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGa                                                                                                                          <  1:2759078/148‑1 (MQ=255)
                              cgcTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGAcccc                                                                                                                      >  1:3265242/1‑148 (MQ=255)
                              cgcTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGAcccc                                                                                                                      >  1:1466829/1‑148 (MQ=255)
                                 ttttCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCAt                                                                                                                   <  2:1762261/148‑1 (MQ=255)
                                       gACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGACCGACCCCCATTACAGc                                                                                                             >  2:878575/1‑148 (MQ=255)
                                       gACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGc                                                                                                             >  2:1371424/1‑148 (MQ=255)
                                       gACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGc                                                                                                             <  1:421046/148‑1 (MQ=255)
                                        aCGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGAGAACGACCCCCATTACAGcc                                                                                                            >  2:3253829/1‑148 (MQ=255)
                                           tAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGca                                                                                                         >  1:1548001/1‑148 (MQ=255)
                                           tAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGca                                                                                                         <  1:179569/148‑1 (MQ=255)
                                              aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcacca                                                                                                      >  2:2346784/1‑148 (MQ=255)
                                              aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaaca                                                                                                      >  2:1360353/1‑148 (MQ=255)
                                              aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaaca                                                                                                      <  2:1080215/148‑1 (MQ=255)
                                              aaCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaaca                                                                                                      >  2:855510/1‑148 (MQ=255)
                                               aCAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGcaacaa                                                                                                     <  1:5937/148‑1 (MQ=255)
                                                cAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAAt                                                                                                    >  2:3756303/1‑148 (MQ=255)
                                                   tCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATAcc                                                                                                 >  1:3796968/1‑148 (MQ=255)
                                                    cGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCa                                                                                                >  2:3214920/1‑148 (MQ=255)
                                                    cGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCa                                                                                                >  2:3213110/1‑148 (MQ=255)
                                                     gCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACGaa                                                                                               <  2:1128836/148‑1 (MQ=255)
                                                     gCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACaaa                                                                                               >  1:1034331/1‑148 (MQ=255)
                                                      cTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCaaa                                                                                              <  2:3706281/148‑1 (MQ=255)
                                                       tAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAAc                                                                                             >  1:1949835/1‑148 (MQ=255)
                                                        aaTGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACg                                                                                            >  2:1677822/1‑148 (MQ=255)
                                                        aaTGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACg                                                                                            >  2:2178558/1‑148 (MQ=255)
                                                          tGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGcc                                                                                          <  1:3114212/148‑1 (MQ=255)
                                                          tGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGcc                                                                                          >  1:1265136/1‑148 (MQ=255)
                                                             ggAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGa                                                                                       >  1:786019/1‑148 (MQ=255)
                                                                aataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGt                                                                                    <  2:856315/148‑1 (MQ=255)
                                                                aataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGt                                                                                    <  2:2584698/148‑1 (MQ=255)
                                                                 ataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTg                                                                                   <  2:602350/148‑1 (MQ=255)
                                                                 ataaatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTg                                                                                   >  2:541703/1‑148 (MQ=255)
                                                                  taaatCCGTAAGCCAATAAAATGCCGAGTAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGa                                                                                  <  1:1712556/148‑1 (MQ=255)
                                                                  taaatCCGTAAGCCAATAAAATGCCGAGGAAAGACCCCACCTTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGa                                                                                  <  1:1712603/148‑1 (MQ=255)
                                                                    aatCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGt                                                                                >  1:833054/1‑148 (MQ=255)
                                                                     atCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAACGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTc                                                                               <  2:3076294/148‑1 (MQ=255)
                                                                       ccGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTACCGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTccc                                                                             <  2:62206/148‑1 (MQ=255)
                                                                       ccGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTccc                                                                             <  2:1993225/148‑1 (MQ=255)
                                                                          tAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGACGTGAGTCCCCGa                                                                          >  1:3065024/1‑148 (MQ=255)
                                                                               cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgc                                                                     >  1:2301677/1‑148 (MQ=255)
                                                                               cAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgc                                                                     >  2:155838/1‑148 (MQ=255)
                                                                                aaTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcg                                                                    >  2:3197800/1‑148 (MQ=255)
                                                                                 aTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcgc                                                                   >  2:2063758/1‑148 (MQ=255)
                                                                                 aTAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAgcgc                                                                   >  2:2063647/1‑148 (MQ=255)
                                                                                    aaaTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAg                                                                >  2:2848572/1‑148 (MQ=255)
                                                                                     aaTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGa                                                               >  1:2114715/1‑148 (MQ=255)
                                                                                     aaTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATAGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGa                                                               >  2:2105496/1‑148 (MQ=255)
                                                                                      aTGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGAc                                                              <  1:3676433/148‑1 (MQ=255)
                                                                                       tGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                             >  2:1896736/1‑148 (MQ=255)
                                                                                       tGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                             >  2:1469411/1‑148 (MQ=255)
                                                                                       tGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                             >  1:2921505/1‑148 (MQ=255)
                                                                                       tGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                             >  2:2713947/1‑148 (MQ=255)
                                                                                       tGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                             >  1:2018906/1‑148 (MQ=255)
                                                                                       tGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACt                                                             <  1:2427359/148‑1 (MQ=255)
                                                                                         ccGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTTAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGt                                                           <  1:1683985/148‑1 (MQ=255)
                                                                                         ccGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGt                                                           >  1:1574595/1‑148 (MQ=255)
                                                                                            aGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCg                                                        >  2:3322822/1‑148 (MQ=255)
                                                                                             ggAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGc                                                       >  1:2401642/1‑148 (MQ=255)
                                                                                               aaaGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCg                                                     >  1:156408/1‑148 (MQ=255)
                                                                                               aaaGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCg                                                     >  2:163949/1‑148 (MQ=255)
                                                                                                aaGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCgg                                                    >  2:2230707/1‑148 (MQ=255)
                                                                                                   tCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGAc                                                 >  2:1712436/1‑148 (MQ=255)
                                                                                                   tCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGAc                                                 >  2:2827702/1‑148 (MQ=255)
                                                                                                    ccccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACt                                                >  1:1849718/1‑148 (MQ=255)
                                                                                                    ccccACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACt                                                >  1:2027738/1‑148 (MQ=255)
                                                                                                      ccaccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTc                                              >  2:396959/1‑148 (MQ=255)
                                                                                                      ccaccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTc                                              >  2:2652886/1‑148 (MQ=255)
                                                                                                         ccaTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGc                                           >  1:2236326/1‑148 (MQ=255)
                                                                                                          caTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGcc                                          <  1:1333636/148‑1 (MQ=255)
                                                                                                           aTCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCt                                         >  2:731001/1‑148 (MQ=255)
                                                                                                            tCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTc                                        >  1:3573449/1‑148 (MQ=255)
                                                                                                             cGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCg                                       >  1:2600682/1‑148 (MQ=255)
                                                                                                               cATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCt                                     <  1:127422/148‑1 (MQ=255)
                                                                                                                 tgtgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCttt                                   <  2:133775/148‑1 (MQ=255)
                                                                                                                   tgCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCg                                 <  2:625554/148‑1 (MQ=255)
                                                                                                                    gCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGt                                >  1:1631846/1‑148 (MQ=255)
                                                                                                                    gCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGt                                <  1:1873581/148‑1 (MQ=255)
                                                                                                                    gCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGt                                <  1:1873683/148‑1 (MQ=255)
                                                                                                                       aTAAGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCg                             >  2:1225126/1‑148 (MQ=255)
                                                                                                                       aTAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCg                             >  1:913346/1‑148 (MQ=255)
                                                                                                                         aaGCGCACCCAGCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGtc                           >  2:3590953/1‑148 (MQ=255)
                                                                                                                            cgcACCCCGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCa                        <  1:2008029/148‑1 (MQ=255)
                                                                                                                                  cAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATctct                  >  1:3545161/1‑148 (MQ=255)
                                                                                                                                  cAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATctct                  <  1:2723993/148‑1 (MQ=255)
                                                                                                                                   aGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCtt                 >  1:51580/1‑148 (MQ=255)
                                                                                                                                    gCTCGGCGGCAGGACTATCGGCTTAACATAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCACACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTc                <  1:2254098/148‑1 (MQ=255)
                                                                                                                                    gCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTc                >  2:3561717/1‑148 (MQ=255)
                                                                                                                                    gCTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTc                <  1:323467/148‑1 (MQ=255)
                                                                                                                                     cTCGGCGGCAGGACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTca               >  1:542285/1‑148 (MQ=255)
                                                                                                                                               ggACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCagag     >  1:2733982/1‑148 (MQ=255)
                                                                                                                                                gACTATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGc    >  1:3813212/1‑148 (MQ=255)
                                                                                                                                                gACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGc    <  2:2512378/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGACGCACTGCATCATTTTGCTGGTTTCGGCGCTTTTCTGACGTAAAACAGTCGCTAATGGGGAAATAAATCCGTAAGCCAATAAAATGCCGAGGAAAGTCCCCACCATCGCATGTGCGATAAGCGCACCCAGCTCGGCGGCAGGACGATCGGCTGAACCTAACGCGTGAACGACCCCCATTACAGCCGCAACAATACCAAACGCCGGAAGTGAGTCCCCGACCAGCGCCAGACTGTTCGCCGGGACTTCTGCCTCGCTTTCGTGCGTCTCAATCTCTTCATCCATCAGAGCTT  >  CP009273/1970673‑1970966

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: