Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,161,800 A→C 7.0% intergenic (+116/‑185) rrfB → / → murB 5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,161,8000AC 7.0% 413.0 / 4.2 168intergenic (+116/‑185)rrfB/murB5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (83/69);  new base C (5/11);  total (88/80)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

TCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAC  >  CP009273/4161653‑4161947
                                                                                                                                                   |                                                                                                                                                   
tcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg                                                                                                                                                    >  1:3270742/1‑149 (MQ=21)
  cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTCTCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGTAGCa                                                                                                                                                  >  1:4334249/1‑149 (MQ=18)
  cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                                                                                                                  >  2:191388/1‑149 (MQ=33)
  cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCa                                                                                                                                                  >  1:1717962/1‑149 (MQ=18)
    cATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc                                                                                                                                                <  2:4647330/149‑1 (MQ=32)
    cATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGGTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCAAc                                                                                                                                                >  2:3455230/1‑149 (MQ=17)
      tGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                                                                                                                              >  2:356790/1‑149 (MQ=32)
      tGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                                                                                                                              <  1:2089908/149‑1 (MQ=32)
       gCGAGAGTAGGTAACTGCCAGGCATAAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGc                                                                                                                                             <  1:2895040/149‑1 (MQ=17)
       gCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGc                                                                                                                                             >  1:4585801/1‑149 (MQ=32)
         gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                           >  2:531713/1‑149 (MQ=32)
           gagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                                                                                                                         >  1:85761/1‑149 (MQ=32)
           gagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAGTTGAACGTTGCGAAGCAACGGCCCgg                                                                                                                                         >  2:2183054/1‑149 (MQ=17)
                gggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                                    <  2:2817663/149‑1 (MQ=32)
                gggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTGTTCTCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGCGCGGATTTGAACGTTGAGAAGCACCGGCCCGGAGGGt                                                                                                                                    >  2:2935886/1‑149 (MQ=11)
                  gAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                                                                  >  1:3361565/1‑149 (MQ=32)
                   aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc                                                                                                                                 >  2:237533/1‑149 (MQ=32)
                       gCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAATAATGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGc                                                                                                                             >  1:2024329/1‑149 (MQ=11)
                       gCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGc                                                                                                                             >  1:4639573/1‑149 (MQ=11)
                       gCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTGTCGTTTTATCTGTTGTTTGTCGGTGACCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGgggcgggc                                                                                                                             >  1:3117772/1‑149 (MQ=2)
                          aggccacAAATAAACCGACAGGCTCCGTCGAACGACTGGGCCTTCCGTTTTTTCTGTGGTGTGTCGGGGAACGCTCCCCTGAGTAGGACAACGCCGCCGGGAGCGGATTTGAAGGTTGCTACGCACCGGCCCGGAGGGTGGCGGGCAgg                                                                                                                          <  2:1834927/143‑1 (MQ=1)
                           gccaacAAATAAAACGAAAGGCTCAGTCGAAAGACGGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGa                                                                                                                         <  2:3761176/144‑1 (MQ=0)
                                                 cctgTCGAAAGACGGGGCCTTTCTTTTTATCTGTTTTTTGTGGGTGAACGCTCTCCTGAGTAGTACAAATCCGCGGGGAGCGTATTTGAACGTTGCTAACCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                                                   <  1:3363616/146‑1 (MQ=11)
                                                  ctggcGAATGAGTGGGCATTTGTTTTTTTGTGTTTTTTGTGGGTGCACGCTCTCCTGAGTTGCAAAACGCCGCGGGGAGGGGTTTTGACGGTGGCGACCCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                  <  1:623665/145‑1 (MQ=2)
                                                  cAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                  >  2:4045836/1‑149 (MQ=31)
                                                  cAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                  >  1:4276378/1‑149 (MQ=31)
                                                  cAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGGGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                  <  2:373039/149‑1 (MQ=14)
                                                  cAGTCGAAAGACTGGGCCTTACGTTTTATCTGTGGTTTGTCGGTGAACGCGCTCCGGAGTAGGACAACTCCGCGGGGAGCGGATTTGAACGTTGCGAAGCCACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                  <  2:3212033/149‑1 (MQ=11)
                                                    ttctaaaGAGTGGGCCTTTGTTTTTATGTGTTGTTGGTGGGTTAACGCTCTCGTGATTTGGAAAAAGCCGCGGGGTGGTTTTTTTAACGTTGCTACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                <  1:4270369/145‑1 (MQ=2)
                                                      cGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaa                                                                                              >  2:4807428/1‑149 (MQ=32)
                                                       gAAAGACTGGGCCTTTCGTTTTATCTGTTTTTTGTGGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                             <  2:4192822/149‑1 (MQ=14)
                                                       gAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                             >  2:4305196/1‑149 (MQ=31)
                                                       gAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                             <  2:4193479/149‑1 (MQ=31)
                                                       gAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                             >  1:303248/1‑149 (MQ=31)
                                                       gAAAGACTGGGCATTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                             <  2:2940098/149‑1 (MQ=14)
                                                        aaaGACTGGCCCTTTCGTTTTATATGTTGTTTGTCGGTGAACGCTCTCCTTAGTAGGACAAATCCGCCGGGTGCGTATTTTAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGACCGCCATAAACTGCCAGGCATCAAAt                                                                                            <  1:863719/149‑1 (MQ=11)
                                                                                              gAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGt                                                      >  2:3580174/1‑149 (MQ=31)
                                                                                               aaCGCTCTCCTGAGTCGGACAAATCCGCCGGGAGCGTTTTTGAACTTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGtt                                                     <  1:2074241/149‑1 (MQ=255)
                                                                                               aaCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCggt                                                     >  2:4477006/1‑147 (MQ=11)
                                                                                               aaCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGtt                                                     <  2:654929/149‑1 (MQ=32)
                                                                                                ccGCTCTCCTGAGTAGGACAAACCCGCCGGGCGCGGTTTTGAACGTGGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                    <  1:2095553/148‑1 (MQ=255)
                                                                                                aCGCTCTCCTGAGTCGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                    <  1:2095729/149‑1 (MQ=17)
                                                                                                aCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                    >  2:1720320/1‑149 (MQ=32)
                                                                                                 cGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTc                                                   <  2:577366/149‑1 (MQ=33)
                                                                                                  gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCAAAAACAGCCAGGCATCAAATTAAGCAGAATGCCAACCAGACGGATGGCCTTTTTTCGTTTCt                                                  >  1:3194815/1‑149 (MQ=11)
                                                                                                  gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGTCGGATGGCCTTTTTGCGTTTCt                                                  >  1:5112482/1‑149 (MQ=21)
                                                                                                  gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAAAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCAGACGGATTGCCTTTTTGCGTTTCt                                                  >  1:2862695/1‑149 (MQ=255)
                                                                                                        cTGCGTAGGCCAAATCCGCCGGTCGCGTTTTTGAAGTTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc                                            <  1:2111340/149‑1 (MQ=14)
                                                                                                         tgtgTAGGACAAATCCGGGGGGTGCGTTTTTGAAGGTTGGGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                                           <  1:3158632/146‑1 (MQ=14)
                                                                                                           aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAActct                                         <  2:2562414/149‑1 (MQ=255)
                                                                                                            gTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCtt                                        >  1:4805069/1‑149 (MQ=255)
                                                                                                            gTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCAGACGGATGGCCTTTTTGCGTTTCTACAAACTCtt                                        >  1:5230871/1‑149 (MQ=255)
                                                                                                             tcggACAAATCCGCCGGGTGCGTATTTTAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTc                                       <  1:5178126/147‑1 (MQ=255)
                                                                                                             tAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCATAAACTGCCAGGCCTCAAATTAAGCAGAAGGCCATCCAGACGGAAGGCCGTTTTGCGTTTCTACAAACTCtta                                       >  1:4538500/1‑148 (MQ=14)
                                                                                                             tAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCTCCTGACGGATGGGCTTTTTGCGGTTGTACAAACTCTTc                                       >  1:220691/1‑149 (MQ=255)
                                                                                                             tAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTc                                       >  2:5106720/1‑149 (MQ=255)
                                                                                                              atgACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGGTTCTACAAACTCTTcc                                      >  2:275201/3‑149 (MQ=255)
                                                                                                              aGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTcc                                      >  2:275270/1‑149 (MQ=255)
                                                                                                              aGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTcc                                      >  2:3253961/1‑149 (MQ=255)
                                                                                                               ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCTCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt                                     >  1:2179513/1‑149 (MQ=255)
                                                                                                               ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt                                     <  2:3361565/149‑1 (MQ=255)
                                                                                                               ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATACTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt                                     >  1:3500835/1‑149 (MQ=255)
                                                                                                               ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAGAAACAGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATTGCCATTTTTCGTTTCTAAAAACTCTTACt                                     >  1:4625722/1‑149 (MQ=14)
                                                                                                                gACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCAGCAAATTAAGCAGAAGGCCCTCCTGACGGATGGCCTTTTTGCGGGTCTACAAACTCTTCCTg                                    >  2:1617946/1‑149 (MQ=255)
                                                                                                                 aCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgt                                   <  1:1312408/149‑1 (MQ=255)
                                                                                                                    aaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgt                                <  2:661286/149‑1 (MQ=255)
                                                                                                                      tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCa                              <  2:4687436/149‑1 (MQ=255)
                                                                                                                      tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCa                              <  1:1174754/149‑1 (MQ=255)
                                                                                                                       ccgccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCat                             <  2:3028143/149‑1 (MQ=255)
                                                                                                                         gccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCatat                           <  2:4082370/149‑1 (MQ=255)
                                                                                                                           ggggAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCt                         <  2:3944866/148‑1 (MQ=255)
                                                                                                                           cgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAGTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCt                         >  2:2867660/1‑149 (MQ=255)
                                                                                                                            gggAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTGTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTa                        >  2:359010/1‑149 (MQ=255)
                                                                                                                             ggcgCGTATTTGACCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAc                       <  1:4175899/146‑1 (MQ=255)
                                                                                                                              gAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCAGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCAGTCGTCAAATCTACa                      >  1:1272527/1‑149 (MQ=255)
                                                                                                                                 cGGATTTTAACGTTGCTAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGc                   <  1:1557761/149‑1 (MQ=255)
                                                                                                                                  ggATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGcc                  <  2:3125958/149‑1 (MQ=255)
                                                                                                                                   gATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCAGACGGATGGCCTTTTTGCGTTTCTACAAACACTTCCGGTCGTCATATCTACAAGCCa                 >  2:4844807/1‑149 (MQ=255)
                                                                                                                                   gATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCCCCATAAACAGCCAGGCATCAAATTAAGCAGAAGGCCATCCGGACGGGTGGCCTTTTTGCGGTTCTTCAAAAACTTGCTGTCGTAATAAATACAAAccc                 >  2:2252207/1‑148 (MQ=11)
                                                                                                                                    aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAt                <  2:5152513/149‑1 (MQ=255)
                                                                                                                                     tttGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCGTTTTGCGTTTCTACAAACTCTTCCTGTCGTCCTATCTACAAGCCATc               >  2:3281166/1‑149 (MQ=255)
                                                                                                                                       tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTGTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc             >  2:1980426/1‑149 (MQ=255)
                                                                                                                                       tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc             <  2:3131134/149‑1 (MQ=255)
                                                                                                                                        gAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcccc            <  2:3799359/149‑1 (MQ=255)
                                                                                                                                         aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccccc           >  2:312758/1‑149 (MQ=255)
                                                                                                                                          aCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCCGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcccccc          >  2:868449/1‑149 (MQ=255)
                                                                                                                                          aCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCTCCTGACGGCTGGCCGTTTTGCGTTTCTACCAACACTTCCTGTCGTCATATCTACAAGCCATcccccc          >  2:4541637/1‑149 (MQ=18)
                                                                                                                                          aCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCCTcccccc          >  2:370249/1‑149 (MQ=255)
                                                                                                                                            gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCTAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac        <  2:4564724/149‑1 (MQ=255)
                                                                                                                                            gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCTAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac        <  2:4490391/149‑1 (MQ=255)
                                                                                                                                            gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac        >  2:4440583/1‑149 (MQ=255)
                                                                                                                                            gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAAGCCCGCCATAAACTGCTAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac        <  2:4490087/149‑1 (MQ=255)
                                                                                                                                               gCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGa     <  2:2169841/149‑1 (MQ=255)
                                                                                                                                                 gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATa   <  2:1630358/149‑1 (MQ=255)
                                                                                                                                                 gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATa   <  2:142673/149‑1 (MQ=255)
                                                                                                                                                 gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATa   <  1:3937535/149‑1 (MQ=255)
                                                                                                                                                 gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATa   <  1:3373338/149‑1 (MQ=255)
                                                                                                                                                  aaGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAc  >  2:851202/1‑149 (MQ=255)
                                                                                                                                                  aaGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATACTGTCGGATGGCCTTTTTTGGGTTCTACAAAATATTCCTGGAGTCATATCTAAAATCCATCCCCCCACCGATAc  >  2:492090/1‑149 (MQ=17)
                                                                                                                                                  aaGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAAAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAAACAGAAGGATGGGCTTTTTTGGTTTCTACAAAAAATTACTGTAGTCATATCTAAAATCCAAACCCCCACACATaa  >  1:1238648/1‑148 (MQ=14)
                                                                                                                                                   |                                                                                                                                                   
TCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAC  >  CP009273/4161653‑4161947

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