Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 628,732 G→A 10.0% R68C (CGC→TGC)  ybdH ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273628,7320GA10.0% 909.7 / 69.7 271R68C (CGC→TGC) ybdHputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (128/116);  new base A (15/12);  total (143/128)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 280 total aligned reads displayed.

CGGAACGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGAT  >  CP009273/628583‑628881
                                                                                                                                                     |                                                                                                                                                     
cgGAACGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACcgcg                                                                                                                                                       <  1:2696417/150‑1 (MQ=255)
cgGAACGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACcgcg                                                                                                                                                       <  1:98579/150‑1 (MQ=255)
  gAACGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGaa                                                                                                                                                     >  2:1388169/1‑150 (MQ=255)
   aaCGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGaaa                                                                                                                                                    >  1:4323821/1‑150 (MQ=255)
    aCGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAAc                                                                                                                                                   >  1:2119709/1‑150 (MQ=255)
       gCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACaaa                                                                                                                                                >  1:1080197/1‑150 (MQ=255)
       gCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACaaa                                                                                                                                                <  2:2916273/150‑1 (MQ=255)
        caacaaACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACaaaa                                                                                                                                               >  1:2574047/1‑150 (MQ=255)
        caacaaACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACaaaa                                                                                                                                               >  1:3332913/1‑150 (MQ=255)
         aacaaACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAat                                                                                                                                              <  1:5126556/150‑1 (MQ=255)
         aacaaACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAat                                                                                                                                              <  2:531353/150‑1 (MQ=255)
           caaACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAatat                                                                                                                                            <  2:1752665/150‑1 (MQ=255)
             aaCGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGc                                                                                                                                          >  1:1874403/1‑150 (MQ=255)
              aCGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCt                                                                                                                                         >  1:696889/1‑150 (MQ=255)
               cGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCtt                                                                                                                                        <  2:1132149/150‑1 (MQ=255)
                ggCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCttt                                                                                                                                       <  2:2636987/150‑1 (MQ=255)
                  cAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGc                                                                                                                                     >  1:2439239/1‑150 (MQ=255)
                  cAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGc                                                                                                                                     <  1:1614608/150‑1 (MQ=255)
                   aGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGcc                                                                                                                                    >  2:4829140/1‑150 (MQ=255)
                   aGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGcc                                                                                                                                    >  2:3047666/1‑150 (MQ=255)
                   aGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTGGTACATCGCTTTCGCTGCAATGACCGCGAACCAAAATATGCTGTGcc                                                                                                                                    >  2:422686/1‑150 (MQ=255)
                       cGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTg                                                                                                                                >  1:3384388/1‑150 (MQ=255)
                       cGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTg                                                                                                                                >  2:2666971/1‑150 (MQ=255)
                        gagaCGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTgg                                                                                                                               >  2:3308690/1‑150 (MQ=255)
                               cgggcgAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCa                                                                                                                        <  1:4759477/150‑1 (MQ=255)
                               cgggcgAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTGTCGCTGCAATGACCGCGACACAAAATATGCTTTGCCCCTGGCAGTCCa                                                                                                                        >  2:2452882/1‑150 (MQ=255)
                                 ggcgAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCaaa                                                                                                                      <  1:3217015/150‑1 (MQ=255)
                                   cgAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACg                                                                                                                    <  2:491981/150‑1 (MQ=255)
                                    gAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGc                                                                                                                   >  2:2731309/1‑150 (MQ=255)
                                        gCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGc                                                                                                               <  1:2379394/150‑1 (MQ=255)
                                         cTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCg                                                                                                              >  2:4421557/1‑150 (MQ=255)
                                           ttCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGa                                                                                                            >  1:111323/1‑150 (MQ=255)
                                           ttCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGa                                                                                                            <  2:1101770/150‑1 (MQ=255)
                                           ttCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGAACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCAGGCAGTCCAAACGCTGGCGGa                                                                                                            >  1:352455/1‑150 (MQ=255)
                                            tCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGaa                                                                                                           >  2:1696658/1‑150 (MQ=255)
                                            tCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGaa                                                                                                           >  1:3844485/1‑150 (MQ=255)
                                                ggTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAgttt                                                                                                       <  1:3281829/150‑1 (MQ=255)
                                                    tCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAgtttggtt                                                                                                   >  1:4844432/1‑150 (MQ=255)
                                                     cGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAgtttggttt                                                                                                  >  1:3997966/1‑150 (MQ=255)
                                                       aGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTgc                                                                                                <  2:2053358/150‑1 (MQ=255)
                                                       aGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTgc                                                                                                <  1:4112075/150‑1 (MQ=255)
                                                        gCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTgcg                                                                                               <  2:2748760/150‑1 (MQ=255)
                                                          aGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGcgcc                                                                                             >  2:216061/1‑150 (MQ=255)
                                                            tGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGcgccga                                                                                           <  2:3555791/150‑2 (MQ=255)
                                                                 cGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCACTGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGACGTTTGGTTTGCGCCGCAGCaa                                                                                      >  2:737170/1‑150 (MQ=255)
                                                                 cGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGACGTTTGGTTTGCGCCGCAGCaa                                                                                      >  2:4952427/1‑150 (MQ=255)
                                                                 cGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCaa                                                                                      >  2:1933581/1‑150 (MQ=255)
                                                                  gCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAAt                                                                                     <  2:557173/150‑1 (MQ=255)
                                                                  gCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAAt                                                                                     >  2:882894/1‑150 (MQ=255)
                                                                   ccACCGACGCCAATCACCCCGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATCTGCTGTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATg                                                                                    <  2:1887589/150‑1 (MQ=255)
                                                                   ccACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATg                                                                                    >  1:5911/1‑150 (MQ=255)
                                                                    cACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGACGTTTGGTTTGCGCCGCAGCAATgg                                                                                   >  1:3097625/1‑150 (MQ=255)
                                                                    cACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATgg                                                                                   >  2:3872006/1‑150 (MQ=255)
                                                                     aCCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGgc                                                                                  <  2:969810/150‑1 (MQ=255)
                                                                      ccGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGgcg                                                                                 >  2:952733/1‑150 (MQ=255)
                                                                        gACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                                               >  2:4046262/1‑150 (MQ=255)
                                                                        gACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGcgcg                                                                               >  2:2824048/1‑150 (MQ=255)
                                                                              aaTCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGcc                                                                         >  2:2056698/1‑150 (MQ=255)
                                                                               aTCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACACAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGGTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCg                                                                        >  1:859383/1‑150 (MQ=255)
                                                                                tCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTGTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGt                                                                       >  1:3666635/1‑150 (MQ=255)
                                                                                         tGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGTTTTGCGCCGCAGCAATGGCGCGGTTGCCGTAGATCcaca                                                              >  2:818949/1‑150 (MQ=255)
                                                                                            ggTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcg                                                           >  2:4132581/1‑150 (MQ=255)
                                                                                             gtcgtcACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgc                                                          >  1:390983/1‑150 (MQ=255)
                                                                                             gtcgtcACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgc                                                          >  1:391137/1‑150 (MQ=255)
                                                                                                gtcACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACcgcgcg                                                       >  1:1564575/1‑150 (MQ=255)
                                                                                                  cACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCgag                                                     <  2:1344566/150‑1 (MQ=255)
                                                                                                   aCCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCgaga                                                    >  2:2515258/1‑150 (MQ=255)
                                                                                                    ccGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGaa                                                   >  2:2200547/1‑150 (MQ=255)
                                                                                                     cGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGaaa                                                  <  1:3739014/150‑1 (MQ=255)
                                                                                                        aCTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgtt                                               <  1:3146514/150‑1 (MQ=255)
                                                                                                        aCTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgtt                                               >  2:613574/1‑150 (MQ=255)
                                                                                                          tCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttgt                                             >  2:4734299/1‑150 (MQ=255)
                                                                                                          tCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAgttgt                                             <  2:4609787/150‑1 (MQ=255)
                                                                                                            aGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTtc                                           <  2:5140461/150‑1 (MQ=255)
                                                                                                             gccgccAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTtca                                          <  2:1638626/150‑1 (MQ=255)
                                                                                                                 ccAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                                      <  1:1814606/150‑1 (MQ=255)
                                                                                                                 ccAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAg                                      >  1:2368814/1‑150 (MQ=255)
                                                                                                                  cAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGt                                     <  1:5093307/150‑1 (MQ=255)
                                                                                                                  cAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGt                                     <  1:129872/150‑1 (MQ=255)
                                                                                                                    gttgttGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGa                                   <  1:4360025/150‑1 (MQ=255)
                                                                                                                     ttgttgTACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaa                                  <  2:43181/150‑1 (MQ=255)
                                                                                                                       gttgtACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACCCCGCGCGAGAAAGTTGTTCATCAGTGaaaa                                >  2:1575350/1‑150 (MQ=255)
                                                                                                                        ttgtACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGaaaaa                               >  1:2904064/1‑150 (MQ=255)
                                                                                                                          gtACATCGCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAAt                             <  2:1245690/150‑1 (MQ=255)
                                                                                                                           tACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATc                            <  1:2256403/150‑1 (MQ=255)
                                                                                                                           tACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATc                            <  1:4697917/150‑1 (MQ=255)
                                                                                                                           tACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATc                            >  1:1983653/1‑150 (MQ=255)
                                                                                                                            aCATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCg                           >  1:2595827/1‑150 (MQ=255)
                                                                                                                                cGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCa                       <  1:4213560/150‑1 (MQ=255)
                                                                                                                                 gCTTTCGCTGCAATGACCGCAAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAg                      >  1:736981/1‑150 (MQ=255)
                                                                                                                                  cTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAgg                     <  2:3499411/150‑1 (MQ=255)
                                                                                                                                      cGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGAt                 >  1:2685152/1‑150 (MQ=255)
                                                                                                                                       gCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGAtt                <  2:3098057/150‑1 (MQ=255)
                                                                                                                                         tGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGa              <  2:696889/150‑1 (MQ=255)
                                                                                                                                          gCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGaa             <  1:5058031/150‑1 (MQ=255)
                                                                                                                                          gCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGaa             <  1:2964270/150‑1 (MQ=255)
                                                                                                                                            aaTGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAAc           <  1:2800429/150‑1 (MQ=255)
                                                                                                                                                 ccGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCt      >  2:3647024/1‑150 (MQ=255)
                                                                                                                                                  cgcgAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTg     <  2:3704852/150‑1 (MQ=255)
                                                                                                                                                  cgcgAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTg     <  2:2382682/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGAACGGCAACAAACGGCAGACCGAGACGGCGGGCGAGGGCTTTCGCGGTGTCGAGCAGTGCACCGCCACCGACGCCAATCACCACGCTGCGGTCGTCACCGGACTCAGCCGCCAGTTGTTGTACATCGCTTTCGCTGCAATGACCGCGAAACAAAATATGCTTTGCCCCTGGCAGTCCAAACGCTGGCGGAAGTTTGGTTTGCGCCGCAGCAATGGCGCGTTTGCCGTAGATCCACACCGCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGAT  >  CP009273/628583‑628881

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: