Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 819,200 G→T 15.4% H252N (CAT→AAT)  ybhP ← endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273819,2000GT15.4% 805.7 / 88.1 277H252N (CAT→AAT) ybhPendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base G (133/101);  new base T (29/14);  total (162/115)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 278 total aligned reads displayed.

CGACGTCATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTC  >  CP009273/819054‑819349
                                                                                                                                                  |                                                                                                                                                     
   cGTCATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAAtctc                                                                                                                                                 >  1:2263201/1‑150 (MQ=255)
    gTCATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTcc                                                                                                                                                <  2:3271059/150‑1 (MQ=255)
    gTCATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTcc                                                                                                                                                >  1:1424644/1‑150 (MQ=255)
      cATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGc                                                                                                                                              >  1:738479/1‑150 (MQ=255)
      cATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGc                                                                                                                                              >  2:5203327/1‑150 (MQ=255)
        tCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCAc                                                                                                                                            <  1:3581990/150‑1 (MQ=255)
          cTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACtt                                                                                                                                          <  1:920360/150‑1 (MQ=255)
              aaCCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAg                                                                                                                                      >  1:1580009/1‑150 (MQ=255)
               aCCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAgg                                                                                                                                     >  1:994281/1‑150 (MQ=255)
                  aGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGACAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGc                                                                                                                                  >  2:3926300/1‑150 (MQ=255)
                   gATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCa                                                                                                                                 >  2:776978/1‑150 (MQ=255)
                     tAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATg                                                                                                                               <  1:2367518/150‑1 (MQ=255)
                       aaCGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCCTGAt                                                                                                                             >  2:2068083/1‑150 (MQ=255)
                         cGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCa                                                                                                                           <  2:2817025/150‑1 (MQ=255)
                          gTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAg                                                                                                                          >  2:2869107/1‑150 (MQ=255)
                          gTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAg                                                                                                                          >  2:4631535/1‑150 (MQ=255)
                              cAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAg                                                                                                                      >  1:4663622/1‑150 (MQ=255)
                              cAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAg                                                                                                                      >  1:2971646/1‑150 (MQ=255)
                               aaGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCACACTTAAAGGGGCATGATCAGAAAgg                                                                                                                     >  2:3044033/1‑150 (MQ=255)
                                aGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGgt                                                                                                                    >  1:58850/1‑150 (MQ=255)
                                 gAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGgtg                                                                                                                   <  1:3431066/150‑1 (MQ=255)
                                 gAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGgtg                                                                                                                   >  1:4551292/1‑150 (MQ=255)
                                 gAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGgtg                                                                                                                   <  1:405376/150‑1 (MQ=255)
                                   atgatgCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTc                                                                                                                 <  1:694555/150‑1 (MQ=255)
                                   atgatgCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTc                                                                                                                 <  2:3418685/150‑1 (MQ=255)
                                    tgatgCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCg                                                                                                                <  1:4280472/150‑1 (MQ=255)
                                    tgatgCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCg                                                                                                                <  1:2890733/150‑1 (MQ=255)
                                     gatgCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGc                                                                                                               >  1:5121306/1‑150 (MQ=255)
                                       ttctttCTTGTGCCTCGCAAATCGCCTTAAACACCGCGGGATAATATTGCTCTCCGTTTTACAGCAACTGGATCTTATTGCCTTCGCGCAATCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGAATGATCAGAAAGGTTTCTCCa                                                                                                             <  1:3624763/146‑1 (MQ=255)
                                          gTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGt                                                                                                          >  1:5255798/1‑150 (MQ=255)
                                          gTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGt                                                                                                          <  1:2431488/150‑1 (MQ=255)
                                          gTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGt                                                                                                          >  2:322660/1‑150 (MQ=255)
                                           ttCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTc                                                                                                         <  2:2296803/150‑1 (MQ=255)
                                            tCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTcc                                                                                                        <  2:1044069/150‑1 (MQ=255)
                                            tCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTcc                                                                                                        <  2:2501562/150‑1 (MQ=255)
                                            tCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTcc                                                                                                        <  2:5146292/150‑1 (MQ=255)
                                               tgtcCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGTAGTCGCCATGTCCTCa                                                                                                     >  1:3271284/1‑150 (MQ=255)
                                                gtgCCTCGCCAATCGACTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTTGATCTTATTGACTTCGCGCCAGCTACATTTCATAAATGAATCTCCTCACTTAAAGGGGCATGATAAGAAAGGTGTCGCCATGTCCGCAg                                                                                                    >  2:3848846/1‑150 (MQ=255)
                                                 tgCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGc                                                                                                   >  1:4928613/1‑150 (MQ=255)
                                                      cGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCaa                                                                                              >  2:3365018/1‑150 (MQ=255)
                                                      cGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCaa                                                                                              <  1:2605760/150‑1 (MQ=255)
                                                      cGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCaa                                                                                              >  1:1833555/1‑150 (MQ=255)
                                                       gCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAAc                                                                                             >  2:1915664/1‑150 (MQ=255)
                                                       gCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAAc                                                                                             <  1:3850944/150‑1 (MQ=255)
                                                         cAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAAcgc                                                                                           >  2:3857593/1‑150 (MQ=255)
                                                         cAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAAcgc                                                                                           >  1:3780554/1‑150 (MQ=255)
                                                         cAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAAcgc                                                                                           >  1:2683187/1‑150 (MQ=255)
                                                          aaTCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAAcgcg                                                                                          >  1:4835677/1‑150 (MQ=255)
                                                          aaTCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTATCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCCTGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCACcgcg                                                                                          >  1:2031564/1‑150 (MQ=255)
                                                           aTCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCgg                                                                                         >  2:554154/1‑150 (MQ=255)
                                                              gCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTg                                                                                      <  2:2137683/150‑1 (MQ=255)
                                                                           cgGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCgcgctggcgct                                                                         >  2:675038/1‑150 (MQ=255)
                                                                            gggATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCgcgctggcgctg                                                                        >  1:4029046/1‑150 (MQ=255)
                                                                            gggATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCgcgctggcgctg                                                                        <  1:4154935/150‑1 (MQ=255)
                                                                              gATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCgctggcgctggc                                                                      >  2:3225415/1‑150 (MQ=255)
                                                                                taataTTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACGTAAAGGGGCATGATCAGAAAGGTGTCTCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCAt                                                                    >  2:2699853/1‑150 (MQ=255)
                                                                                  atatTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGcgctggcgctgg                                                                       >  2:1989246/1‑145 (MQ=255)
                                                                                  atatTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCAttt                                                                  <  1:5223904/150‑1 (MQ=255)
                                                                                       gCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGAc                                                             <  1:824384/150‑1 (MQ=255)
                                                                                       gCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGAc                                                             >  2:3738543/1‑150 (MQ=255)
                                                                                        cTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACg                                                            <  2:2450711/150‑1 (MQ=255)
                                                                                        cTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACg                                                            <  2:887039/150‑1 (MQ=255)
                                                                                         tCGCCGTTTTCCAGCAACTGGATCTTAGTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGt                                                           <  2:1487774/150‑1 (MQ=255)
                                                                                           gCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAg                                                         <  2:404878/150‑1 (MQ=255)
                                                                                            ccGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGa                                                        >  2:3285494/1‑150 (MQ=255)
                                                                                              gTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATc                                                      <  2:2263201/150‑1 (MQ=255)
                                                                                              gTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATc                                                      >  1:659843/1‑150 (MQ=255)
                                                                                               ttttCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATcc                                                     >  2:4740205/1‑150 (MQ=255)
                                                                                               tttgCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACTTTTCATAAATGTATCTCCGCACTTATAGGGTCTTGATCAGATTGGTGTCGCCATGTGCGCAGCGGCATCTCGGTTGGCGCGCTGGCGCTGGCATTTTTGTAGTTGATcc                                                     <  2:2908667/150‑1 (MQ=255)
                                                                                               atttCCAGCACCTGGATGTTATAGCCGACGCGCCAGCTACATTTCATAACTGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATcc                                                     <  2:4357018/149‑1 (MQ=255)
                                                                                                 ttCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTg                                                   <  2:5295759/150‑1 (MQ=255)
                                                                                                 ttCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTg                                                   <  2:4333055/150‑1 (MQ=255)
                                                                                                    cAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTcc                                                <  1:1991415/150‑1 (MQ=255)
                                                                                                     aGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCa                                               >  2:4977983/1‑150 (MQ=255)
                                                                                                     aGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCa                                               <  2:4957067/150‑1 (MQ=255)
                                                                                                        aaCTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTc                                            >  2:1140967/1‑150 (MQ=255)
                                                                                                         aCTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCg                                           >  2:4437493/1‑150 (MQ=255)
                                                                                                         aCTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCg                                           >  1:1724553/1‑150 (MQ=255)
                                                                                                          cTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCgt                                          <  2:2674888/150‑1 (MQ=255)
                                                                                                           tGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCgta                                         >  2:4224052/1‑150 (MQ=255)
                                                                                                           tGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCgta                                         <  1:1162732/150‑1 (MQ=255)
                                                                                                                      gCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAAt                              >  1:1063247/1‑150 (MQ=255)
                                                                                                                       ccTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGACAtt                             >  2:1268009/1‑150 (MQ=255)
                                                                                                                        cTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTg                            <  1:4874139/150‑1 (MQ=255)
                                                                                                                                 aGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGaaa                   >  1:3938451/1‑150 (MQ=255)
                                                                                                                                  gCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAAc                  >  1:770406/1‑150 (MQ=255)
                                                                                                                                   cTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACg                 >  1:4367675/1‑150 (MQ=255)
                                                                                                                                   cTACATTTCATAAATGAATCTCCACACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACg                 <  1:3044033/150‑1 (MQ=255)
                                                                                                                                     aCATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTg               <  2:1580009/150‑1 (MQ=255)
                                                                                                                                     aCATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTg               >  2:2974723/1‑150 (MQ=255)
                                                                                                                                     aCATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTg               >  2:3148674/1‑150 (MQ=255)
                                                                                                                                     aCATTTCATAAATGAATCTCCGCACTGAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTg               >  2:4927283/1‑150 (MQ=255)
                                                                                                                                      cATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTgc              >  1:4989901/1‑150 (MQ=255)
                                                                                                                                        tttCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTgcgc            >  1:3143241/1‑150 (MQ=255)
                                                                                                                                        tttCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTgcgc            <  1:3365026/150‑1 (MQ=255)
                                                                                                                                            aTAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTGCGCGCCg        <  1:4320517/150‑1 (MQ=255)
                                                                                                                                              aaaTTAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATGTTTGACGTAGATCCTGTCCCGTCGGAGTAGAGGACAGTGCACCGGAAACGTGCGCGCCggg      >  1:4102698/1‑150 (MQ=255)
                                                                                                                                              aaaTGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTGCGCGCCggg      >  1:3550231/1‑150 (MQ=255)
                                                                                                                                               aaTGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTGCGCGCCGGGc     <  2:1789707/150‑1 (MQ=255)
                                                                                                                                                aTGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTGCGCGCCGGGCg    >  1:3664370/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGACGTCATCCTCAAACCAGATAAACGTTTCAAGAATGATGCGTTCTTGTGCCTCGCCAATCGCCTTAAACACCGCGGGATAATATTGCTCGCCGTTTTCCAGCAACTGGATCTTATTGCCTTCGCGCCAGCTACATTTCATAAATGAATCTCCGCACTTAAAGGGGCATGATCAGAAAGGTGTCGCCATGTCCGCAGCGGCAACGCGGTTGGCGCGCTGGCGCTGGCATTTTTGACGTAGATCCTGTCCAGTCGTAGTAGAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTC  >  CP009273/819054‑819349

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: