Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,034 G→T 13.8% intergenic (+15/‑71) miaA → / → hfq delta(2)‑isopentenylpyrophosphate tRNA‑adenosine transferase/global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,390,0340GT13.8% 810.9 / 73.4 254intergenic (+15/‑71)miaA/hfqdelta(2)‑isopentenylpyrophosphate tRNA‑adenosine transferase/global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base G (86/132);  new base T (18/17);  total (104/150)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

TGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGG  >  CP009273/4389885‑4390183
                                                                                                                                                     |                                                                                                                                                     
      aCGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTAt                                                                                                                                                 <  1:255218/150‑1 (MQ=255)
      aCGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTAt                                                                                                                                                 <  2:1891816/150‑1 (MQ=255)
       cGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATc                                                                                                                                                >  1:1063473/1‑150 (MQ=255)
       cGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATc                                                                                                                                                <  2:2030408/150‑1 (MQ=255)
       cGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATc                                                                                                                                                <  1:312561/150‑1 (MQ=255)
          gaCAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTg                                                                                                                                             >  1:3179804/1‑150 (MQ=255)
           aCAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTgc                                                                                                                                            <  2:1766855/150‑1 (MQ=255)
            cAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTgcg                                                                                                                                           >  2:3953948/1‑150 (MQ=255)
            cAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTGGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTAGCGTgcg                                                                                                                                           >  2:2929724/1‑150 (MQ=255)
             aGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTgcgc                                                                                                                                          >  2:3825938/1‑150 (MQ=255)
             aGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTgcgc                                                                                                                                          >  1:729029/1‑150 (MQ=255)
                    gAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAAttttt                                                                                                                                   >  2:3691784/1‑150 (MQ=255)
                      aGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTc                                                                                                                                 <  2:3638810/150‑1 (MQ=255)
                        cGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAg                                                                                                                               >  1:3782269/1‑150 (MQ=255)
                           cAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAAt                                                                                                                            <  1:1523194/150‑1 (MQ=255)
                            aGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATc                                                                                                                           <  2:2145910/150‑1 (MQ=255)
                             gATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCg                                                                                                                          <  1:4118041/150‑1 (MQ=255)
                               tAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGaa                                                                                                                        <  2:2573356/150‑1 (MQ=255)
                               tAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGaa                                                                                                                        <  1:2525809/150‑1 (MQ=255)
                               tAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGaa                                                                                                                        <  1:2178745/150‑1 (MQ=255)
                                aaCCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGaaa                                                                                                                       >  2:1865156/1‑150 (MQ=255)
                                  ccTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAgg                                                                                                                     <  2:2059875/150‑1 (MQ=255)
                                  ccTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAgg                                                                                                                     <  1:204100/150‑1 (MQ=255)
                                     ggCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTc                                                                                                                  >  2:1647877/1‑150 (MQ=255)
                                       cTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCaa                                                                                                                >  1:4028555/1‑150 (MQ=255)
                                         gCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAg                                                                                                              <  1:1027248/150‑1 (MQ=255)
                                                    aaGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGc                                                                                                   >  1:656695/1‑150 (MQ=255)
                                                    aaGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGc                                                                                                   <  2:506640/150‑1 (MQ=255)
                                                      ggggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCat                                                                                                 <  1:3408573/150‑1 (MQ=255)
                                                       gggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGTCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCata                                                                                                >  1:4312352/1‑150 (MQ=255)
                                                        ggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACTAATAAGCatat                                                                                               >  2:2439691/1‑150 (MQ=255)
                                                        ggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCatat                                                                                               >  2:1172792/1‑150 (MQ=255)
                                                        ggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCatat                                                                                               <  2:2244970/150‑1 (MQ=255)
                                                         gTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCAtata                                                                                              >  2:3070855/1‑150 (MQ=255)
                                                          ttCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCTTGCGCAATTTTTTCAGAATCGAAAGTTTCAAAGTACAAATAAGCATATaa                                                                                             >  2:4281193/1‑150 (MQ=255)
                                                           tCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGCCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAg                                                                                            <  2:4211910/150‑1 (MQ=255)
                                                           tCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAg                                                                                            <  2:4212458/150‑1 (MQ=255)
                                                             aCTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGa                                                                                          >  2:1259805/1‑150 (MQ=255)
                                                             aCTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATAGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCTAAAGGTTCAAAGTACAAATAAGCATATAAGGa                                                                                          >  2:3416803/1‑150 (MQ=255)
                                                              cTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGaa                                                                                         <  2:2464891/150‑1 (MQ=255)
                                                                ggCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGaaaa                                                                                       <  2:656695/150‑1 (MQ=255)
                                                                  cTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAaga                                                                                     <  1:2693678/150‑1 (MQ=255)
                                                                    tGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAagaga                                                                                   <  2:4543517/150‑1 (MQ=255)
                                                                    tGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAagaga                                                                                   <  2:2850849/150‑1 (MQ=255)
                                                                     gACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAagagag                                                                                  >  1:3649524/1‑150 (MQ=255)
                                                                     gACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAagagag                                                                                  <  1:3471055/150‑1 (MQ=255)
                                                                      aCAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAgagaga                                                                                 <  1:3200355/150‑1 (MQ=255)
                                                                      aCAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAgagaga                                                                                 <  1:1617092/150‑1 (MQ=255)
                                                                        aGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAAt                                                                               <  1:1534832/150‑1 (MQ=255)
                                                                          tGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATgg                                                                             <  1:3746127/150‑1 (MQ=255)
                                                                            aaaaaCCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCt                                                                           <  1:312650/150‑1 (MQ=255)
                                                                            aaaaaCCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCt                                                                           >  2:1362489/1‑150 (MQ=255)
                                                                            aaaaaCCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCt                                                                           >  2:1362658/1‑150 (MQ=255)
                                                                             aaaaCCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTa                                                                          <  2:2900227/150‑1 (MQ=255)
                                                                             aaaaCCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTa                                                                          >  2:2530448/1‑150 (MQ=255)
                                                                                 ccAGAACAGGCGCGTGACGAATTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTTAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAggg                                                                      <  2:4398818/150‑1 (MQ=255)
                                                                                        aGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCt                                                               <  1:3669279/150‑1 (MQ=255)
                                                                                         ggCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGATAATGGCTAAGGGGCAATCtt                                                              <  1:4279866/150‑1 (MQ=255)
                                                                                         ggCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCtt                                                              <  2:4438267/150‑1 (MQ=255)
                                                                                           cgcgTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTa                                                            <  1:1580905/150‑1 (MQ=255)
                                                                                           cgcgTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTa                                                            <  1:3845907/150‑1 (MQ=255)
                                                                                             cgTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACa                                                          <  1:3065749/150‑1 (MQ=255)
                                                                                              gTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACaa                                                         >  1:3304138/1‑150 (MQ=255)
                                                                                              gTGACGAAGTATTACAGGGTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACaa                                                         >  1:1479344/1‑150 (MQ=255)
                                                                                                gACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGa                                                       <  2:4054225/150‑1 (MQ=255)
                                                                                                  cGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTAACGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATc                                                     >  2:1260406/1‑150 (MQ=255)
                                                                                                   gAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTTAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATcc                                                    <  2:4187040/150‑1 (MQ=255)
                                                                                                   gAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATcc                                                    <  1:834761/150‑1 (MQ=255)
                                                                                                   gAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAAGGGCTAAGGGGCAATCTTTACAAGATcc                                                    <  1:1295337/150‑1 (MQ=255)
                                                                                                   gAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATcc                                                    >  2:2181238/1‑150 (MQ=255)
                                                                                                    aaGTATTACAGGTTGTTGGTGCTCTCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGATAGGTTCAAAGTACAATTAAGCATATAAGGAAAAGAGAGAATGGCTCAGGTGCAATCTTTACAAGATCCg                                                   >  2:604904/1‑150 (MQ=255)
                                                                                                     aGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGt                                                  >  2:580682/1‑150 (MQ=255)
                                                                                                      gTATTACATGTTTTTGGTGATATCGCAGGCTGAATGTGTACAATTGAGAAGTATCGTGCGCTATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAATCATATAATGAAAAGAGAGAATGGCTAAGGGGCAAACTTTACAAGATCCGtt                                                 <  2:225926/150‑1 (MQ=255)
                                                                                                      gTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGtt                                                 >  2:2422005/1‑150 (MQ=255)
                                                                                                         ttACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCt                                              <  1:1980629/150‑1 (MQ=255)
                                                                                                             aGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAAc                                          >  2:4209091/1‑150 (MQ=255)
                                                                                                              ggTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                         >  2:1591994/1‑150 (MQ=255)
                                                                                                              ggTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                         >  1:896602/1‑150 (MQ=255)
                                                                                                              ggTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                         >  2:1542165/1‑150 (MQ=255)
                                                                                                              ggTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                                         >  2:1543867/1‑150 (MQ=255)
                                                                                                                  gttgGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACt                                     <  1:593725/150‑1 (MQ=255)
                                                                                                                        gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCg                               <  2:270775/150‑1 (MQ=255)
                                                                                                                        gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCg                               <  1:851235/150‑1 (MQ=255)
                                                                                                                        gCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCg                               <  2:3079574/150‑1 (MQ=255)
                                                                                                                            tCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGaa                           <  1:1778039/150‑1 (MQ=255)
                                                                                                                             cGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAAc                          <  1:4532272/150‑1 (MQ=255)
                                                                                                                             cGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAAc                          <  1:1175702/150‑1 (MQ=255)
                                                                                                                              gCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACg                         <  1:3363453/150‑1 (MQ=255)
                                                                                                                                aGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACgtg                       >  1:4488072/1‑150 (MQ=255)
                                                                                                                                   cTGAATGTGTACAATTGATACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTc                    <  1:3953948/150‑1 (MQ=255)
                                                                                                                                   cTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTc                    <  1:3622233/150‑1 (MQ=255)
                                                                                                                                    tGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                   >  2:1578289/1‑150 (MQ=255)
                                                                                                                                        tgtgTACAATTGAGACGTCTCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGtt               >  1:944639/1‑150 (MQ=255)
                                                                                                                                           gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAATGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCt            >  1:1797577/1‑150 (MQ=255)
                                                                                                                                            tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta           >  2:292435/1‑150 (MQ=255)
                                                                                                                                            tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta           >  1:2835120/1‑150 (MQ=255)
                                                                                                                                            tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta           >  2:2186799/1‑150 (MQ=255)
                                                                                                                                              cAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtatt         <  1:1544326/150‑1 (MQ=255)
                                                                                                                                                     gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGATAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTgg  <  1:2931378/150‑1 (MQ=255)
                                                                                                                                                     gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTgg  <  2:2512282/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGG  >  CP009273/4389885‑4390183

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: