Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,971,078 C→A 5.8% A65E (GCG→GAG)  wecF → TDP‑Fuc4NAc:lipidIIFuc4NAc transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,971,0780CA 5.8% 1127.7 / 28.3 299A65E (GCG→GAG) wecFTDP‑Fuc4NAc:lipidIIFuc4NAc transferase
Reads supporting (aligned to +/- strand):  ref base C (148/133);  new base A (8/10);  total (156/143)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

TAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTGTATG  >  CP009273/3970929‑3971227
                                                                                                                                                     |                                                                                                                                                     
tAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGc                                                                                                                                                       <  1:1918054/150‑1 (MQ=255)
 aaccgaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCg                                                                                                                                                      <  2:208042/150‑1 (MQ=255)
 aaccgaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCg                                                                                                                                                      >  1:4308235/1‑150 (MQ=255)
 aaccgaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCg                                                                                                                                                      >  1:1997577/1‑150 (MQ=255)
 aaccgaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCg                                                                                                                                                      >  1:1997121/1‑150 (MQ=255)
 aaccgaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCg                                                                                                                                                      <  1:1866140/150‑1 (MQ=255)
   ccgaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGa                                                                                                                                                    >  1:1468500/1‑150 (MQ=255)
    cgacccgTTGTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGaa                                                                                                                                                   <  2:4184044/150‑1 (MQ=255)
     gaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAg                                                                                                                                                  <  1:3562676/150‑1 (MQ=255)
     gaaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAg                                                                                                                                                  >  1:3093783/1‑150 (MQ=255)
      aaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGc                                                                                                                                                 >  2:1459944/1‑150 (MQ=255)
      aaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGc                                                                                                                                                 >  2:398428/1‑150 (MQ=255)
      aaccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGc                                                                                                                                                 >  2:1456757/1‑150 (MQ=255)
      aaccgTTGTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGTCGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCATTTTTTCCCTGGGAATAAATCGCTGGCGGTAGc                                                                                                                                                 >  2:398664/1‑150 (MQ=255)
       accgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCg                                                                                                                                                >  1:1349742/1‑150 (MQ=255)
        ccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCgg                                                                                                                                               <  1:826981/150‑1 (MQ=255)
        ccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCgg                                                                                                                                               <  1:3307806/150‑1 (MQ=255)
            tttGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCAt                                                                                                                                           >  2:4137695/1‑150 (MQ=255)
            tttGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCAt                                                                                                                                           <  1:2996353/150‑1 (MQ=255)
             ttGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATc                                                                                                                                          <  1:2498650/150‑1 (MQ=255)
                  gTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaa                                                                                                                                     <  1:3054841/150‑1 (MQ=255)
                  gTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaa                                                                                                                                     <  1:1515060/150‑1 (MQ=255)
                   tttttCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGGTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaaa                                                                                                                                    >  2:2421824/1‑150 (MQ=255)
                    ttttCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCGAAAg                                                                                                                                   <  2:675312/150‑1 (MQ=255)
                      ttCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCa                                                                                                                                 <  1:1664896/150‑1 (MQ=255)
                         aaTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCaaaa                                                                                                                              <  1:1153323/150‑1 (MQ=255)
                         aaTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCaaaa                                                                                                                              >  1:3734331/1‑150 (MQ=255)
                         aaTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCGAAAGCaaaa                                                                                                                              >  2:3181227/1‑150 (MQ=255)
                            gACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCt                                                                                                                           >  1:2597283/1‑150 (MQ=255)
                                cgcTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAAcc                                                                                                                       >  2:2372619/1‑150 (MQ=255)
                                cgcTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAAcc                                                                                                                       >  2:1008006/1‑150 (MQ=255)
                                cgcTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAAcc                                                                                                                       <  2:2936243/150‑1 (MQ=255)
                                  cTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGt                                                                                                                     <  1:403374/150‑1 (MQ=255)
                                          gacgaGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCg                                                                                                             >  2:4342644/1‑150 (MQ=255)
                                            cgagcgagCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGtt                                                                                                           <  1:1784662/150‑1 (MQ=255)
                                              agcgagcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGtttt                                                                                                         <  1:55146/150‑1 (MQ=255)
                                                cgagcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGtttttt                                                                                                       <  1:1100506/150‑1 (MQ=255)
                                                 gagcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTttc                                                                                                      >  2:322040/1‑150 (MQ=255)
                                                 gagcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTttc                                                                                                      >  1:487520/1‑150 (MQ=255)
                                                   gcACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTttctt                                                                                                    <  2:3751690/150‑1 (MQ=255)
                                                     aCGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTcc                                                                                                  <  1:1451678/150‑1 (MQ=255)
                                                             gAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTGGTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTCACCGTCAGCAGCGTTTTTTCTTCCACGGTCAg                                                                                          <  2:174919/150‑1 (MQ=255)
                                                                tttATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTc                                                                                       <  2:355143/150‑1 (MQ=255)
                                                                   aTGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAAt                                                                                    <  2:2101557/150‑1 (MQ=255)
                                                                        tgttgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCcac                                                                               <  2:1837287/150‑1 (MQ=255)
                                                                           tgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACt                                                                            <  1:2859093/150‑1 (MQ=255)
                                                                                aGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggc                                                                       <  2:3458566/150‑1 (MQ=255)
                                                                                  gacgacGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggctg                                                                     >  1:42439/1‑150 (MQ=255)
                                                                                  gacgacGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggctg                                                                     <  2:2424636/150‑1 (MQ=255)
                                                                                   acgacgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGtggctgg                                                                    <  2:2833805/150‑1 (MQ=255)
                                                                                     gacgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTggctggct                                                                  <  2:2975544/150‑1 (MQ=255)
                                                                                      acgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGctc                                                                 >  1:213320/1‑150 (MQ=255)
                                                                                        ggCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctg                                                               <  2:2597283/150‑1 (MQ=255)
                                                                                          cTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctgct                                                             >  2:3384350/1‑150 (MQ=255)
                                                                                            tAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGa                                                           >  1:97037/1‑150 (MQ=255)
                                                                                            tAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGa                                                           >  1:3826820/1‑150 (MQ=255)
                                                                                              aGTGATAGCTGTCCGGCGCTTACTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCTAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAgt                                                         >  1:3687149/1‑150 (MQ=255)
                                                                                                    aGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGc                                                   <  1:1402534/150‑1 (MQ=255)
                                                                                                    aGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGc                                                   >  1:2659841/1‑150 (MQ=255)
                                                                                                     gCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCa                                                  <  2:1511313/150‑1 (MQ=255)
                                                                                                        gTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTa                                               >  2:4459160/1‑150 (MQ=255)
                                                                                                        gTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTa                                               <  2:807459/150‑1 (MQ=255)
                                                                                                        gTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTa                                               <  2:2659841/150‑1 (MQ=255)
                                                                                                         tCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTaa                                              >  1:1330524/1‑150 (MQ=255)
                                                                                                           cGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGc                                            >  1:4207907/1‑150 (MQ=255)
                                                                                                            ggCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGcc                                           >  1:460359/1‑150 (MQ=255)
                                                                                                             gcgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGccc                                          >  2:3856625/1‑150 (MQ=255)
                                                                                                              cgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCcca                                         <  1:2734114/150‑1 (MQ=255)
                                                                                                              cgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCcca                                         <  2:3124509/150‑1 (MQ=255)
                                                                                                              cgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCcca                                         <  1:2760700/150‑1 (MQ=255)
                                                                                                                cTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCccagc                                       >  2:551266/1‑150 (MQ=255)
                                                                                                                 tttCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTcagcccagcc                                      >  2:1163370/1‑150 (MQ=255)
                                                                                                                  ttCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCccagcca                                     >  2:181597/1‑150 (MQ=255)
                                                                                                                   tCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGACAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCccagccag                                    <  2:4276217/150‑1 (MQ=255)
                                                                                                                   tCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCccagccag                                    >  1:3229399/1‑150 (MQ=255)
                                                                                                                    cTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGt                                   >  1:2326326/1‑150 (MQ=255)
                                                                                                                     tgtgCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGtt                                  >  2:528963/1‑150 (MQ=255)
                                                                                                                      gtgCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGttt                                 <  2:3039999/150‑1 (MQ=255)
                                                                                                                      gtgCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGttt                                 <  2:2716895/150‑1 (MQ=255)
                                                                                                                      gtgCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGttt                                 >  2:1644025/1‑150 (MQ=255)
                                                                                                                       tgCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGtttt                                <  2:1346644/150‑1 (MQ=255)
                                                                                                                        gCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGttttt                               <  1:1008006/150‑1 (MQ=255)
                                                                                                                          aaTTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCt                             <  1:4469878/150‑1 (MQ=255)
                                                                                                                             ttttCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGc                          <  2:3735419/150‑1 (MQ=255)
                                                                                                                              tttCCCTGGGAAAAAATCGCTGGAGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCa                         >  2:1759571/1‑150 (MQ=255)
                                                                                                                                   cTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCt                    <  2:4053781/150‑1 (MQ=255)
                                                                                                                                    tGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTg                   >  1:3741378/1‑150 (MQ=255)
                                                                                                                                    tGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTg                   >  1:1476570/1‑150 (MQ=255)
                                                                                                                                    tGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTg                   >  1:4343875/1‑150 (MQ=255)
                                                                                                                                           aaaTCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCa            <  2:1981123/150‑1 (MQ=255)
                                                                                                                                            aaTCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAg           >  2:1130893/1‑150 (MQ=255)
                                                                                                                                            aaTCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAg           >  2:2684902/1‑150 (MQ=255)
                                                                                                                                             aTCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTTTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGa          >  2:1586767/1‑150 (MQ=255)
                                                                                                                                             aTCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGa          >  2:3846924/1‑150 (MQ=255)
                                                                                                                                                gCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCt       <  1:4013831/150‑1 (MQ=255)
                                                                                                                                                 cTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTg      <  2:3724188/150‑1 (MQ=255)
                                                                                                                                                  tGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTGt     >  1:2898991/1‑150 (MQ=255)
                                                                                                                                                  tGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTGt     <  2:3914694/150‑1 (MQ=255)
                                                                                                                                                   ggCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTGTa    <  2:1630252/150‑1 (MQ=255)
                                                                                                                                                     cGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTGTATg  >  2:3929829/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGCAGACCTGTATG  >  CP009273/3970929‑3971227

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: