Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,510,311 C→A 10.0% G119C (GGT→TGT)  damX ← cell division protein that binds to the septal ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,510,3110CA10.0% 1005.0 / 64.3 285G119C (GGT→TGT) damXcell division protein that binds to the septal ring
Reads supporting (aligned to +/- strand):  ref base C (118/138);  new base A (9/20);  total (127/158)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 315 total aligned reads displayed.

CTGAGTATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTT  >  CP009273/3510162‑3510460
                                                                                                                                                     |                                                                                                                                                     
ctgAGTATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAaccgatacc                                                                                                                                                       >  1:3902574/1‑150 (MQ=255)
   aGTATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAaccgataccgat                                                                                                                                                    <  2:3643730/150‑1 (MQ=255)
    gTATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatg                                                                                                                                                   <  2:2950643/150‑1 (MQ=255)
    gTATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatg                                                                                                                                                   <  1:3062462/150‑1 (MQ=255)
     tATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatga                                                                                                                                                  >  1:3876209/1‑150 (MQ=255)
     tATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTATCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatga                                                                                                                                                  >  1:1926188/1‑150 (MQ=255)
      aTTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatgat                                                                                                                                                 >  2:2472508/1‑150 (MQ=255)
      aTTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatgat                                                                                                                                                 <  2:2973562/150‑1 (MQ=255)
      aTTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCgatgat                                                                                                                                                 >  2:3587624/1‑150 (MQ=255)
       ttttCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATc                                                                                                                                                <  2:3674165/150‑1 (MQ=255)
        tttCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCa                                                                                                                                               <  1:2697536/150‑1 (MQ=255)
        tttCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCa                                                                                                                                               <  1:3341070/150‑1 (MQ=255)
           ccGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAAca                                                                                                                                            >  2:3179417/1‑150 (MQ=255)
           ccGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAAca                                                                                                                                            <  2:1254388/150‑1 (MQ=255)
           ccGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAAca                                                                                                                                            >  2:2263435/1‑150 (MQ=255)
             gCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAAcagc                                                                                                                                          >  1:169441/1‑150 (MQ=255)
                 aCGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTa                                                                                                                                      >  2:2271429/1‑150 (MQ=255)
                   gTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGa                                                                                                                                    >  1:3633238/1‑150 (MQ=255)
                      gTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAAcc                                                                                                                                 <  2:829159/150‑1 (MQ=255)
                       ttCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCa                                                                                                                                <  1:585031/150‑1 (MQ=255)
                       ttCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCa                                                                                                                                >  2:3554931/1‑150 (MQ=255)
                       ttCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCa                                                                                                                                <  1:4494814/150‑1 (MQ=255)
                       ttCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCa                                                                                                                                <  1:2045009/150‑1 (MQ=255)
                           cGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAAt                                                                                                                            <  2:1474173/150‑1 (MQ=255)
                           cGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAAt                                                                                                                            <  2:1938626/150‑1 (MQ=255)
                            ggCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATg                                                                                                                           <  1:1804019/150‑1 (MQ=255)
                             gcgcTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGc                                                                                                                          <  2:1152568/150‑1 (MQ=255)
                             gcgcTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGc                                                                                                                          >  2:539663/1‑150 (MQ=255)
                              cgcTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGcc                                                                                                                         >  1:1261494/1‑150 (MQ=255)
                              cgcTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGcc                                                                                                                         >  1:2312291/1‑150 (MQ=255)
                               gctggctgCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCg                                                                                                                        <  2:846642/150‑1 (MQ=255)
                                     tgCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGccc                                                                                                                  <  1:308668/150‑1 (MQ=255)
                                       cACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCcat                                                                                                                >  2:2605575/1‑150 (MQ=255)
                                        aCGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCcatc                                                                                                               <  1:3760755/150‑1 (MQ=255)
                                        aCGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCcatc                                                                                                               <  1:99/150‑1 (MQ=255)
                                           ccATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCcatcatc                                                                                                            <  2:3428926/150‑1 (MQ=255)
                                           ccATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCcatcatc                                                                                                            <  2:3940841/150‑1 (MQ=255)
                                             aTTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCcatcatcat                                                                                                          <  1:3465074/150‑1 (MQ=255)
                                                cGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCAtata                                                                                                       <  1:4189815/150‑1 (MQ=255)
                                                cGCCTGATCGGTCGCAGTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCAtata                                                                                                       <  2:3480995/150‑1 (MQ=255)
                                                  ccTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACt                                                                                                     >  1:2550825/1‑150 (MQ=255)
                                                  ccTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACt                                                                                                     <  1:398158/150‑1 (MQ=255)
                                                  ccTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACt                                                                                                     <  1:398696/150‑1 (MQ=255)
                                                   cTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTg                                                                                                    >  2:2516663/1‑150 (MQ=255)
                                                          gTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGaa                                                                                             <  1:677594/150‑1 (MQ=255)
                                                          gTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGaa                                                                                             >  2:4514763/1‑150 (MQ=255)
                                                           tCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAg                                                                                            <  2:1079011/150‑1 (MQ=255)
                                                                  gCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGttt                                                                                     >  2:2578687/1‑150 (MQ=255)
                                                                    cAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTAc                                                                                   >  2:3877364/1‑150 (MQ=255)
                                                                    cAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTAc                                                                                   <  2:2479826/150‑1 (MQ=255)
                                                                    cAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTGTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTAc                                                                                   >  1:2941464/1‑150 (MQ=255)
                                                                       cAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTgg                                                                                >  1:1675658/1‑150 (MQ=255)
                                                                        aaGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGgc                                                                               <  1:2575612/150‑1 (MQ=255)
                                                                           aTCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCCTCATATACTGACGAGAAGCGGGGTTACTGgctgc                                                                            >  2:3026438/1‑150 (MQ=255)
                                                                           aTCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGgctgc                                                                            >  2:1591474/1‑150 (MQ=255)
                                                                            tCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGgctgct                                                                           >  2:599570/1‑150 (MQ=255)
                                                                               aTACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCtttt                                                                        >  1:3207315/1‑150 (MQ=255)
                                                                                  ctctTCTCGCCAGACGCGGTTTGCTCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTg                                                                     <  2:3079414/150‑1 (MQ=255)
                                                                                  ctctTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGAGACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTg                                                                     <  2:817485/150‑1 (MQ=255)
                                                                                    cttctCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTgcg                                                                   >  2:675992/1‑150 (MQ=255)
                                                                                     ttctCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTgcgc                                                                  >  2:860783/1‑150 (MQ=255)
                                                                                       ctcGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCtt                                                                <  1:779745/150‑1 (MQ=255)
                                                                                       ctcGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCtt                                                                <  2:4295514/150‑1 (MQ=255)
                                                                                        tcGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTa                                                               >  2:2030776/1‑150 (MQ=255)
                                                                                        tcGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTa                                                               >  2:1772249/1‑150 (MQ=255)
                                                                                           ccAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTAcgc                                                            >  2:4244257/1‑150 (MQ=255)
                                                                                             aGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCgg                                                          >  2:3893533/1‑150 (MQ=255)
                                                                                             aGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCgg                                                          >  1:438119/1‑150 (MQ=255)
                                                                                              gACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGga                                                         <  2:1410816/150‑1 (MQ=255)
                                                                                               aCGCGGTGTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGgac                                                        <  1:2836965/150‑1 (MQ=255)
                                                                                                  cgGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGgacgac                                                     <  2:2057747/150‑1 (MQ=255)
                                                                                                    gTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGc                                                   >  2:3833282/1‑150 (MQ=255)
                                                                                                       tGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGctct                                                <  1:2454209/150‑1 (MQ=255)
                                                                                                        gATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTa                                               >  1:4252979/1‑150 (MQ=255)
                                                                                                          tCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTAcg                                             <  1:906279/150‑1 (MQ=255)
                                                                                                           cGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTAcgc                                            >  2:3237848/1‑150 (MQ=255)
                                                                                                            gCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGCTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTAcgcg                                           >  2:616892/1‑150 (MQ=255)
                                                                                                                 aaGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCtt                                      <  1:1536601/150‑1 (MQ=255)
                                                                                                                 aaGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCtt                                      >  1:100965/1‑150 (MQ=255)
                                                                                                                     ggTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCttcttc                                  >  2:1536457/1‑150 (MQ=255)
                                                                                                                       tCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCa                                <  2:1840127/150‑1 (MQ=255)
                                                                                                                           gggggCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGt                            >  2:4093972/1‑150 (MQ=255)
                                                                                                                            ggggCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGtt                           >  1:1438475/1‑150 (MQ=255)
                                                                                                                                cTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCAt                       <  1:1796300/150‑1 (MQ=255)
                                                                                                                                 ttttAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATc                      >  1:4148238/1‑150 (MQ=255)
                                                                                                                                    tAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTc                   <  2:2133744/150‑1 (MQ=255)
                                                                                                                                     aGCGCAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCg                  >  2:512434/1‑150 (MQ=255)
                                                                                                                                      gcgcAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCgg                 <  1:3421169/150‑1 (MQ=255)
                                                                                                                                        gcAGAACCGATACAGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGAt               <  2:2320196/150‑1 (MQ=255)
                                                                                                                                          aGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTc             >  2:2852048/1‑150 (MQ=255)
                                                                                                                                           gAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCg            >  2:4264022/1‑150 (MQ=255)
                                                                                                                                            aaCCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGt           <  1:1560229/150‑1 (MQ=255)
                                                                                                                                               cgataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCaa        >  2:1718313/1‑150 (MQ=255)
                                                                                                                                                 ataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAAtt      <  1:1981950/150‑1 (MQ=255)
                                                                                                                                                 ataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAAtt      <  1:2618727/150‑1 (MQ=255)
                                                                                                                                                   accgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATttc    <  1:2495142/150‑1 (MQ=255)
                                                                                                                                                   accgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATttc    >  2:688233/1‑150 (MQ=255)
                                                                                                                                                    ccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATttct   <  2:4445225/150‑1 (MQ=255)
                                                                                                                                                    ccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATttct   >  2:237537/1‑150 (MQ=255)
                                                                                                                                                     cgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATttctt  <  1:2158306/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAGTATTTTCCGCAGACGTGGTTCCCGGCGCTGGCTGCACGCCATTCGCCTGATCGGTCGCATTGCCAGCAAGATCAATACTCTTCTCGCCAGACGCGGTTTGATCGCTGGAAGTGGTCGAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTT  >  CP009273/3510162‑3510460

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: