Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,152,523 C→A 6.3% A661A (GCG→GCT ygiQ ← Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,152,5230CA 6.3% 978.3 / 24.5 254A661A (GCG→GCTygiQRadical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (104/134);  new base A (10/6);  total (114/140)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 261 total aligned reads displayed.

GGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTA  >  CP009273/3152374‑3152672
                                                                                                                                                     |                                                                                                                                                     
ggCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCATTCAGTTGCTTTTTGCCAATc                                                                                                                                                       <  2:3278257/150‑1 (MQ=255)
ggCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATc                                                                                                                                                       <  2:1445382/150‑1 (MQ=255)
  cGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGc                                                                                                                                                     >  2:385561/1‑150 (MQ=255)
   gTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCt                                                                                                                                                    >  1:923980/1‑150 (MQ=255)
    tCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCtt                                                                                                                                                   <  2:403004/150‑1 (MQ=255)
    tCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCtt                                                                                                                                                   <  2:2854514/150‑1 (MQ=255)
      tGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTcc                                                                                                                                                 <  1:1069846/150‑1 (MQ=255)
       gACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCa                                                                                                                                                <  1:4041344/150‑1 (MQ=255)
          gCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAgcg                                                                                                                                             <  1:2845741/150‑1 (MQ=255)
            tGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGcc                                                                                                                                           >  1:1686646/1‑150 (MQ=255)
             gggTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCt                                                                                                                                          <  2:3171541/150‑1 (MQ=255)
              ggTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATAGCTTCCAGCGCCTg                                                                                                                                         <  2:3097670/150‑1 (MQ=255)
               gTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTgg                                                                                                                                        <  1:357428/150‑1 (MQ=255)
               gTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTgg                                                                                                                                        <  1:357351/150‑1 (MQ=255)
                tCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGc                                                                                                                                       <  2:430644/150‑1 (MQ=255)
                 cGCCATCGGCGTATGATACGTCAACGCAGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGAAGGCCCTAAGCAATCGCGACGGCTGCAAATCAGATGCTATTTGCCCATAGCTTCCAGCGCCTGGCg                                                                                                                                      >  2:3873817/1‑150 (MQ=255)
                    cATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGAt                                                                                                                                   <  2:3429102/150‑1 (MQ=255)
                     aTCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGAtt                                                                                                                                  >  1:1492852/1‑150 (MQ=255)
                       cGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTaa                                                                                                                                >  2:1727035/1‑150 (MQ=255)
                        ggCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAAc                                                                                                                               <  2:2996993/150‑1 (MQ=255)
                          cGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACgg                                                                                                                             <  2:198410/150‑1 (MQ=255)
                           gTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGc                                                                                                                            >  2:802376/1‑150 (MQ=255)
                           gTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGc                                                                                                                            >  1:135417/1‑150 (MQ=255)
                            tATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGcc                                                                                                                           >  1:3963020/1‑150 (MQ=255)
                            tATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGcc                                                                                                                           >  1:1078574/1‑150 (MQ=255)
                            tATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGcc                                                                                                                           >  2:2540310/1‑150 (MQ=255)
                            tATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGcc                                                                                                                           >  1:333705/1‑150 (MQ=255)
                                  tCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTg                                                                                                                     <  2:537998/150‑1 (MQ=255)
                                    gTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGcc                                                                                                                   >  2:2909083/1‑150 (MQ=255)
                                    gTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGcc                                                                                                                   <  1:1396981/150‑1 (MQ=255)
                                     tCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCg                                                                                                                  <  2:2555143/150‑1 (MQ=255)
                                         cGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATc                                                                                                              >  2:3581470/1‑150 (MQ=255)
                                         cGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATc                                                                                                              >  1:2639011/1‑150 (MQ=255)
                                         cGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATc                                                                                                              >  1:2639465/1‑150 (MQ=255)
                                         cGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATc                                                                                                              <  1:2709292/150‑1 (MQ=255)
                                            cGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTg                                                                                                           >  2:2154424/1‑150 (MQ=255)
                                            cGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTg                                                                                                           <  2:3245919/150‑1 (MQ=255)
                                             ggACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTgg                                                                                                          <  1:1619524/150‑1 (MQ=255)
                                             ggACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTgg                                                                                                          <  2:1434256/150‑1 (MQ=255)
                                              gACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGt                                                                                                         >  1:2518903/1‑150 (MQ=255)
                                                    tATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGca                                                                                                   <  1:710612/150‑1 (MQ=255)
                                                     aTTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGcaa                                                                                                  >  2:1472217/1‑150 (MQ=255)
                                                     aTTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGcaa                                                                                                  >  2:1404787/1‑150 (MQ=255)
                                                         cGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCaacaac                                                                                              <  2:924038/150‑1 (MQ=255)
                                                          ggTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACg                                                                                             <  1:1112310/150‑1 (MQ=255)
                                                            ttCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACTGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCt                                                                                           <  2:1589875/150‑1 (MQ=255)
                                                             tCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCtt                                                                                          <  1:3774668/150‑1 (MQ=255)
                                                                ggCGACGAGCTTCACGCGTCTCTTCAAGGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTACCGCAACAACGCTTTAt                                                                                       >  2:920251/1‑150 (MQ=255)
                                                                ggCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTAt                                                                                       <  1:3549749/150‑1 (MQ=255)
                                                                  cgacgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGc                                                                                     <  1:2358067/150‑1 (MQ=255)
                                                                  cgacgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGc                                                                                     <  1:3218086/150‑1 (MQ=255)
                                                                   gacgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCa                                                                                    >  2:3631423/1‑150 (MQ=255)
                                                                     cgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCaaa                                                                                  >  2:2777440/1‑150 (MQ=255)
                                                                     cgaGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCaaa                                                                                  >  2:689709/1‑150 (MQ=255)
                                                                         cTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAacgac                                                                              <  2:1844089/150‑1 (MQ=255)
                                                                              cGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGt                                                                         >  2:113785/1‑150 (MQ=255)
                                                                              cGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGt                                                                         >  2:2058707/1‑150 (MQ=255)
                                                                                cATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGttt                                                                       <  2:1340195/150‑1 (MQ=255)
                                                                                cATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGttt                                                                       >  1:3844193/1‑150 (MQ=255)
                                                                                cATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGttt                                                                       <  1:1948517/150‑1 (MQ=255)
                                                                                cATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGttt                                                                       >  1:297420/1‑150 (MQ=255)
                                                                                cATCTCTTCAAGGGTTGGCGCAGGCACTACGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCTGTTTGCCGGATCGTGGTACCGCAACAACGCTTTTTGCAAACGACGCTGttt                                                                       <  1:1049009/150‑1 (MQ=255)
                                                                                  tctcTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGt                                                                     <  2:431746/150‑1 (MQ=255)
                                                                                  tctcTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGt                                                                     <  2:131653/150‑1 (MQ=255)
                                                                                    tctTCATTGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCTGCGCCTGGCGGTTTAACGGCCAGTTTGCCGGATCGTGTTAACGCAAAAACGCTTTATGCAAACGACGCTGTTTGTCg                                                                   <  1:1323314/150‑1 (MQ=255)
                                                                                    tctTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCg                                                                   <  2:3478840/150‑1 (MQ=255)
                                                                                        cAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCt                                                               >  1:1865054/1‑150 (MQ=255)
                                                                                          aTGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTc                                                             <  1:2338722/150‑1 (MQ=255)
                                                                                           tGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCg                                                            <  1:2525491/150‑1 (MQ=255)
                                                                                            ggTGGGCGCAGGCAATAAGCAATCGCGCCGGGTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCTGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCgg                                                           <  1:3385160/150‑1 (MQ=255)
                                                                                              ttGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTa                                                         <  2:2378824/150‑1 (MQ=255)
                                                                                               tGGCGCAGGCACTACGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGTATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTAc                                                        <  1:2482010/150‑1 (MQ=255)
                                                                                               tGGCGCAGGCACTAAGAAATCGCGACGGCTGCCAAACAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGTATTAACGGCCCGTTTGCCGTATCGTGGTAACGCAACACCGCTTTATGCACACGACGCTGTTTGTCGCCCTTCGGTAc                                                        <  1:3912634/150‑1 (MQ=255)
                                                                                                ggCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACg                                                       <  1:4064537/150‑1 (MQ=255)
                                                                                                  cgcAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGaa                                                     >  2:2202530/1‑150 (MQ=255)
                                                                                                  cgcAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGaa                                                     <  2:1209540/150‑1 (MQ=255)
                                                                                                          cTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTGGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCtt                                             <  2:1951009/150‑1 (MQ=255)
                                                                                                          cTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATAGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCtt                                             <  1:3265212/150‑1 (MQ=255)
                                                                                                              gCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATAGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGct                                         >  2:942169/1‑150 (MQ=255)
                                                                                                                aaTCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGctct                                       >  2:3727041/1‑150 (MQ=255)
                                                                                                                   cgcgACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTtat                                    <  1:423180/150‑1 (MQ=255)
                                                                                                                     cgACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATaa                                  >  1:1270336/1‑150 (MQ=255)
                                                                                                                     cgACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATaa                                  <  2:1354549/150‑1 (MQ=255)
                                                                                                                     cgACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATaa                                  <  1:2466302/150‑1 (MQ=255)
                                                                                                                     cgACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATaa                                  <  2:4115638/150‑1 (MQ=255)
                                                                                                                      gACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAAc                                 >  1:793295/1‑150 (MQ=255)
                                                                                                                       aCGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAAcc                                >  1:1159782/1‑150 (MQ=255)
                                                                                                                       aCGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAAcc                                >  2:2277433/1‑150 (MQ=255)
                                                                                                                             gCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCtt                          <  1:2173159/150‑1 (MQ=255)
                                                                                                                                 aTCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGcc                      <  1:1794189/150‑1 (MQ=255)
                                                                                                                                  tCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCa                     <  2:2445628/150‑1 (MQ=255)
                                                                                                                                  tCAGATGCTTTTTGCCCATAGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCa                     >  2:1566009/1‑150 (MQ=255)
                                                                                                                                   cAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCGGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAg                    <  1:748441/150‑1 (MQ=255)
                                                                                                                                           ttttGCCCATCGCTTCCAGCGCCTGGCGGATTACCGGCCAGTGTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGtttt            <  1:728606/150‑1 (MQ=255)
                                                                                                                                           ttttGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGtttt            >  2:1644714/1‑150 (MQ=255)
                                                                                                                                            tttGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGttttt           <  2:3233101/150‑1 (MQ=255)
                                                                                                                                                 ccATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAATACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCgg      >  2:1841494/1‑150 (MQ=255)
                                                                                                                                                 ccATAGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCgg      >  1:1565759/1‑150 (MQ=255)
                                                                                                                                                   aTCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGgtg    <  1:3992522/150‑1 (MQ=255)
                                                                                                                                                     cGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCCGCGGGTTTTTCCCGGTGTa  >  2:4056090/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGTCTGACGCTGGGTCGCCATCGGCGTATGTTTCGTCAACGCCGGACGGGTATTGCGGTTCTGGCGACGAGCTTCACGCATCTCTTCAATGGTTGGCGCAGGCACTAAGCAATCGCGACGGCTGCCAATCAGATGCTTTTTGCCCATCGCTTCCAGCGCCTGGCGGATTAACGGCCAGTTTGCCGGATCGTGGTAACGCAACAACGCTTTATGCAAACGACGCTGTTTGTCGCCCTTCGGTACGAAGACGTCTTCGCTCTTATAACCAATCTTCGCCAGCGGGTTTTTCCCGGTGTA  >  CP009273/3152374‑3152672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: