Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,482,005 G→T 100% D172Y (GAT→TAT)  ydcF → DUF218 superfamily protein, SAM‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,482,0050GT98.6% 962.0 / 2.4 296D172Y (GAT→TAT) ydcFDUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (148/143);  total (151/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 296 total aligned reads displayed.

TCAACAAACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACGC  >  CP009273/1481858‑1482151
                                                                                                                                                   |                                                                                                                                                  
tcatcaAACTGCGGTGAAAACGCACTCTTTAGCTTCGAGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCATGACCCCACCATGCAGCGGCGCACGATGTCGACGTTCCGCCGTATGACTGGTGACAATCCCTAt                                                                                                                                                  <  2:224475/150‑1 (MQ=255)
tCAACAAACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGAt                                                                                                                                                  <  1:4406929/150‑1 (MQ=255)
 caacaaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATg                                                                                                                                                 >  1:3697092/1‑150 (MQ=255)
 caacaaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATg                                                                                                                                                 >  2:26923/1‑150 (MQ=255)
  aacaaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGc                                                                                                                                                >  2:2106698/1‑150 (MQ=255)
    caaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcac                                                                                                                                              >  1:494615/1‑150 (MQ=255)
     aaaCTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcacc                                                                                                                                             >  1:3228111/1‑150 (MQ=255)
      aaCTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcacca                                                                                                                                            >  2:2553539/1‑150 (MQ=255)
       aCTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcaccac                                                                                                                                           <  1:4881170/150‑1 (MQ=255)
       aCTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcaccac                                                                                                                                           <  1:2275944/150‑1 (MQ=255)
          gCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCt                                                                                                                                        <  2:696789/150‑1 (MQ=255)
           cGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTg                                                                                                                                       >  2:2068162/1‑150 (MQ=255)
             gTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGt                                                                                                                                     <  1:1996012/150‑1 (MQ=255)
             gTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGt                                                                                                                                     <  1:1393363/150‑1 (MQ=255)
             gTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGt                                                                                                                                     <  2:1802849/150‑1 (MQ=255)
               gAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTa                                                                                                                                   >  1:1604313/1‑150 (MQ=255)
               gAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTa                                                                                                                                   >  1:3299058/1‑150 (MQ=255)
               gAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTa                                                                                                                                   >  1:4094645/1‑150 (MQ=255)
                 aaaCGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAg                                                                                                                                 >  1:1268980/1‑150 (MQ=255)
                 aaaCGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAg                                                                                                                                 <  2:3914694/150‑1 (MQ=255)
                 aaaCGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAg                                                                                                                                 >  2:4139945/1‑150 (MQ=255)
                    cgcacgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTa                                                                                                                              <  2:4916534/150‑1 (MQ=255)
                    cgcacgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTa                                                                                                                              <  1:2699317/150‑1 (MQ=255)
                        cgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATccc                                                                                                                          >  2:66198/1‑150 (MQ=255)
                        cgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATccc                                                                                                                          <  1:1206500/150‑1 (MQ=255)
                         gcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCg                                                                                                                         <  2:1413795/150‑1 (MQ=255)
                         gcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCg                                                                                                                         >  1:3559511/1‑150 (MQ=255)
                           tttAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGa                                                                                                                       >  1:2535910/1‑150 (MQ=255)
                             tAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGAtt                                                                                                                     >  2:710126/1‑150 (MQ=255)
                              aGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTc                                                                                                                    <  2:1243821/150‑1 (MQ=255)
                              aGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTc                                                                                                                    >  2:114913/1‑150 (MQ=255)
                               gCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCg                                                                                                                   <  1:3664894/150‑1 (MQ=255)
                                cATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGt                                                                                                                  >  2:588229/1‑150 (MQ=255)
                                   cgcgCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTcc                                                                                                               <  2:1743324/150‑1 (MQ=255)
                                    gcgcTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCt                                                                                                              >  1:4636983/1‑150 (MQ=255)
                                    gcgcTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCt                                                                                                              <  1:1864796/150‑1 (MQ=255)
                                       cTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTCTCCCGGATTCGTTCCTCAg                                                                                                           >  2:4711175/1‑150 (MQ=255)
                                       cTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCAGTATGACTGGGGACAATCCCTATGCACCACGCTTGTTAAGTTATCCCGGATTCGTTCCTCAg                                                                                                           <  1:2774862/150‑1 (MQ=255)
                                         aTTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGtt                                                                                                         >  2:2912176/1‑150 (MQ=255)
                                         aTTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGtt                                                                                                         <  1:2404017/150‑1 (MQ=255)
                                          ttGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTa                                                                                                        <  1:4854774/150‑1 (MQ=255)
                                           tGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAg                                                                                                       >  2:1860684/1‑150 (MQ=255)
                                              aTCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGaa                                                                                                    >  1:2189713/1‑150 (MQ=255)
                                                cAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAAt                                                                                                  >  1:2055916/1‑150 (MQ=255)
                                                 aGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATa                                                                                                 >  1:4596684/1‑150 (MQ=255)
                                                  ggCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATaa                                                                                                >  2:4203082/1‑150 (MQ=255)
                                                     cGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGc                                                                                             >  2:5153634/1‑150 (MQ=255)
                                                        aGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGa                                                                                          >  1:4281806/1‑150 (MQ=255)
                                                          aaCGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACa                                                                                        >  2:3540996/1‑150 (MQ=255)
                                                            cGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgt                                                                                      <  1:3073619/150‑1 (MQ=255)
                                                            cGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgt                                                                                      <  2:124691/150‑1 (MQ=255)
                                                            cGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgt                                                                                      <  1:3862247/150‑1 (MQ=255)
                                                             gAGTTCAGACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgtg                                                                                     <  2:3455294/150‑1 (MQ=255)
                                                                ttCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTaa                                                                                  <  2:2802482/150‑1 (MQ=255)
                                                                ttCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTaa                                                                                  >  2:4000414/1‑150 (MQ=255)
                                                                   aTACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCt                                                                               <  1:749377/150‑1 (MQ=255)
                                                                   aTACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCt                                                                               >  1:3998001/1‑150 (MQ=255)
                                                                   aTACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCt                                                                               >  2:368141/1‑150 (MQ=255)
                                                                    tACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCtt                                                                              <  2:3415773/150‑1 (MQ=255)
                                                                      cGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTtta                                                                            >  1:469244/1‑150 (MQ=255)
                                                                            tCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATc                                                                      >  1:4533738/1‑150 (MQ=255)
                                                                             cGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCa                                                                     <  2:1696185/150‑1 (MQ=255)
                                                                                tgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGtt                                                                  >  1:1096144/1‑150 (MQ=255)
                                                                                tgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGtt                                                                  <  1:3620479/150‑1 (MQ=255)
                                                                                 gttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTa                                                                 <  2:2813101/150‑1 (MQ=255)
                                                                                 gttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTa                                                                 <  1:65997/150‑1 (MQ=255)
                                                                                    cAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACaa                                                              >  2:4289632/1‑150 (MQ=255)
                                                                                      ggACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAgg                                                            >  1:3504069/1‑150 (MQ=255)
                                                                                       gACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGa                                                           <  1:3002595/150‑1 (MQ=255)
                                                                                       gACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGa                                                           >  1:5126632/1‑150 (MQ=255)
                                                                                          cccACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTa                                                        <  2:2034340/150‑1 (MQ=255)
                                                                                          cccACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTa                                                        <  1:41704/150‑1 (MQ=255)
                                                                                           ccaccaTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTAt                                                       <  2:1639845/150‑1 (MQ=255)
                                                                                                aTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCa                                                  >  1:3833576/1‑150 (MQ=255)
                                                                                                  gCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGt                                                >  2:856678/1‑150 (MQ=255)
                                                                                                  gCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGt                                                >  2:4521368/1‑150 (MQ=255)
                                                                                                    aGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTg                                              >  2:3374073/1‑150 (MQ=255)
                                                                                                    aGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTg                                              >  2:3831935/1‑150 (MQ=255)
                                                                                                       ggcgCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGc                                           <  1:1108737/150‑1 (MQ=255)
                                                                                                          gcACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGtt                                        >  2:243494/1‑150 (MQ=255)
                                                                                                           cACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTa                                       <  1:2350478/150‑1 (MQ=255)
                                                                                                           cACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTa                                       >  1:3808583/1‑150 (MQ=255)
                                                                                                           cACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTa                                       >  2:4432159/1‑150 (MQ=255)
                                                                                                                 ggCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTAtctctc                                 >  1:1481810/1‑150 (MQ=255)
                                                                                                                      cGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCactac                            <  1:2220494/150‑1 (MQ=255)
                                                                                                                         tCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCa                         <  1:2832327/150‑1 (MQ=255)
                                                                                                                          ccgccgTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCAc                        >  2:3202777/1‑150 (MQ=255)
                                                                                                                             ccgTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTgg                     <  2:56464/150‑1 (MQ=255)
                                                                                                                               gTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCg                   >  1:4819697/1‑150 (MQ=255)
                                                                                                                                tATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGa                  <  2:2825214/150‑1 (MQ=255)
                                                                                                                                   gACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCt               <  1:2584150/150‑1 (MQ=255)
                                                                                                                                      tGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGcc            >  2:3632552/1‑150 (MQ=255)
                                                                                                                                       ggggACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACACTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCcg           >  2:4778962/1‑150 (MQ=255)
                                                                                                                                        gggACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCcgc          <  2:2515725/150‑1 (MQ=255)
                                                                                                                                          gACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGt        >  2:1913455/1‑150 (MQ=255)
                                                                                                                                            cAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGttt      <  1:564731/150‑1 (MQ=255)
                                                                                                                                            cAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGttt      >  1:1765147/1‑150 (MQ=255)
                                                                                                                                             aaTCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTa     <  2:1432663/150‑1 (MQ=255)
                                                                                                                                              aTCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTAc    >  1:99759/1‑150 (MQ=255)
                                                                                                                                                cccTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTAcgc  <  2:4616613/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCAACAAACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACGC  >  CP009273/1481858‑1482151

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: