Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,672,064 T→A 34.9% T162S (ACC→TCC)  yhjH ← cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,672,0640TA34.9% 329.6 / inf 284T162S (ACC→TCC) yhjHcyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated
Reads supporting (aligned to +/- strand):  ref base T (92/93);  new base A (54/45);  total (146/138)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 285 total aligned reads displayed.

CCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAGGG  >  CP009273/3671915‑3672210
                                                                                                                                                     |                                                                                                                                                  
cccGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGt                                                                                                                                                    <  2:4148522/150‑1 (MQ=255)
   gcgACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTAcc                                                                                                                                                 >  2:1149083/1‑150 (MQ=255)
    cgACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCa                                                                                                                                                >  1:2864384/1‑150 (MQ=255)
     gACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCaa                                                                                                                                               >  1:4272524/1‑150 (MQ=255)
        aaTAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAat                                                                                                                                            >  1:3717743/1‑150 (MQ=255)
              gATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCa                                                                                                                                      >  1:2045039/1‑150 (MQ=255)
              gATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCa                                                                                                                                      <  1:1371620/150‑1 (MQ=255)
                  cATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCca                                                                                                                                  <  1:976979/150‑1 (MQ=255)
                  cATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCca                                                                                                                                  <  1:995778/150‑1 (MQ=255)
                  cATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCca                                                                                                                                  >  2:512445/1‑150 (MQ=255)
                    ttAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCcaca                                                                                                                                <  1:2127485/150‑1 (MQ=255)
                      aGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGc                                                                                                                              <  1:3517267/150‑1 (MQ=255)
                          gTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGc                                                                                                                          >  1:3608809/1‑150 (MQ=255)
                           tAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGcc                                                                                                                         <  2:1544988/150‑1 (MQ=255)
                             aaaGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCaa                                                                                                                       <  2:2263490/150‑1 (MQ=255)
                               aGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAAt                                                                                                                     <  1:638612/150‑1 (MQ=255)
                               aGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAAt                                                                                                                     >  1:2751411/1‑150 (MQ=255)
                               aGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAAt                                                                                                                     >  2:1804184/1‑150 (MQ=255)
                                gCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAAtt                                                                                                                    <  1:3029640/150‑1 (MQ=255)
                                    agaAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTcaca                                                                                                                >  1:1101268/1‑150 (MQ=255)
                                     gaAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTCCCAAAATCATCCAGCCACAGCGGGCCAAATTcacac                                                                                                               >  2:2149182/1‑150 (MQ=255)
                                      aaaaGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTcacaca                                                                                                              >  1:2530546/1‑150 (MQ=255)
                                       aaaGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACAt                                                                                                             >  1:1732229/1‑150 (MQ=255)
                                           agTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAg                                                                                                         >  1:3604532/1‑150 (MQ=255)
                                           agTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAg                                                                                                         >  1:2727145/1‑150 (MQ=255)
                                            gtgtGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAgg                                                                                                        >  1:3693320/1‑150 (MQ=255)
                                            gtgtGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAgg                                                                                                        >  1:3211277/1‑150 (MQ=255)
                                            gtgtGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAgg                                                                                                        <  2:3906082/150‑1 (MQ=255)
                                             tgtgCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGc                                                                                                       <  2:1204081/150‑1 (MQ=255)
                                                 cGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAg                                                                                                   >  1:1213549/1‑150 (MQ=255)
                                                  gTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAgg                                                                                                  >  2:4480693/1‑150 (MQ=255)
                                                   tCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGt                                                                                                 >  2:1987425/1‑150 (MQ=255)
                                                       tCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGaa                                                                                             >  2:3343998/1‑150 (MQ=255)
                                                       tCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGaa                                                                                             <  1:79752/150‑1 (MQ=255)
                                                       tCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGaa                                                                                             >  2:695718/1‑150 (MQ=255)
                                                       tCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGaa                                                                                             <  2:2475928/150‑1 (MQ=255)
                                                          ggCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCt                                                                                          >  2:3425960/1‑150 (MQ=255)
                                                           gCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCtt                                                                                         <  1:891348/150‑1 (MQ=255)
                                                           gCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCtt                                                                                         <  2:1903687/150‑1 (MQ=255)
                                                             gACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTc                                                                                       <  2:2975168/150‑1 (MQ=255)
                                                              aCTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCg                                                                                      <  1:2064735/150‑1 (MQ=255)
                                                               cTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCgg                                                                                     >  1:1276588/1‑150 (MQ=255)
                                                                tGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGc                                                                                    >  2:675395/1‑150 (MQ=255)
                                                                tGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGc                                                                                    <  2:365445/150‑1 (MQ=255)
                                                                tGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGc                                                                                    >  1:3249870/1‑150 (MQ=255)
                                                                  aCGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAg                                                                                  <  1:806448/150‑1 (MQ=255)
                                                                  aCGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAg                                                                                  <  2:2345430/150‑1 (MQ=255)
                                                                    gcagcaTCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGAc                                                                                <  1:3776272/150‑1 (MQ=255)
                                                                    gcagcaTCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGAc                                                                                <  1:3069420/150‑1 (MQ=255)
                                                                      agcaTCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACgg                                                                              >  1:1687341/1‑150 (MQ=255)
                                                                               aaCAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTcc                                                                     <  1:1511416/150‑1 (MQ=255)
                                                                                 cAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTccac                                                                   <  1:1886866/150‑1 (MQ=255)
                                                                                 cAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTccac                                                                   <  2:1441778/150‑1 (MQ=255)
                                                                                      ctcGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGtt                                                              >  1:3643690/1‑150 (MQ=255)
                                                                                       tcGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTc                                                             <  1:3659779/150‑1 (MQ=255)
                                                                                        cgcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCg                                                            >  2:2817743/1‑150 (MQ=255)
                                                                                        cgcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCg                                                            >  2:3315747/1‑150 (MQ=255)
                                                                                        cgcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGAGATGCTCCACCAGTTCg                                                            >  1:3301858/1‑150 (MQ=255)
                                                                                         gcgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGa                                                           <  2:3248201/150‑1 (MQ=255)
                                                                                          cgcgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGaa                                                          >  2:4011753/1‑150 (MQ=255)
                                                                                            cgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAAc                                                        <  1:573480/150‑1 (MQ=255)
                                                                                            cgATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAAc                                                        >  1:1482700/1‑150 (MQ=255)
                                                                                             gATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACg                                                       <  1:4614001/150‑1 (MQ=255)
                                                                                               ttttGATGTAGTCATAACGCTCTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCa                                                     <  1:743868/150‑1 (MQ=255)
                                                                                                 ttGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGc                                                   >  2:641209/1‑150 (MQ=255)
                                                                                                  tGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGcc                                                  >  1:1172024/1‑150 (MQ=255)
                                                                                                        aGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGaa                                            <  1:2204389/150‑1 (MQ=255)
                                                                                                         gTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAg                                           >  2:1009772/1‑150 (MQ=255)
                                                                                                           cATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGAc                                         >  1:4222418/1‑150 (MQ=255)
                                                                                                            aTAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACg                                        >  2:4273263/1‑150 (MQ=255)
                                                                                                              aaCGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCt                                      <  2:378260/150‑1 (MQ=255)
                                                                                                                  cacttcacttAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAAt                                  <  2:3717997/150‑1 (MQ=255)
                                                                                                                    cttcacttAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCt                                <  2:2499795/150‑1 (MQ=255)
                                                                                                                    cttcacttAGCGCAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCt                                <  2:4246721/150‑1 (MQ=255)
                                                                                                                       cacttAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGgc                             <  1:2928177/150‑1 (MQ=255)
                                                                                                                         cttAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGgcgc                           >  2:3139764/1‑150 (MQ=255)
                                                                                                                           tAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAg                         >  2:36326/1‑150 (MQ=255)
                                                                                                                           tAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAg                         <  2:2461522/150‑1 (MQ=255)
                                                                                                                           tAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAg                         >  2:2072925/1‑150 (MQ=255)
                                                                                                                            aGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAgg                        <  2:1228889/150‑1 (MQ=255)
                                                                                                                             gcgcAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGa                       >  2:322034/1‑150 (MQ=255)
                                                                                                                             gcgcAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGa                       <  2:2727145/150‑1 (MQ=255)
                                                                                                                               gcAGAGAAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAtt                     >  2:165426/1‑150 (MQ=255)
                                                                                                                                cAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAttt                    <  1:1037967/150‑1 (MQ=255)
                                                                                                                                cAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAttt                    >  1:3952872/1‑150 (MQ=255)
                                                                                                                                 agagAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAtttt                   >  1:2218359/1‑150 (MQ=255)
                                                                                                                                 agagAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGAtttt                   <  2:2737019/150‑1 (MQ=255)
                                                                                                                                   agaAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTg                 <  2:1308443/150‑1 (MQ=255)
                                                                                                                                   agaAATTTGCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTg                 >  1:4334410/1‑150 (MQ=255)
                                                                                                                                    gaAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTgg                <  1:1476430/150‑1 (MQ=255)
                                                                                                                                           gCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGa         >  1:2623032/1‑150 (MQ=255)
                                                                                                                                           gCCATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGa         >  2:2873848/1‑150 (MQ=255)
                                                                                                                                            ccATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGAc        >  1:675610/1‑150 (MQ=255)
                                                                                                                                            ccATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGAc        <  2:1719114/150‑1 (MQ=255)
                                                                                                                                             cATCCCGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACg       >  1:3882974/1‑150 (MQ=255)
                                                                                                                                              aTCCCGGTACCAAAATCATCCAGCCCCAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGc      >  2:1969622/1‑150 (MQ=255)
                                                                                                                                                cccGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAg    >  1:1953934/1‑150 (MQ=255)
                                                                                                                                                cccGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAg    >  1:4498665/1‑150 (MQ=255)
                                                                                                                                                cccGGAACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGCTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAg    >  1:4498578/1‑150 (MQ=255)
                                                                                                                                                  cGGTACCAAACTCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCATTGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAATCGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAggg  <  1:3730746/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATTTGCCATCCCGGTACCAAAATCATCCAGCCACAGCGGGCCAAATTCACACATCGAGGCAAAGGTTGAATCTTTCGGCAGACGGATATGCTCCACCAGTTCGAAACGCAGCCAGGGAAGACGCTCAATCTGGCGCAGGATTTTTGGTTGCTGACGCAGGG  >  CP009273/3671915‑3672210

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: