Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,653,563 C→A 5.0% S413* (TCG→TAG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,653,5630CA 5.0% 1001.5 / 19.9 258S413* (TCG→TAG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (111/134);  new base A (8/5);  total (119/139)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 260 total aligned reads displayed.

AATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATG  >  CP009273/1653414‑1653712
                                                                                                                                                     |                                                                                                                                                     
 aTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCg                                                                                                                                                      <  1:4625189/150‑1 (MQ=255)
     gCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCt                                                                                                                                                  <  2:1168368/150‑1 (MQ=255)
          gCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtacc                                                                                                                                             <  1:2877909/150‑1 (MQ=255)
            aaaaCCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtaccat                                                                                                                                           >  1:2680050/1‑150 (MQ=255)
            aaaaCCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtaccat                                                                                                                                           <  1:666325/150‑1 (MQ=255)
            aaaaCCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtaccat                                                                                                                                           <  1:666361/150‑1 (MQ=255)
              aaCCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtaccatta                                                                                                                                         <  1:3961267/150‑1 (MQ=255)
               aCCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtaccattac                                                                                                                                        <  2:3248623/150‑1 (MQ=255)
                cccGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtaccattacc                                                                                                                                       >  2:2613594/1‑150 (MQ=255)
                   gCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATa                                                                                                                                    <  1:2762275/150‑1 (MQ=255)
                   gCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATa                                                                                                                                    <  1:4173998/150‑1 (MQ=255)
                   gCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATa                                                                                                                                    <  1:788696/150‑1 (MQ=255)
                      tataTCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATAGACCTATACCATTACCATAGaa                                                                                                                                 >  1:3897496/1‑150 (MQ=255)
                       atatCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAAc                                                                                                                                >  1:2848712/1‑150 (MQ=255)
                         atCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGc                                                                                                                              >  2:4546386/1‑150 (MQ=255)
                           cTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCt                                                                                                                            <  1:1018959/150‑1 (MQ=255)
                           cTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCt                                                                                                                            <  1:1130601/150‑1 (MQ=255)
                            tGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTg                                                                                                                           >  1:2594513/1‑150 (MQ=255)
                            tGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTg                                                                                                                           <  2:613918/150‑1 (MQ=255)
                             gCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGc                                                                                                                          >  2:1873767/1‑150 (MQ=255)
                                 gggCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGt                                                                                                                      <  2:2086849/150‑1 (MQ=255)
                                  ggCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTg                                                                                                                     <  2:1679785/150‑1 (MQ=255)
                                    cTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCt                                                                                                                   <  2:2776291/150‑1 (MQ=255)
                                     tGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTg                                                                                                                  <  1:1580179/150‑1 (MQ=255)
                                       ggggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGa                                                                                                                >  1:232347/1‑150 (MQ=255)
                                       ggggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGa                                                                                                                >  2:2515400/1‑150 (MQ=255)
                                          gCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATaa                                                                                                             >  2:697595/1‑150 (MQ=255)
                                           ccACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAAt                                                                                                            <  2:2869904/150‑1 (MQ=255)
                                           ccACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAAt                                                                                                            <  2:3896802/150‑1 (MQ=255)
                                             acaACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATcc                                                                                                          >  1:1091018/1‑150 (MQ=255)
                                               aaCGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCgg                                                                                                        >  1:1293856/1‑150 (MQ=255)
                                               aaCGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCgg                                                                                                        <  1:1395127/150‑1 (MQ=255)
                                               aaCGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCgg                                                                                                        >  1:4061364/1‑150 (MQ=255)
                                                aCGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGt                                                                                                       <  2:1419569/150‑1 (MQ=255)
                                                aCGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGt                                                                                                       <  2:2560176/150‑1 (MQ=255)
                                                  gCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCa                                                                                                     <  1:2078599/150‑1 (MQ=255)
                                                   ccAGGCATACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGTCGGGCAATGTCGTCATCAGCGGCGGAAATAGTGGCGCTCGTGAATATACCTATACCATTACCATAGATCGCCTGCCGGTGCTGGATAATCCGGTCaa                                                                                                    <  2:3263858/150‑1 (MQ=255)
                                                     aGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCaaaa                                                                                                  <  2:457663/150‑1 (MQ=255)
                                                      ggCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAAc                                                                                                 <  2:731966/150‑1 (MQ=255)
                                                         aaaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTc                                                                                              >  1:120863/1‑150 (MQ=255)
                                                          aaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCa                                                                                             <  2:2295861/150‑1 (MQ=255)
                                                          aaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCa                                                                                             <  1:1830418/150‑1 (MQ=255)
                                                          aaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCa                                                                                             <  1:1438583/150‑1 (MQ=255)
                                                          aaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCa                                                                                             <  2:1056007/150‑1 (MQ=255)
                                                            cGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAt                                                                                           >  1:890433/1‑150 (MQ=255)
                                                            cGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATAGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAt                                                                                           <  1:1345953/150‑1 (MQ=255)
                                                              gCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAtct                                                                                         <  2:290082/150‑1 (MQ=255)
                                                               cGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAtctc                                                                                        <  1:2163047/150‑1 (MQ=255)
                                                               cGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAtctc                                                                                        >  1:651200/1‑150 (MQ=255)
                                                                gAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTcc                                                                                       >  1:835736/1‑150 (MQ=255)
                                                                 aaCTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCt                                                                                      >  1:2725620/1‑150 (MQ=255)
                                                                    tgactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCt                                                                                   <  1:3245631/150‑1 (MQ=255)
                                                                    tgactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATAGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCt                                                                                   >  2:4232153/1‑150 (MQ=255)
                                                                     gactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCtt                                                                                  <  1:1646552/150‑1 (MQ=255)
                                                                          cACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCt                                                                             <  1:2646637/150‑1 (MQ=255)
                                                                          cACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCt                                                                             <  1:2252571/150‑1 (MQ=255)
                                                                              cgcTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGAc                                                                         <  1:2472105/150‑1 (MQ=255)
                                                                                     gCTATGTTACCTATTTTGACGGGCAATGTCGGCATCATCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCt                                                                  >  2:826704/1‑150 (MQ=255)
                                                                                     gCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCt                                                                  <  1:953991/150‑1 (MQ=255)
                                                                                            tACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATc                                                           >  2:4063059/1‑150 (MQ=255)
                                                                                               cTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATg                                                        >  1:754497/1‑150 (MQ=255)
                                                                                               cTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATg                                                        >  1:3374987/1‑150 (MQ=255)
                                                                                               cTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATg                                                        <  1:3709950/150‑1 (MQ=255)
                                                                                                tATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATgg                                                       >  1:3798689/1‑150 (MQ=255)
                                                                                                   tttGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTcc                                                    >  1:279694/1‑150 (MQ=255)
                                                                                                      gACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCa                                                 >  1:4553565/1‑150 (MQ=255)
                                                                                                       aCGGGCAATGTCGGCATCAGCGGCGGAACTAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCaa                                                >  2:4592511/1‑150 (MQ=255)
                                                                                                         gggCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAAt                                              <  1:4370937/150‑1 (MQ=255)
                                                                                                         gggCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAAt                                              <  1:2696270/150‑1 (MQ=255)
                                                                                                            cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGAc                                           >  2:3383184/1‑150 (MQ=255)
                                                                                                            cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGAc                                           >  1:2010453/1‑150 (MQ=255)
                                                                                                            cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGAc                                           >  1:4248154/1‑150 (MQ=255)
                                                                                                             aaTGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACg                                          >  2:1635378/1‑150 (MQ=255)
                                                                                                             aaTGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACg                                          <  2:4567219/150‑1 (MQ=255)
                                                                                                                 tCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCaa                                      >  2:4014811/1‑150 (MQ=255)
                                                                                                                    gCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATcc                                   >  1:548371/1‑150 (MQ=255)
                                                                                                                      aTCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCACTCcgc                                 >  2:542178/1‑150 (MQ=255)
                                                                                                                      aTCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCcgc                                 <  1:2369723/150‑1 (MQ=255)
                                                                                                                       tCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCcgcg                                >  2:4178884/1‑150 (MQ=255)
                                                                                                                         aGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGAc                              <  1:3532689/150‑1 (MQ=255)
                                                                                                                          gcggcgGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACg                             <  2:3588088/150‑1 (MQ=255)
                                                                                                                            ggcggAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGc                           >  1:3269539/1‑150 (MQ=255)
                                                                                                                              cggACATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCACATGACGGCAATCCGCGACGGCGt                         >  2:2394056/1‑150 (MQ=255)
                                                                                                                              cggAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGt                         >  2:2367437/1‑150 (MQ=255)
                                                                                                                               ggAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTc                        <  1:4463476/150‑1 (MQ=255)
                                                                                                                               ggAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTc                        >  1:4544986/1‑150 (MQ=255)
                                                                                                                               ggAAATAGTGGCGCGCGTGAATAGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTc                        >  1:3837311/1‑150 (MQ=255)
                                                                                                                                gAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTcc                       <  2:879124/150‑1 (MQ=255)
                                                                                                                                 aaaTAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcg                      <  2:126102/150‑1 (MQ=255)
                                                                                                                                  aaTAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgc                     <  2:1317109/150‑1 (MQ=255)
                                                                                                                                   aTAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgcg                    <  2:342836/150‑1 (MQ=255)
                                                                                                                                   aTAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgcg                    >  1:1858948/1‑150 (MQ=255)
                                                                                                                                       tGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCac                >  2:2767380/1‑149 (MQ=255)
                                                                                                                                        ggCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCaaa               <  1:2631467/150‑1 (MQ=255)
                                                                                                                                           gcgcgTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGAt            >  2:2108049/1‑150 (MQ=255)
                                                                                                                                              cgTGAATAGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATaaa         >  1:3429202/1‑150 (MQ=255)
                                                                                                                                               gTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAAc        <  1:3023288/150‑1 (MQ=255)
                                                                                                                                                 gAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTg      <  2:1922963/150‑1 (MQ=255)
                                                                                                                                                 gAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTg      <  1:77844/150‑1 (MQ=255)
                                                                                                                                                  aaTCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTgg     <  2:3791352/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATG  >  CP009273/1653414‑1653712

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: