Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,681,855 C→A 30.4% D227Y (GAC→TAC)  bcsC ← cellulose synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,681,8550CA30.4% 656.3 / 251.6 315D227Y (GAC→TAC) bcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base C (119/100);  new base A (40/56);  total (159/156)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

ACCTTGCGCACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTG  >  CP009273/3681706‑3682004
                                                                                                                                                     |                                                                                                                                                     
  cTTGCGCACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCtt                                                                                                                                                     <  1:1943652/150‑1 (MQ=255)
    tGCGCACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTa                                                                                                                                                   <  1:3518679/150‑1 (MQ=255)
    tGCGCACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTa                                                                                                                                                   <  2:4080499/150‑1 (MQ=255)
      cgcACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAAt                                                                                                                                                 <  1:3666054/150‑1 (MQ=255)
      cgcACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACGAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAAt                                                                                                                                                 <  1:4114059/150‑1 (MQ=255)
       gcACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATc                                                                                                                                                <  1:1473025/150‑1 (MQ=255)
           gAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGcc                                                                                                                                            >  2:1661982/1‑150 (MQ=255)
                 ggAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTAcc                                                                                                                                      >  2:4010104/1‑150 (MQ=255)
                 ggAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTAcc                                                                                                                                      >  2:1524251/1‑150 (MQ=255)
                  gAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCa                                                                                                                                     >  1:4433544/1‑150 (MQ=255)
                      gCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGAtt                                                                                                                                 <  2:4231899/150‑1 (MQ=255)
                      gCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGAtt                                                                                                                                 <  2:3406704/150‑1 (MQ=255)
                        aGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGAtttt                                                                                                                               >  2:1954541/1‑150 (MQ=255)
                        aGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGAtttt                                                                                                                               <  1:810798/150‑1 (MQ=255)
                         ggATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTa                                                                                                                              >  1:4414521/1‑150 (MQ=255)
                          gATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTag                                                                                                                             >  2:4377484/1‑150 (MQ=255)
                            tCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTagag                                                                                                                           <  2:219435/150‑1 (MQ=255)
                                  aGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCTAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCt                                                                                                                     <  1:3841804/150‑1 (MQ=255)
                                    cTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTc                                                                                                                   <  2:1300456/150‑1 (MQ=255)
                                    cTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTc                                                                                                                   <  1:901695/150‑1 (MQ=255)
                                     tGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCAATAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGTCGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTcg                                                                                                                  <  1:212577/150‑1 (MQ=255)
                                       tttttGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTcgcg                                                                                                                <  2:2669913/150‑1 (MQ=255)
                                         tttGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGcc                                                                                                              >  1:4588843/1‑150 (MQ=255)
                                           tGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCg                                                                                                            >  1:1541353/1‑150 (MQ=255)
                                             cTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGc                                                                                                          >  1:2271808/1‑150 (MQ=255)
                                             cTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGc                                                                                                          <  1:637888/150‑1 (MQ=255)
                                                  cTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCg                                                                                                     <  1:878227/150‑1 (MQ=255)
                                                     ccAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAAACTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGGTCGAtt                                                                                                  >  2:2305712/1‑150 (MQ=255)
                                                       aGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGAtttt                                                                                                >  1:2937771/1‑150 (MQ=255)
                                                         ttGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGc                                                                                              >  2:88434/1‑150 (MQ=255)
                                                           gCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCa                                                                                            >  1:3613710/1‑150 (MQ=255)
                                                           gCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCa                                                                                            >  1:1289570/1‑150 (MQ=255)
                                                              aTTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCt                                                                                         <  2:2760430/150‑1 (MQ=255)
                                                              aTTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCt                                                                                         <  1:1869917/150‑1 (MQ=255)
                                                              aTTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCt                                                                                         >  2:1383673/1‑150 (MQ=255)
                                                              aTTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCt                                                                                         <  1:2383105/150‑1 (MQ=255)
                                                                tGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCAGGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGt                                                                                       >  2:1125791/1‑150 (MQ=255)
                                                                    cAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCa                                                                                   <  1:1751163/150‑1 (MQ=255)
                                                                        cgccACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGAc                                                                               <  2:1568675/150‑1 (MQ=255)
                                                                           cACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGc                                                                            >  2:1967300/1‑150 (MQ=255)
                                                                             cGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCaa                                                                          <  1:3898730/150‑1 (MQ=255)
                                                                              gCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCaaa                                                                         <  2:1827718/150‑1 (MQ=255)
                                                                                  tCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCt                                                                     <  2:4596078/150‑1 (MQ=255)
                                                                                   cGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCtt                                                                    <  1:1346058/150‑1 (MQ=255)
                                                                                    gCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTc                                                                   <  1:2227758/150‑1 (MQ=255)
                                                                                     ccATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCg                                                                  >  1:733505/1‑150 (MQ=255)
                                                                                      cATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGt                                                                 <  2:4143967/150‑1 (MQ=255)
                                                                                        tCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCa                                                               >  2:3691829/1‑150 (MQ=255)
                                                                                        tCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCa                                                               <  1:2053838/150‑1 (MQ=255)
                                                                                        tCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCa                                                               <  1:3729153/150‑1 (MQ=255)
                                                                                         cACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCAc                                                              >  2:4024681/1‑150 (MQ=255)
                                                                                          aCTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACg                                                             <  2:3539045/150‑1 (MQ=255)
                                                                                              aaGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGa                                                         >  2:560222/1‑150 (MQ=255)
                                                                                              aaGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGa                                                         >  2:3315878/1‑150 (MQ=255)
                                                                                               aGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGAt                                                        >  2:4555091/1‑150 (MQ=255)
                                                                                                gATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATc                                                       <  1:98671/150‑1 (MQ=255)
                                                                                                   cGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTATAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGct                                                    >  1:254218/1‑150 (MQ=255)
                                                                                                   cGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGct                                                    <  2:305715/150‑1 (MQ=255)
                                                                                                    gagagaTATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGcta                                                   >  1:77970/1‑150 (MQ=255)
                                                                                                      gagaTATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGctact                                                 <  2:2411497/150‑1 (MQ=255)
                                                                                                       agaTATTTTTTCAGCGCCGACACACGGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGctacta                                                <  1:3254859/150‑1 (MQ=255)
                                                                                                         atatTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTaaa                                              <  1:2366323/150‑1 (MQ=255)
                                                                                                         atatTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTaaa                                              <  1:984366/150‑1 (MQ=255)
                                                                                                          tatTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAAc                                             >  2:1456822/1‑150 (MQ=255)
                                                                                                           atTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACa                                            >  1:1096390/1‑150 (MQ=255)
                                                                                                           atTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACa                                            >  1:1709561/1‑150 (MQ=255)
                                                                                                                 tCAGCGCCGACACACTGGCATCACTGTCGGGCATGTCTTTAATCTGCCCGTACCTGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAATACCGTCGTCTCGGCGATCGGTCCTAAACAGCAATa                                      >  1:316450/1‑150 (MQ=255)
                                                                                                                  cAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAg                                     >  1:2692026/1‑150 (MQ=255)
                                                                                                                  cAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAg                                     >  1:612533/1‑150 (MQ=255)
                                                                                                                     cgcCGACACACTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAgcgc                                  >  2:2954876/1‑150 (MQ=255)
                                                                                                                      gcCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGcc                                 <  2:294885/150‑1 (MQ=255)
                                                                                                                        cGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAg                               >  1:3443264/1‑150 (MQ=255)
                                                                                                                         gACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGa                              >  1:1425264/1‑150 (MQ=255)
                                                                                                                         gACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGa                              >  2:34760/1‑150 (MQ=255)
                                                                                                                          acacacTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGAt                             <  2:2062231/150‑1 (MQ=255)
                                                                                                                          acacacTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGAt                             >  2:2621100/1‑150 (MQ=255)
                                                                                                                           cacacTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGAtt                            >  2:62514/1‑150 (MQ=255)
                                                                                                                            acacTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTg                           >  2:649648/1‑150 (MQ=255)
                                                                                                                              acTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGtt                         <  2:3972791/150‑1 (MQ=255)
                                                                                                                              acTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGtt                         <  1:4329390/150‑1 (MQ=255)
                                                                                                                              acTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGtt                         >  1:3397059/1‑150 (MQ=255)
                                                                                                                               cTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGttt                        <  1:448313/150‑1 (MQ=255)
                                                                                                                               cTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGttt                        <  1:615982/150‑1 (MQ=255)
                                                                                                                               cTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGttt                        >  2:3778415/1‑150 (MQ=255)
                                                                                                                               cTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGttt                        <  1:2934475/150‑1 (MQ=255)
                                                                                                                               cTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGttt                        <  1:2311515/150‑1 (MQ=255)
                                                                                                                               cTGGCATCACTGACGGGCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGttt                        <  2:1289570/150‑1 (MQ=255)
                                                                                                                                tGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGtttt                       >  2:3162772/1‑150 (MQ=255)
                                                                                                                                 ggCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTg                      >  2:1017019/1‑150 (MQ=255)
                                                                                                                                     tCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAgg                  <  1:1019486/150‑1 (MQ=255)
                                                                                                                                         tGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCg              <  2:4454526/150‑1 (MQ=255)
                                                                                                                                             gggCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTAtt          <  1:4377484/150‑1 (MQ=255)
                                                                                                                                             gggCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTAtt          <  2:1800270/150‑1 (MQ=255)
                                                                                                                                             gggCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTAtt          <  2:3201637/150‑1 (MQ=255)
                                                                                                                                              ggCATGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTg         <  2:4278548/150‑1 (MQ=255)
                                                                                                                                                 aTGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGccc      >  1:825461/1‑150 (MQ=255)
                                                                                                                                                 aTGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGccc      <  1:4179380/150‑1 (MQ=255)
                                                                                                                                                 aTGTATTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGccc      <  2:4409602/150‑1 (MQ=255)
                                                                                                                                                  tGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCg     <  1:2618381/150‑1 (MQ=255)
                                                                                                                                                  tGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCg     <  1:400893/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCTTGCGCACGAGCGCGGAAAGCAGGATCGGCCAGCTGTTTTTGCTGTTCTGCCAGTTGCGATTGCGCAGCCGCCACGCTATCGCCATCACTAAAGATCGAGAGATATTTTTTCAGCGCCGACACACTGGCATCACTGACGGGCATGTCTTTAATCTGCCCGTACCAGATTTTAGAGGCCCCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTG  >  CP009273/3681706‑3682004

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: ‑