Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,666,805 A→C 78.1% S197R (AGT→CGT)  ydgD → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,666,8050AC78.1% 852.7 / 177.0 354S197R (AGT→CGT) ydgDputative peptidase
Reads supporting (aligned to +/- strand):  ref base A (47/31);  new base C (139/137);  total (186/168)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 356 total aligned reads displayed.

ATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCC  >  CP009273/1666657‑1666954
                                                                                                                                                    |                                                                                                                                                     
attgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAg                                                                                                                                                      <  2:2421898/150‑1 (MQ=255)
  tgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGtc                                                                                                                                                    <  2:1682566/150‑1 (MQ=255)
  tgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGtc                                                                                                                                                    <  2:2359022/150‑1 (MQ=255)
   tattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGtca                                                                                                                                                   >  2:3498661/2‑150 (MQ=255)
   gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGtca                                                                                                                                                   >  1:4189088/1‑150 (MQ=255)
        tgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCaaa                                                                                                                                              >  2:709130/1‑150 (MQ=255)
         gCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCaaaa                                                                                                                                             <  1:2581302/150‑1 (MQ=255)
          cTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAAc                                                                                                                                            <  2:2095036/150‑1 (MQ=255)
           tACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACt                                                                                                                                           <  1:4138865/150‑1 (MQ=255)
           tACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACt                                                                                                                                           >  1:4257539/1‑150 (MQ=255)
           tACGTAATCCCCCGTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCCGGTCGTACAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATCCGTGGTACCGTCATCAAAACt                                                                                                                                           >  2:2366341/1‑150 (MQ=255)
               tAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGa                                                                                                                                       >  2:369138/1‑150 (MQ=255)
                aaTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGaa                                                                                                                                      >  2:4398165/1‑150 (MQ=255)
                   cccccTTCTGGCATTACGCCGTTGCCGTTATTTGATGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTg                                                                                                                                   <  2:4276191/150‑1 (MQ=255)
                   cccccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTg                                                                                                                                   >  1:3638855/1‑150 (MQ=255)
                    ccccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGa                                                                                                                                  >  1:2270906/1‑150 (MQ=255)
                      ccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGAct                                                                                                                                <  2:2509683/150‑1 (MQ=255)
                      ccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGAct                                                                                                                                <  1:4304185/150‑1 (MQ=255)
                      ccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGAct                                                                                                                                <  2:3960640/150‑1 (MQ=255)
                      ccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGAct                                                                                                                                >  2:75853/1‑150 (MQ=255)
                           tGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGActggctg                                                                                                                           >  2:2228889/1‑150 (MQ=255)
                           tGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGActggctg                                                                                                                           <  2:4594233/150‑1 (MQ=255)
                            gtcATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGActggctgg                                                                                                                          <  2:32350/148‑1 (MQ=255)
                                ttACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGgcgc                                                                                                                      >  2:1992079/1‑150 (MQ=255)
                                   cGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCaaa                                                                                                                   <  1:775423/150‑1 (MQ=255)
                                      cgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGt                                                                                                                >  1:53265/1‑150 (MQ=255)
                                      cgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGt                                                                                                                >  1:2838424/1‑150 (MQ=255)
                                       gttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTc                                                                                                               <  2:4100741/150‑1 (MQ=255)
                                       gttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTc                                                                                                               <  2:135956/150‑1 (MQ=255)
                                        ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                                                                                              <  2:2021331/150‑1 (MQ=255)
                                        ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                                                                                              <  2:3111683/150‑1 (MQ=255)
                                        ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                                                                                              >  1:468046/1‑150 (MQ=255)
                                           ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACGGGCTGGGCGCAAACGTCGGTg                                                                                                           <  2:4835756/150‑1 (MQ=255)
                                            cgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGa                                                                                                          >  2:572580/1‑150 (MQ=255)
                                                  ttGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTcac                                                                                                    >  1:4026960/1‑150 (MQ=255)
                                                     aGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATc                                                                                                 >  2:1254973/1‑150 (MQ=255)
                                                      gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCa                                                                                                >  1:147716/1‑150 (MQ=255)
                                                      gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCa                                                                                                <  1:4940098/150‑1 (MQ=255)
                                                      gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCa                                                                                                >  1:524445/1‑150 (MQ=255)
                                                       ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAg                                                                                               <  1:1678693/150‑1 (MQ=255)
                                                       ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAg                                                                                               <  2:3778458/150‑1 (MQ=255)
                                                          gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGc                                                                                            <  2:4447930/150‑1 (MQ=255)
                                                            tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGa                                                                                          >  1:65608/1‑150 (MQ=255)
                                                            tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGa                                                                                          >  1:2857109/1‑150 (MQ=255)
                                                            tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGa                                                                                          >  2:1150309/1‑150 (MQ=255)
                                                              aaGCCGCGCTTACCGCCGAATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGTGCGCAAACGTCGGTGATGTCACATCAGTGCGATa                                                                                        <  1:3074326/150‑1 (MQ=255)
                                                               aGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATAc                                                                                       >  1:2039894/1‑150 (MQ=255)
                                                                 ccGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGAGGTCACATCAGTGCGATACCt                                                                                     >  2:2778915/1‑150 (MQ=255)
                                                                  cgcgCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCtt                                                                                    <  2:2460451/150‑1 (MQ=255)
                                                                  cgcgCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCtt                                                                                    >  1:329026/1‑150 (MQ=255)
                                                                    cgcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGc                                                                                  >  1:241314/1‑150 (MQ=255)
                                                                      cTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCg                                                                                >  1:5044090/1‑150 (MQ=255)
                                                                      cTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCg                                                                                <  1:3334418/150‑1 (MQ=255)
                                                                       ttACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCgg                                                                               >  1:1980555/1‑150 (MQ=255)
                                                                           cgccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGa                                                                           >  2:651680/1‑150 (MQ=255)
                                                                           cgccgaATTAAAAGCGGCAGGTCGTATAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGa                                                                           <  1:2538743/150‑1 (MQ=255)
                                                                                cATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCg                                                                      >  1:1057298/1‑150 (MQ=255)
                                                                                cATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCg                                                                      >  1:2055852/1‑150 (MQ=255)
                                                                                   tAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGtt                                                                   >  1:4599881/1‑150 (MQ=255)
                                                                                    aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTc                                                                  <  1:1584443/150‑1 (MQ=255)
                                                                                     aaaGCTGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCg                                                                 <  1:3498517/150‑1 (MQ=255)
                                                                                         cGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCctc                                                             <  2:4599856/150‑1 (MQ=255)
                                                                                           gCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTg                                                           >  2:1084875/1‑150 (MQ=255)
                                                                                            cAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGa                                                          >  1:5074592/1‑150 (MQ=255)
                                                                                             aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGAt                                                         >  1:1297520/1‑150 (MQ=255)
                                                                                             aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGAt                                                         >  1:3877785/1‑150 (MQ=255)
                                                                                              ggTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATg                                                        >  1:1142274/1‑150 (MQ=255)
                                                                                                 cGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTg                                                     <  1:711652/150‑1 (MQ=255)
                                                                                                  gTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGc                                                    >  2:3568229/1‑150 (MQ=255)
                                                                                                      aGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATAc                                                <  1:99170/150‑1 (MQ=255)
                                                                                                       gTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATAcc                                               <  2:2847172/150‑1 (MQ=255)
                                                                                                        tGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCg                                              >  2:4489921/1‑150 (MQ=255)
                                                                                                        tGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCg                                              >  1:4507360/1‑150 (MQ=255)
                                                                                                         gACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGa                                             <  2:3293952/150‑1 (MQ=255)
                                                                                                         gACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGa                                             <  2:3868803/150‑1 (MQ=255)
                                                                                                             caggcaggCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGAc                                         <  2:929599/150‑1 (MQ=255)
                                                                                                                gcaggcTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGc                                      <  2:4736611/150‑1 (MQ=255)
                                                                                                                gcaggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCAGCTGATGTTGCATACCTATGACGGc                                      >  2:254399/1‑150 (MQ=255)
                                                                                                                   ggcTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTTGATGGGCGCAAACTTAGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTgg                                   <  1:2174008/150‑1 (MQ=255)
                                                                                                                   ggcTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTgg                                   >  1:3717246/1‑150 (MQ=255)
                                                                                                                     cTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCa                                 <  2:3633421/150‑1 (MQ=255)
                                                                                                                      tACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaa                                >  2:2297512/1‑150 (MQ=255)
                                                                                                                         ccTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaatta                             <  1:181224/150‑1 (MQ=255)
                                                                                                                          cTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaa                            >  2:1580127/1‑150 (MQ=255)
                                                                                                                           tGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaat                           >  1:234936/1‑150 (MQ=255)
                                                                                                                            gAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaatt                          <  2:1492519/150‑1 (MQ=255)
                                                                                                                             aaGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTg                         >  1:2729298/1‑150 (MQ=255)
                                                                                                                               gATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTggg                       >  1:2990405/1‑150 (MQ=255)
                                                                                                                                atcatcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTgggg                      <  2:381908/150‑1 (MQ=255)
                                                                                                                                 tcatcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGt                     <  1:1213905/150‑1 (MQ=255)
                                                                                                                                    tctcGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCa                  <  2:1859167/150‑1 (MQ=255)
                                                                                                                                     ctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCaa                 >  1:921804/1‑150 (MQ=255)
                                                                                                                                           aCGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCg           <  2:3393587/150‑1 (MQ=255)
                                                                                                                                            cGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCgg          >  2:3152170/1‑150 (MQ=255)
                                                                                                                                             gTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc         >  2:1897917/1‑150 (MQ=255)
                                                                                                                                              ttGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCt        >  2:5001742/1‑150 (MQ=255)
                                                                                                                                                gTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTcc      <  2:817271/150‑1 (MQ=255)
                                                                                                                                                 tACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCt     <  1:4893837/150‑1 (MQ=255)
                                                                                                                                                 tACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCt     >  2:848757/1‑150 (MQ=255)
                                                                                                                                                    aGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGcc  <  2:2659252/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCC  >  CP009273/1666657‑1666954

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: ‑