Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 278,630 G→T 5.5% A214A (GCG→GCT yagE → 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273278,6300GT 5.5% 1499.6 / 39.2 391A214A (GCG→GCTyagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (194/175);  new base T (13/9);  total (207/184)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 397 total aligned reads displayed.

GGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGA  >  CP009273/278481‑278779
                                                                                                                                                     |                                                                                                                                                     
ggCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCg                                                                                                                                                       <  2:2156327/150‑1 (MQ=255)
 gCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGa                                                                                                                                                      <  1:3402118/150‑1 (MQ=255)
  cATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAg                                                                                                                                                     <  2:4513607/150‑1 (MQ=255)
    tCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCg                                                                                                                                                   >  1:1369631/1‑150 (MQ=255)
    tCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCg                                                                                                                                                   <  1:2615431/150‑1 (MQ=255)
     cAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCTAGCgg                                                                                                                                                  >  2:4110914/1‑150 (MQ=255)
     cAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCgg                                                                                                                                                  >  1:4064289/1‑150 (MQ=255)
     cAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCgg                                                                                                                                                  >  1:3667664/1‑150 (MQ=255)
      aaaGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGc                                                                                                                                                 >  1:5533659/1‑150 (MQ=255)
      aaaGACACCATCCACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCACAGGTGCCCAGCCGCACTTCACAGTGCACTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGc                                                                                                                                                 >  2:1775055/1‑150 (MQ=255)
          acacCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACt                                                                                                                                             >  1:5011562/1‑150 (MQ=255)
          acacCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACt                                                                                                                                             >  1:2451344/1‑150 (MQ=255)
            acCCTCGACTCCGTCGCCCACCGGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCGGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACttt                                                                                                                                           <  1:4318957/150‑1 (MQ=255)
            acCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACttt                                                                                                                                           <  1:4006841/150‑1 (MQ=255)
                 cGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGcccc                                                                                                                                      <  2:1839212/150‑1 (MQ=255)
                  gACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCg                                                                                                                                     >  1:2298796/1‑150 (MQ=255)
                  gACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCg                                                                                                                                     >  2:5266268/1‑150 (MQ=255)
                  gACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCg                                                                                                                                     >  1:3215876/1‑150 (MQ=255)
                   aCTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGGTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCTAGCGGCAACTTTGCCCCGc                                                                                                                                    >  1:2388525/1‑150 (MQ=255)
                    cTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCa                                                                                                                                   <  1:3416550/150‑1 (MQ=255)
                     tCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAg                                                                                                                                  >  2:1053103/1‑150 (MQ=255)
                        gTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCTAGCGGCAACTTTGCCCCGCAGgtg                                                                                                                               <  1:1020587/150‑1 (MQ=255)
                          cGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTc                                                                                                                             <  1:1220886/150‑1 (MQ=255)
                           gCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCg                                                                                                                            <  2:1241496/150‑1 (MQ=255)
                             ccACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGt                                                                                                                          >  2:5065539/1‑150 (MQ=255)
                             ccACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGt                                                                                                                          >  1:829217/1‑150 (MQ=255)
                              cACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTg                                                                                                                         >  1:4891164/1‑150 (MQ=255)
                               aCCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGa                                                                                                                        >  2:1652865/1‑150 (MQ=255)
                               aCCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGa                                                                                                                        >  2:3698819/1‑150 (MQ=255)
                                ccTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGGGaa                                                                                                                       >  1:2954010/1‑150 (MQ=255)
                                      cagcgTGATCCATCCCGTCCAAGGTGCCCCTCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAAtcttct                                                                                                                 <  2:3055661/150‑1 (MQ=255)
                                       agcaTGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTg                                                                                                                <  2:2843338/150‑1 (MQ=255)
                                         caTGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGaa                                                                                                              <  1:4676117/150‑1 (MQ=255)
                                            gATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCTAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGc                                                                                                           >  2:3861878/1‑150 (MQ=255)
                                            gATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGc                                                                                                           >  1:1962138/1‑150 (MQ=255)
                                             aTCCAGACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGcc                                                                                                          <  1:4113258/150‑1 (MQ=255)
                                              tCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCt                                                                                                         <  1:2079868/150‑1 (MQ=255)
                                              tCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCt                                                                                                         <  1:3507383/150‑1 (MQ=255)
                                                cATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTgg                                                                                                       <  1:5639198/150‑1 (MQ=255)
                                                 aTACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGgc                                                                                                      <  1:4184257/150‑1 (MQ=255)
                                                      gTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGAc                                                                                                 <  1:1620047/150‑1 (MQ=255)
                                                      gTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGAc                                                                                                 <  1:4862054/150‑1 (MQ=255)
                                                      gTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGAc                                                                                                 >  1:3560426/1‑150 (MQ=255)
                                                             gTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACg                                                                                          >  1:1444216/1‑150 (MQ=255)
                                                               gCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTg                                                                                        >  1:3342132/1‑150 (MQ=255)
                                                                  cATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCg                                                                                     <  2:2778952/150‑1 (MQ=255)
                                                                     ccGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGCCGTGGCGaaa                                                                                  >  1:5218887/1‑150 (MQ=255)
                                                                     ccGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGaaa                                                                                  >  1:3253790/1‑150 (MQ=255)
                                                                         aCTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCgg                                                                              >  2:4330812/1‑150 (MQ=255)
                                                                           ttCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGcc                                                                            <  2:2996136/150‑1 (MQ=255)
                                                                                cGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTa                                                                       <  2:3404807/150‑1 (MQ=255)
                                                                                 gTGCTCTGCGGCTACGACGATCATCTGTTCAATCCCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTat                                                                      <  1:3756244/150‑1 (MQ=255)
                                                                                  tGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTatc                                                                     <  2:1907873/150‑1 (MQ=255)
                                                                                    ctctGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTatcat                                                                   >  2:2710224/1‑150 (MQ=255)
                                                                                    ctctGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTatcat                                                                   <  2:2515408/150‑1 (MQ=255)
                                                                                    ctctGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTatcat                                                                   <  2:636404/150‑1 (MQ=255)
                                                                                     tctGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTatcatc                                                                  <  2:4059335/150‑1 (MQ=255)
                                                                                       tGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAg                                                                >  2:1477632/1‑150 (MQ=255)
                                                                                       tGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAg                                                                >  2:4280637/1‑150 (MQ=255)
                                                                                         cGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGAc                                                              <  2:862081/150‑1 (MQ=255)
                                                                                           gCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCt                                                            >  1:876523/1‑150 (MQ=255)
                                                                                           gCTACGACGATCATCTGTTCAATACCCGGCGGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCt                                                            <  1:1497015/150‑1 (MQ=255)
                                                                                            cTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCtt                                                           >  1:2779819/1‑150 (MQ=255)
                                                                                             tACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTTAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTtg                                                          <  1:735621/150‑1 (MQ=255)
                                                                                             tACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTGTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTtg                                                          <  2:1162364/150‑1 (MQ=255)
                                                                                                 acgaTCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTtgctgc                                                      >  2:1962749/1‑150 (MQ=255)
                                                                                                  cgaTCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCa                                                     >  2:4219386/1‑150 (MQ=255)
                                                                                                     tcatcTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAAt                                                  >  2:4861408/1‑150 (MQ=255)
                                                                                                     tcatcTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAAt                                                  <  1:1863801/150‑1 (MQ=255)
                                                                                                       atcTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTc                                                >  1:5405474/1‑150 (MQ=255)
                                                                                                            ttCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAg                                           >  1:219431/1‑150 (MQ=255)
                                                                                                                aTACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGt                                       >  1:1134879/1‑150 (MQ=255)
                                                                                                                 tACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTa                                      <  2:2229306/150‑1 (MQ=255)
                                                                                                                 tACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTa                                      <  2:4768515/150‑1 (MQ=255)
                                                                                                                 tACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTa                                      <  2:4598845/150‑1 (MQ=255)
                                                                                                                  aCCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGtct                                     >  2:3808481/1‑148 (MQ=255)
                                                                                                                   cccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATc                                    >  2:3532644/1‑150 (MQ=255)
                                                                                                                   cccTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATc                                    <  2:1293028/150‑1 (MQ=255)
                                                                                                                        ctgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTg                               <  2:3215876/150‑1 (MQ=255)
                                                                                                                         tgctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTgg                              <  2:5559976/150‑1 (MQ=255)
                                                                                                                         tgctCGGCGGAGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTgg                              <  1:1949856/150‑1 (MQ=255)
                                                                                                                         tgctCGGCGGAGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTgg                              <  1:1049923/150‑1 (MQ=255)
                                                                                                                           ctCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGAt                            <  2:5324647/150‑1 (MQ=255)
                                                                                                                            tCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCACATTCCGCAGATGTATCAGCTGGATa                           >  2:1820625/1‑150 (MQ=255)
                                                                                                                            tCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATa                           <  2:3616989/150‑1 (MQ=255)
                                                                                                                                cggcgACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGcc                       <  1:1733125/150‑1 (MQ=255)
                                                                                                                                     aCGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTtg                  >  2:4410337/1‑150 (MQ=255)
                                                                                                                                       ggggCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCACATTCCGCAGATGTATCAGCTGGATACGCCGTTtgtg                >  2:4822392/1‑150 (MQ=255)
                                                                                                                                       ggggCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTtgtg                >  1:4421413/1‑150 (MQ=255)
                                                                                                                                       ggggCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTtgtg                >  1:45568/1‑150 (MQ=255)
                                                                                                                                        gggCGATATCGGCTAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGa               <  1:1901586/150‑1 (MQ=255)
                                                                                                                                         ggCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGTTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGaa              <  2:5637987/150‑1 (MQ=255)
                                                                                                                                          gCGATATCGGCGAGCGGCAACTGTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGCCGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAAc             <  2:1596566/150‑1 (MQ=255)
                                                                                                                                           cGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACg            >  2:2341645/1‑150 (MQ=255)
                                                                                                                                            gATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGt           <  2:3667418/150‑1 (MQ=255)
                                                                                                                                             atatCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTg          <  1:4403733/150‑1 (MQ=255)
                                                                                                                                                  ggCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTaa     >  1:3862776/1‑150 (MQ=255)
                                                                                                                                                  ggCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTaa     <  1:602899/150‑1 (MQ=255)
                                                                                                                                                    cGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAag   >  1:5439590/1‑150 (MQ=255)
                                                                                                                                                     gAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAaga  >  1:4754125/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCATCAAAGACACCATCGACTCCGTCGCCCACCTGCGCAGCATGATCCATACCGTCAAAGGTGCCCATCCGCACTTCACCGTGCTCTGCGGCTACGACGATCATCTGTTCAATACCCTGCTGCTCGGCGGCGACGGGGCGATATCGGCGAGCGGCAACTTTGCCCCGCAGGTGTCGGTGAATCTTCTGAAAGCCTGGCGCGACGGGGACGTGGCGAAAGCGGCCGGGTATCATCAGACCTTGCTGCAAATTCCGCAGATGTATCAGCTGGATACGCCGTTTGTGAACGTGATTAAAGA  >  CP009273/278481‑278779

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: