Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,482,005 G→T 87.7% 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,0050GT87.7% 873.2 / 129.9 349D172Y (GAT→TAT) ydcFDUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base G (29/14);  new base T (139/167);  total (168/181)
Fisher's exact test for biased strand distribution p-value = 8.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 349 total aligned reads displayed.

AACAAACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACG  >  CP009273/1481860‑1482150
                                                                                                                                                 |                                                                                                                                                 
 acaaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGc                                                                                                                                               >  1:393021/1‑149 (MQ=255)
 acaaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGc                                                                                                                                               >  1:3803830/1‑149 (MQ=255)
  caaACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGca                                                                                                                                              >  2:298317/1‑149 (MQ=255)
    aaCTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcacc                                                                                                                                            >  1:903508/1‑149 (MQ=255)
     aCTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGcacca                                                                                                                                           >  2:647896/1‑149 (MQ=255)
          ggTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTg                                                                                                                                      <  2:1470456/149‑1 (MQ=255)
             gAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGtt                                                                                                                                   <  1:41775/149‑1 (MQ=255)
               aaaCGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTaa                                                                                                                                 >  1:862791/1‑149 (MQ=255)
               aaaCGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTaa                                                                                                                                 >  1:3569143/1‑149 (MQ=255)
                  cgcacgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGACAATCCCTATGCACCACGCTGGTTAAGTTa                                                                                                                             <  2:2862938/149‑1 (MQ=255)
                   gcacgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTa                                                                                                                             >  2:3953489/1‑149 (MQ=255)
                     acgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATc                                                                                                                           <  1:4739130/149‑1 (MQ=255)
                     acgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATc                                                                                                                           >  1:3389560/1‑149 (MQ=255)
                      cgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATcc                                                                                                                          <  2:906238/149‑1 (MQ=255)
                      cgcTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATcc                                                                                                                          >  1:4147514/1‑149 (MQ=255)
                         tttAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCgg                                                                                                                       >  1:2497555/1‑149 (MQ=255)
                          ttAGCATCGCGCTATTGAATCAGGCCGTAGAAAGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATTCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCACgta                                                                                                                      >  2:190021/1‑147 (MQ=255)
                            aGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGAtt                                                                                                                    <  2:3239066/149‑1 (MQ=255)
                               aTCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGt                                                                                                                 >  1:797413/1‑149 (MQ=255)
                                  gcgcTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTcc                                                                                                              <  1:1606141/149‑1 (MQ=255)
                                     cTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCa                                                                                                           <  1:4593640/149‑1 (MQ=255)
                                      tATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCAGGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACACTCCCTATGCACCTCGCTGGTGAAGTTATCCCGGAGTCGTTCCTCAg                                                                                                          >  2:2879329/1‑149 (MQ=255)
                                         tGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTa                                                                                                       <  1:2860351/149‑1 (MQ=255)
                                           aaTCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAgg                                                                                                     <  1:1034665/149‑1 (MQ=255)
                                              cAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGaaa                                                                                                  <  1:2865045/149‑1 (MQ=255)
                                                ggCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAACTa                                                                                                >  2:4592687/1‑149 (MQ=255)
                                                   cGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACg                                                                                             >  1:357889/1‑149 (MQ=255)
                                                   cGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACg                                                                                             >  1:3667841/1‑149 (MQ=255)
                                                     tAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCa                                                                                           <  1:3533652/149‑1 (MQ=255)
                                                     tAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCa                                                                                           >  1:1522812/1‑149 (MQ=255)
                                                       gAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGa                                                                                         <  2:404870/149‑1 (MQ=255)
                                                       gAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGa                                                                                         <  1:3901677/149‑1 (MQ=255)
                                                        aaCGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGAc                                                                                        <  2:2501596/149‑1 (MQ=255)
                                                          cGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAg                                                                                      >  2:2674958/1‑149 (MQ=255)
                                                          cGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAg                                                                                      <  1:3844044/149‑1 (MQ=255)
                                                            aggtCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTCTGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgtg                                                                                    >  1:1383603/4‑149 (MQ=255)
                                                             gTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgtgt                                                                                   >  1:717834/1‑149 (MQ=255)
                                                             gTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAgtgt                                                                                   <  2:4050045/149‑1 (MQ=255)
                                                              ttCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTa                                                                                  <  1:2166498/149‑1 (MQ=255)
                                                                cATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAAt                                                                                >  1:974648/1‑149 (MQ=255)
                                                                   aCGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCtt                                                                             <  2:392718/149‑1 (MQ=255)
                                                                   aCGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCtt                                                                             <  1:4154848/149‑1 (MQ=255)
                                                                     ggCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTtta                                                                           >  2:794032/1‑149 (MQ=255)
                                                                     ggCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTtta                                                                           >  2:3739620/1‑149 (MQ=255)
                                                                      gCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTttat                                                                          >  2:2143161/1‑149 (MQ=255)
                                                                      gCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTttat                                                                          >  1:3997028/1‑149 (MQ=255)
                                                                      gCTATCGTTGGTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTttat                                                                          <  2:4709575/149‑1 (MQ=255)
                                                                       cTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTttatt                                                                         >  1:2935099/1‑149 (MQ=255)
                                                                         aTCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTaa                                                                       >  1:2665494/1‑149 (MQ=255)
                                                                         aTCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTaa                                                                       >  1:1863467/1‑149 (MQ=255)
                                                                         aTCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTaa                                                                       <  2:3177697/149‑1 (MQ=255)
                                                                         aTCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTaa                                                                       >  2:2633830/1‑149 (MQ=255)
                                                                         aTCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTaa                                                                       <  2:1172278/149‑1 (MQ=255)
                                                                            gttgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCa                                                                    <  1:4562293/149‑1 (MQ=255)
                                                                            gttgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCa                                                                    <  1:4630847/149‑1 (MQ=255)
                                                                            gttgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCa                                                                    >  2:1728489/1‑149 (MQ=255)
                                                                             ttgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTCGGAACTAACGCAGACAGTGTAATCTTTATTCATCAg                                                                   >  2:3313632/1‑149 (MQ=255)
                                                                             ttgttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAg                                                                   >  1:2805588/1‑149 (MQ=255)
                                                                                ttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTa                                                                <  2:2366714/149‑1 (MQ=255)
                                                                                ttCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTa                                                                >  1:1000020/1‑149 (MQ=255)
                                                                                    ggACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAg                                                            >  1:4210991/1‑149 (MQ=255)
                                                                                     gACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAgg                                                           <  2:1801831/149‑1 (MQ=255)
                                                                                       ccccACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGAt                                                         >  1:2097784/1‑149 (MQ=255)
                                                                                       ccccACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGAt                                                         <  1:1257707/149‑1 (MQ=255)
                                                                                       ccccACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGAt                                                         >  1:1669383/1‑149 (MQ=255)
                                                                                       ccccACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGAt                                                         <  1:1006050/149‑1 (MQ=255)
                                                                                         ccaccaTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTa                                                       <  1:1702947/149‑1 (MQ=255)
                                                                                         ccaccaTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTa                                                       >  1:3015654/1‑149 (MQ=255)
                                                                                          caccaTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTAt                                                      <  2:421630/149‑1 (MQ=255)
                                                                                           accaTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATAATg                                                     >  2:466529/1‑149 (MQ=255)
                                                                                             caTGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGc                                                   <  2:55771/149‑1 (MQ=255)
                                                                                                 cAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGt                                               <  1:1192898/149‑1 (MQ=255)
                                                                                                  aGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGtt                                              <  2:4285617/149‑1 (MQ=255)
                                                                                                   gcggcgCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTg                                             <  2:1651060/149‑1 (MQ=255)
                                                                                                    cggcgCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGa                                            <  2:2090921/149‑1 (MQ=255)
                                                                                                     ggcgCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTATGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAg                                           <  1:1086091/149‑1 (MQ=255)
                                                                                                     ggcgCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAg                                           <  2:1255242/149‑1 (MQ=255)
                                                                                                      gcgcACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGc                                          <  1:3726578/149‑1 (MQ=255)
                                                                                                         cACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGtt                                       >  1:2754477/1‑149 (MQ=255)
                                                                                                           cGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTAt                                     <  1:2039458/149‑1 (MQ=255)
                                                                                                            gATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTCATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTAtc                                    >  2:4480158/1‑149 (MQ=255)
                                                                                                               ggCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTAtctct                                 <  1:3347427/149‑1 (MQ=255)
                                                                                                                 cGACGTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTACTCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCa                               <  2:3662196/149‑1 (MQ=255)
                                                                                                                     gTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCactac                           >  2:3081263/1‑149 (MQ=255)
                                                                                                                     gTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCactac                           >  1:4541135/1‑149 (MQ=255)
                                                                                                                     gTTCCGCCGTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCactac                           >  2:2424228/1‑149 (MQ=255)
                                                                                                                        ccgccgTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCa                        <  1:4245339/149‑1 (MQ=255)
                                                                                                                         cgccgTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCAc                       >  2:2252957/1‑149 (MQ=255)
                                                                                                                             gTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGc                   >  2:3869158/1‑149 (MQ=255)
                                                                                                                             gTATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGc                   >  2:4366837/1‑149 (MQ=255)
                                                                                                                              tttgtcTGTGGACAATCCCTAAGCACCACGCTGGTTAAGTTATCCCGGATTCGTTACTAAGTTAGGATATTACGCAGACAGTGTAATCTTTATTATTAAGTTACAAGGATTATGGCCAGTTGTGCGTTATCTCTCACTACTCACTGGCg                  <  1:2267533/144‑1 (MQ=255)
                                                                                                                              tATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCg                  >  2:3380088/1‑149 (MQ=255)
                                                                                                                              tATGACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCg                  >  1:2371159/1‑149 (MQ=255)
                                                                                                                                 gACTGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGc               <  1:3516452/149‑1 (MQ=255)
                                                                                                                                    tGGGGACAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGc            >  2:3500258/1‑149 (MQ=255)
                                                                                                                                       ggACAATCCCGATGCACCCCGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGTGCGTTATCTCTCACTACTCACTGGCGAGCTGCcgc         >  2:2924438/1‑149 (MQ=255)
                                                                                                                                          cAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGtt      <  2:4114728/149‑1 (MQ=255)
                                                                                                                                          cAATCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGtt      >  2:1609333/1‑149 (MQ=255)
                                                                                                                                            aTCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTa    <  2:1323552/149‑1 (MQ=255)
                                                                                                                                             tCCCTATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCTCTGGCGAGCTGCCGCGTTTAc   >  2:3580798/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AACAAACTGCGGTGAAAACGCACGCTTTAGCATCGCGCTATTGAATCAGGCCGTAGAACGAGTTCATACGGCTATCGTTGTTCAGGACCCCACCATGCAGCGGCGCACGATGGCGACGTTCCGCCGTATGACTGGGGACAATCCCGATGCACCACGCTGGTTAAGTTATCCCGGATTCGTTCCTCAGTTAGGAAATAACGCAGACAGTGTAATCTTTATTAATCAGTTACAAGGATTATGGCCAGTTGAGCGTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACG  >  CP009273/1481860‑1482150

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: