Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,340,820 G→T 39.8% I60I (ATC→ATA gmr ← cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,340,8200GT39.8% 307.5 / 279.4 246I60I (ATC→ATAgmrcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base G (76/72);  new base T (54/44);  total (130/116)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

CCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTT  >  CP009273/1340673‑1340966
                                                                                                                                                   |                                                                                                                                                  
 ccGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGttatg                                                                                                                                                 >  2:10233/1‑150 (MQ=255)
 ccGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtgatg                                                                                                                                                 <  1:2188044/150‑1 (MQ=255)
  cGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtgatga                                                                                                                                                <  2:1107302/150‑1 (MQ=255)
    tgtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGAcc                                                                                                                                              <  1:1244105/150‑1 (MQ=255)
     gtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGAcct                                                                                                                                             <  2:216844/150‑2 (MQ=255)
     gtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCg                                                                                                                                             >  2:1456824/1‑150 (MQ=255)
          aaaTCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAtt                                                                                                                                        >  1:3913359/1‑150 (MQ=255)
            aTCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTCCCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTc                                                                                                                                      >  1:2072804/1‑150 (MQ=255)
            aTCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTc                                                                                                                                      >  2:3067376/1‑150 (MQ=255)
             tCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTcc                                                                                                                                     >  1:525396/1‑150 (MQ=255)
               aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCt                                                                                                                                   >  1:3046261/1‑150 (MQ=255)
               aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCt                                                                                                                                   <  1:538591/150‑1 (MQ=255)
               aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCt                                                                                                                                   <  1:2905231/150‑1 (MQ=255)
               aCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGACCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCt                                                                                                                                   <  1:3654518/150‑1 (MQ=255)
                 gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGa                                                                                                                                 >  2:3842463/1‑150 (MQ=255)
                  aGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGAt                                                                                                                                <  2:3184279/150‑1 (MQ=255)
                   gCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATg                                                                                                                               <  1:101421/150‑1 (MQ=255)
                      aCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGa                                                                                                                            <  2:1929887/150‑1 (MQ=255)
                        cgccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGagc                                                                                                                          >  2:1220512/1‑150 (MQ=255)
                         gccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGTTGATGACCGTCATTTCCCTGATGCGagca                                                                                                                         >  2:2059289/1‑150 (MQ=255)
                         gccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGagca                                                                                                                         <  1:2811171/150‑1 (MQ=255)
                           cgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGagcagc                                                                                                                       >  2:4057669/1‑150 (MQ=255)
                             gtgtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCt                                                                                                                     >  1:254799/1‑150 (MQ=255)
                                 gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTc                                                                                                                 >  1:2034009/1‑150 (MQ=255)
                                   catgccatgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCgg                                                                                                               >  2:409621/1‑150 (MQ=255)
                                          tgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATaa                                                                                                        >  1:2888120/1‑150 (MQ=255)
                                          tgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATaa                                                                                                        >  2:813813/1‑150 (MQ=255)
                                           gcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAAc                                                                                                       >  1:3367070/1‑150 (MQ=255)
                                           gcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAAc                                                                                                       <  1:3650673/150‑1 (MQ=255)
                                               aCGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCa                                                                                                   <  2:2683507/150‑1 (MQ=255)
                                               aCGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCa                                                                                                   <  2:1225877/150‑1 (MQ=255)
                                                cgcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCaa                                                                                                  <  1:4035925/150‑1 (MQ=255)
                                                cgcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCaa                                                                                                  >  1:4082769/1‑150 (MQ=255)
                                                     tGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTAc                                                                                             >  1:1825272/1‑150 (MQ=255)
                                                      gCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACg                                                                                            >  2:400354/1‑150 (MQ=255)
                                                      gCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACg                                                                                            <  2:1713498/150‑1 (MQ=255)
                                                       ccAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGa                                                                                           >  2:642829/1‑150 (MQ=255)
                                                        cAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGAt                                                                                          <  2:498449/150‑1 (MQ=255)
                                                          gCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCa                                                                                        <  2:1188477/150‑1 (MQ=255)
                                                               ttCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTc                                                                                   <  1:2300049/150‑1 (MQ=255)
                                                               ttCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTc                                                                                   <  1:1817444/150‑1 (MQ=255)
                                                                tCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCt                                                                                  <  1:1560983/150‑1 (MQ=255)
                                                                 cccGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCtg                                                                                 <  1:1698478/150‑1 (MQ=255)
                                                                   cGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCtgtg                                                                               <  2:749086/150‑1 (MQ=255)
                                                                   cGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCtgtg                                                                               >  1:1801114/1‑150 (MQ=255)
                                                                    gTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTgg                                                                              >  2:2290136/1‑150 (MQ=255)
                                                                     tttATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGt                                                                             >  1:137878/1‑150 (MQ=255)
                                                                      ttATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGtt                                                                            >  2:1345196/1‑150 (MQ=255)
                                                                      ttATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGtt                                                                            >  1:830338/1‑150 (MQ=255)
                                                                      ttATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGtt                                                                            >  2:219945/1‑150 (MQ=255)
                                                                       tATTGATTTTTCGTCCTACCAGATGCACAGTAAGATAGATTTCATCGAAGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCATGATGCGAGCAGCCTGTTCGGCAGATATCTGCTATGTACTATCAGGTTCTGTGGttt                                                                           <  1:3633327/150‑1 (MQ=255)
                                                                        aTTGATTTTTCGTCCTACCAGATGCACAGANAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTc                                                                          >  1:3853290/1‑150 (MQ=255)
                                                                        aTTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTc                                                                          <  2:2857613/150‑1 (MQ=255)
                                                                        aTTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTc                                                                          <  2:2467473/150‑1 (MQ=255)
                                                                         ttGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGGGGTTTCg                                                                         >  2:1513528/1‑150 (MQ=255)
                                                                            aTTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCggtg                                                                      <  1:2115368/150‑1 (MQ=255)
                                                                            aTTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCggtg                                                                      >  2:480053/1‑150 (MQ=255)
                                                                             tttttCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCggtgg                                                                     >  1:3114949/1‑150 (MQ=255)
                                                                              ttttCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGTTTTCggtggt                                                                    >  2:2037983/1‑150 (MQ=255)
                                                                                 tCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGa                                                                 <  1:2729749/150‑1 (MQ=255)
                                                                                  cGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGaa                                                                >  2:3500055/1‑150 (MQ=255)
                                                                                  cGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGaa                                                                >  1:4147986/1‑150 (MQ=255)
                                                                                     ccTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAg                                                             >  1:2351653/1‑150 (MQ=255)
                                                                                      cTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAgc                                                            <  2:1114683/150‑1 (MQ=255)
                                                                                       tACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAgcg                                                           <  1:2003637/150‑1 (MQ=255)
                                                                                          cAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCaa                                                        >  2:3399144/1‑150 (MQ=255)
                                                                                            gATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCaaaa                                                      >  2:2488534/1‑150 (MQ=255)
                                                                                               gaacaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGATGAAAAGCGCAAAACAt                                                   >  2:1060049/3‑150 (MQ=255)
                                                                                               gCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACAt                                                   <  1:1775486/150‑1 (MQ=255)
                                                                                                cacaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATc                                                  <  2:1656467/150‑1 (MQ=255)
                                                                                                 acaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCg                                                 <  2:3297418/150‑1 (MQ=255)
                                                                                                   aGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgct                                               >  2:2409965/1‑150 (MQ=255)
                                                                                                    gAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCgctg                                              <  2:3264034/150‑1 (MQ=255)
                                                                                                     aaaGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctgc                                             >  2:891357/1‑150 (MQ=255)
                                                                                                      aaGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctgct                                            >  2:3875506/1‑150 (MQ=255)
                                                                                                      aaGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctgct                                            >  2:364310/1‑150 (MQ=255)
                                                                                                          tCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTcc                                        >  2:3546952/1‑150 (MQ=255)
                                                                                                           cGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCg                                       <  2:2686573/150‑1 (MQ=255)
                                                                                                           cGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCg                                       <  1:1178733/150‑1 (MQ=255)
                                                                                                            gCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGt                                      >  1:2946072/1‑150 (MQ=255)
                                                                                                             cTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTc                                     <  1:1754665/150‑1 (MQ=255)
                                                                                                              tttCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCa                                    <  1:4135204/150‑1 (MQ=255)
                                                                                                               ttCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCaa                                   <  1:679364/150‑1 (MQ=255)
                                                                                                               ttCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCaa                                   >  2:341265/1‑150 (MQ=255)
                                                                                                                tCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAAc                                  >  2:1903803/1‑150 (MQ=255)
                                                                                                                  aTCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCg                                <  1:10233/150‑1 (MQ=255)
                                                                                                                  aTCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCg                                >  2:558565/1‑150 (MQ=255)
                                                                                                                  aTCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCg                                <  2:3010552/150‑1 (MQ=255)
                                                                                                                           cAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATgg                       <  1:2630602/150‑1 (MQ=255)
                                                                                                                            aGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGa                      >  2:7592/1‑150 (MQ=255)
                                                                                                                            aGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGa                      <  2:3231527/150‑1 (MQ=255)
                                                                                                                              aCTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAtt                    <  1:2872151/150‑1 (MQ=255)
                                                                                                                                    catcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGa              >  1:1229261/1‑150 (MQ=255)
                                                                                                                                     atcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAt             >  1:3286623/1‑150 (MQ=255)
                                                                                                                                     atcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAt             <  1:3337852/150‑1 (MQ=255)
                                                                                                                                      tcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGCAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATc            <  2:2086563/150‑1 (MQ=255)
                                                                                                                                       caGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATcc           >  2:1982548/1‑150 (MQ=255)
                                                                                                                                        aGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCg          <  1:716126/150‑1 (MQ=255)
                                                                                                                                        aGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCg          >  2:614634/1‑150 (MQ=255)
                                                                                                                                             ggAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGaaa     <  1:3597072/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTT  >  CP009273/1340673‑1340966

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: