Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,321,041 G→T 9.0% F194L (TTC→TTA yfaP ← DUF2135 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,321,0410GT 9.0% 1285.5 / 65.7 345F194L (TTC→TTAyfaPDUF2135 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base G (151/163);  new base T (16/15);  total (167/178)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 348 total aligned reads displayed.

ATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAG  >  CP009273/2320892‑2321190
                                                                                                                                                     |                                                                                                                                                     
   ggTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGaaa                                                                                                                                                    >  1:2982907/1‑150 (MQ=255)
    gTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGaaaa                                                                                                                                                   >  1:160994/1‑150 (MQ=255)
     tACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAAt                                                                                                                                                  <  2:4938891/150‑1 (MQ=255)
     tACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATAAACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAAt                                                                                                                                                  >  1:4333037/1‑150 (MQ=255)
           aaaTGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCggg                                                                                                                                            >  1:5088592/1‑150 (MQ=255)
           aaaTGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCggg                                                                                                                                            >  1:5034534/1‑150 (MQ=255)
           aaaTGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCggg                                                                                                                                            >  2:2056775/1‑150 (MQ=255)
             aTGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTc                                                                                                                                          <  1:3914397/150‑1 (MQ=255)
             aTGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTc                                                                                                                                          <  1:720827/150‑1 (MQ=255)
              tGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTcc                                                                                                                                         <  2:3740825/150‑1 (MQ=255)
               gTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCg                                                                                                                                        <  1:204203/150‑1 (MQ=255)
                tttCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGt                                                                                                                                       <  1:4787626/150‑1 (MQ=255)
                tttCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGt                                                                                                                                       >  2:2867196/1‑150 (MQ=255)
                 ataatGTTTCTCATTGACCGATCATTCATCGGTGATCATCGTCAGTTGGGCGGTGGTTAATTCATATTCGCTTCTCCAAACATATTAGATGATATACACCAGATATATCCCGAGAAATGTCTCGGGAATTGCGAAAATCTCGGGTCCGTa                                                                                                                                      >  2:2280081/5‑150 (MQ=255)
                  tCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTAc                                                                                                                                     >  1:119305/1‑150 (MQ=255)
                  tCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTAc                                                                                                                                     >  1:245577/1‑150 (MQ=255)
                   ccTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTAcc                                                                                                                                    <  1:4534683/150‑1 (MQ=255)
                     tGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTAcccc                                                                                                                                  >  2:1372995/1‑150 (MQ=255)
                       tttCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCcgt                                                                                                                                >  1:4196582/1‑150 (MQ=255)
                        ttCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCcgtc                                                                                                                               <  1:4951043/150‑1 (MQ=255)
                         tctcATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCcgtcg                                                                                                                              >  2:3033730/1‑150 (MQ=255)
                          ctcATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCcgtcgt                                                                                                                             >  2:2144386/1‑150 (MQ=255)
                            cATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGa                                                                                                                           >  2:4234630/1‑150 (MQ=255)
                             aTTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGAc                                                                                                                          >  1:1837523/1‑150 (MQ=255)
                               tGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACAt                                                                                                                        >  2:1678020/1‑150 (MQ=255)
                                  ccGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCa                                                                                                                     >  1:2701989/1‑150 (MQ=255)
                                  ccGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCa                                                                                                                     <  2:2212559/150‑1 (MQ=255)
                                       ccTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATcc                                                                                                                >  1:990693/1‑150 (MQ=255)
                                         ttCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAg                                                                                                              <  1:3882750/150‑1 (MQ=255)
                                           cATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTg                                                                                                            >  1:4633463/1‑150 (MQ=255)
                                             tCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGACTTGGCGCGGGCCTGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACCTCCATATCCAGTGca                                                                                                          >  2:2607202/1‑150 (MQ=255)
                                               ggTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGcacc                                                                                                        >  1:1098791/1‑150 (MQ=255)
                                                    tCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGt                                                                                                   >  2:110534/1‑150 (MQ=255)
                                                    tCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGt                                                                                                   <  1:3263769/150‑1 (MQ=255)
                                                      aGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGttt                                                                                                 <  2:1526425/150‑1 (MQ=255)
                                                      aGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGttt                                                                                                 <  1:3039809/150‑1 (MQ=255)
                                                        cGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGttttt                                                                                               <  1:3442359/150‑1 (MQ=255)
                                                         gTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTc                                                                                              >  2:802045/1‑150 (MQ=255)
                                                          tCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCa                                                                                             >  2:3576243/1‑150 (MQ=255)
                                                          tCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCa                                                                                             <  1:838743/150‑1 (MQ=255)
                                                          tCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCa                                                                                             <  1:1489817/150‑1 (MQ=255)
                                                            aGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGc                                                                                           >  1:873483/1‑150 (MQ=255)
                                                              ttGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCAc                                                                                         >  1:2310067/1‑150 (MQ=255)
                                                               tGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACg                                                                                        <  1:2133160/150‑1 (MQ=255)
                                                                  gCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGgtg                                                                                     >  2:3642594/1‑150 (MQ=255)
                                                                      tggtTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTAc                                                                                 >  1:4466298/1‑150 (MQ=255)
                                                                         ttAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGt                                                                              <  1:433223/150‑1 (MQ=255)
                                                                           aaTTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTAc                                                                            >  2:59815/1‑150 (MQ=255)
                                                                           aaTTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTAc                                                                            >  2:2771930/1‑150 (MQ=255)
                                                                                 gTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCg                                                                      >  1:672636/1‑150 (MQ=255)
                                                                                   ttCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCgtg                                                                    <  2:5034534/150‑1 (MQ=255)
                                                                                    tCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCgtgt                                                                   >  2:3985492/1‑150 (MQ=255)
                                                                                    tCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCgtgt                                                                   >  1:3140187/1‑150 (MQ=255)
                                                                                          cgcCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCa                                                             <  1:2280081/150‑1 (MQ=255)
                                                                                          cgcCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCa                                                             >  1:3365526/1‑150 (MQ=255)
                                                                                           gcCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCAt                                                            >  2:520378/1‑150 (MQ=255)
                                                                                             ccaccaTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCg                                                          >  1:3493270/1‑150 (MQ=255)
                                                                                              caccaTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCgg                                                         >  2:2926130/1‑150 (MQ=255)
                                                                                              caccaTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCgg                                                         <  1:1386381/150‑1 (MQ=255)
                                                                                                ccaTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGc                                                       <  1:2862167/150‑1 (MQ=255)
                                                                                                 caTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCg                                                      >  2:4849387/1‑150 (MQ=255)
                                                                                                  aTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGt                                                     >  2:4964889/1‑150 (MQ=255)
                                                                                                  aTAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGt                                                     >  2:917698/1‑150 (MQ=255)
                                                                                                  aTAGTAGTTGATATACACCTGATAGCGCCCGTGAATGGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGt                                                     <  2:3355494/150‑1 (MQ=255)
                                                                                                    agtagtTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTaa                                                   >  1:2233433/1‑150 (MQ=255)
                                                                                                      tagtTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTaaca                                                 <  1:1254596/150‑1 (MQ=255)
                                                                                                        gtTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTaacaac                                               <  1:4110158/150‑1 (MQ=255)
                                                                                                          tGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGt                                             <  2:4430181/150‑1 (MQ=255)
                                                                                                            atatACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGa                                           >  1:966528/1‑150 (MQ=255)
                                                                                                            atatACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGa                                           <  1:4769339/150‑1 (MQ=255)
                                                                                                             tataCACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGaa                                          <  2:4783585/150‑1 (MQ=255)
                                                                                                             tataCACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGaa                                          <  2:3468965/150‑1 (MQ=255)
                                                                                                               taCACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCTGGCGTGTTCGCCATCGGGCGTAACAACGTGAagg                                        <  2:719869/150‑1 (MQ=255)
                                                                                                                 cacCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtc                                      <  2:2559535/150‑1 (MQ=255)
                                                                                                                  acCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcg                                     >  1:3652339/1‑150 (MQ=255)
                                                                                                                  acCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcg                                     >  2:187401/1‑150 (MQ=255)
                                                                                                                   ccTGATAGCGCCCGTGAATTGGCGCGGGCATGGCTAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcga                                    >  1:1813973/1‑150 (MQ=255)
                                                                                                                   ccTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcga                                    >  2:1545331/1‑150 (MQ=255)
                                                                                                                   ccTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGGGAaggtcga                                    >  2:2769454/1‑150 (MQ=255)
                                                                                                                    cTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCCCTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcgag                                   <  2:4223429/150‑1 (MQ=255)
                                                                                                                     tGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcgagg                                  <  2:3111233/150‑1 (MQ=255)
                                                                                                                       aTAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcgaggtc                                <  2:4667314/150‑1 (MQ=255)
                                                                                                                       aTAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcgaggtc                                <  1:3429064/150‑1 (MQ=255)
                                                                                                                       aTAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAaggtcgaggtc                                >  2:4363683/1‑150 (MQ=255)
                                                                                                                            gcCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTAt                           >  1:1071169/1‑150 (MQ=255)
                                                                                                                             cccGTGAAGTGGCGCGGGCAGGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTAtt                          <  2:4588330/150‑1 (MQ=255)
                                                                                                                              ccGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTg                         >  1:3175009/1‑150 (MQ=255)
                                                                                                                               cGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGt                        >  2:3233683/1‑150 (MQ=255)
                                                                                                                                  gAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCg                     >  1:3661373/1‑150 (MQ=255)
                                                                                                                                     ttGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCTCGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTAt                  <  2:2309905/150‑1 (MQ=255)
                                                                                                                                         cgcgGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCa              >  1:3735355/1‑150 (MQ=255)
                                                                                                                                         cgcgGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCa              >  2:1013472/1‑150 (MQ=255)
                                                                                                                                          gcgGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCAc             <  1:1150272/150‑1 (MQ=255)
                                                                                                                                           cgGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACg            >  1:3606952/1‑150 (MQ=255)
                                                                                                                                            gggCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACga           >  2:2107848/1‑150 (MQ=255)
                                                                                                                                                    cGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCa   <  2:3958986/150‑1 (MQ=255)
                                                                                                                                                    cGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCa   <  1:2813715/150‑1 (MQ=255)
                                                                                                                                                    cGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCa   <  1:919938/150‑1 (MQ=255)
                                                                                                                                                    cGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCa   <  2:1337230/150‑1 (MQ=255)
                                                                                                                                                     gAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAg  >  2:3876396/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGATCAGCGTCAGTTGGGCGGTGGTTAATTCCGTTTCGCTGCGCCCACCATAGTAGTTGATATACACCTGATAGCGCCCGTGAATTGGCGCGGGCATGGCGAAAATCTCGGGTCCGTACCCCGTCGTGACATCCATATCCAGTGCACCACTGTTTTTCAGCACGGTGTTACCGTACCAGGCGTGTTCGCCATCGGGCGTAACAACGTGAAGGTCGAGGTCGGTATTGTCCGTATCCCACGAGAGAACCAG  >  CP009273/2320892‑2321190

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: