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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,666,8050AC64.2% 514.0 / inf 355S197R (AGT→CGT) ydgDputative peptidase
Reads supporting (aligned to +/- strand):  ref base A (68/59);  new base C (117/111);  total (185/170)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 357 total aligned reads displayed.

TGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCC  >  CP009273/1666659‑1666954
                                                                                                                                                  |                                                                                                                                                     
 gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGtca                                                                                                                                                   <  2:2128604/150‑1 (MQ=255)
    tgtgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGtcatca                                                                                                                                                >  1:3457029/1‑150 (MQ=255)
     gtgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCaa                                                                                                                                               <  2:3901728/150‑1 (MQ=255)
      tgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCaac                                                                                                                                              >  1:5001024/1‑149 (MQ=255)
      tgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCaaa                                                                                                                                              <  1:588881/150‑1 (MQ=255)
        cTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAAc                                                                                                                                            <  2:3175970/150‑1 (MQ=255)
         tACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACt                                                                                                                                           >  2:3853518/1‑150 (MQ=255)
           cGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACtgt                                                                                                                                         >  1:3153492/1‑150 (MQ=255)
             tAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGa                                                                                                                                       >  2:2896011/1‑150 (MQ=255)
              aaTCCCCCTTCTGGCATTCCGCCGTGGCCTTTTTTTGTGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGaa                                                                                                                                      <  1:4611291/150‑1 (MQ=255)
              aaTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGaa                                                                                                                                      >  2:694434/1‑150 (MQ=255)
                  ccccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGAGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTTACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGa                                                                                                                                  <  2:3517331/150‑1 (MQ=255)
                     cTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGActg                                                                                                                               >  1:3991072/1‑150 (MQ=255)
                     cTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActg                                                                                                                               >  1:1928873/1‑150 (MQ=255)
                         tGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGActggctg                                                                                                                           >  1:3839767/1‑150 (MQ=255)
                         tGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGActggctg                                                                                                                           <  1:3105888/150‑1 (MQ=255)
                          ggCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActggctgg                                                                                                                          >  2:4211527/1‑150 (MQ=255)
                            cATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGgc                                                                                                                        >  2:4609210/1‑150 (MQ=255)
                            cATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGgc                                                                                                                        >  2:2790776/1‑150 (MQ=255)
                                  gccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAAc                                                                                                                  <  2:4888267/150‑1 (MQ=255)
                                   ccgttgcctttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACg                                                                                                                 <  1:464348/150‑1 (MQ=255)
                                           gttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGAt                                                                                                         >  1:2470737/1‑150 (MQ=255)
                                           gttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGAGACTGGCTGGGCGCACACGTCGGTGAt                                                                                                         >  2:269928/1‑150 (MQ=255)
                                           gttATTGGCGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGAt                                                                                                         <  1:3885446/150‑1 (MQ=255)
                                            ttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATg                                                                                                        >  1:2744319/1‑150 (MQ=255)
                                                     ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAg                                                                                               >  1:4250824/1‑150 (MQ=255)
                                                       agaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTg                                                                                             >  1:4864634/1‑150 (MQ=255)
                                                           aaaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGAt                                                                                         >  2:3098900/1‑150 (MQ=255)
                                                           aaaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGAt                                                                                         >  2:3755697/1‑150 (MQ=255)
                                                              gCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATAcc                                                                                      <  2:4011058/150‑1 (MQ=255)
                                                               ccGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCt                                                                                     <  1:2967404/150‑1 (MQ=255)
                                                               ccGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCt                                                                                     <  2:3440367/150‑1 (MQ=255)
                                                               ccGCGCTTACCGCCGAATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCt                                                                                     <  1:779214/150‑1 (MQ=255)
                                                                cgcgCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCtt                                                                                    <  1:2384552/150‑1 (MQ=255)
                                                                 gcgcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTg                                                                                   <  1:3775355/150‑1 (MQ=255)
                                                                  cgcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGc                                                                                  >  1:2040922/1‑150 (MQ=255)
                                                                  agcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGc                                                                                  <  2:1311998/149‑1 (MQ=255)
                                                                   gcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGcc                                                                                 <  1:2156473/150‑1 (MQ=255)
                                                                   gcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGcc                                                                                 <  2:4120287/150‑1 (MQ=255)
                                                                   gcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGcc                                                                                 <  2:1210360/150‑1 (MQ=255)
                                                                   gcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGcc                                                                                 <  2:4528990/150‑1 (MQ=255)
                                                                     ttACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGAGCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCgg                                                                               >  2:1889873/1‑150 (MQ=255)
                                                                      tACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCggg                                                                              <  2:2041149/150‑1 (MQ=255)
                                                                       aCCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGt                                                                             <  2:3029916/150‑1 (MQ=255)
                                                                       aCCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGt                                                                             <  1:1011764/150‑1 (MQ=255)
                                                                        ccgccgCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTg                                                                            <  2:5002631/150‑1 (MQ=255)
                                                                           ccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACa                                                                         <  1:1372495/150‑1 (MQ=255)
                                                                            cgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAg                                                                        >  2:4348041/1‑150 (MQ=255)
                                                                             gcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGc                                                                       >  2:4400977/1‑150 (MQ=255)
                                                                             gcATTAAAAGCGGCAGGTCGTAAAGTGAATCAGGCAGGCTACCCTGAAGATCATCTCTATACTTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGc                                                                       >  2:748361/1‑150 (MQ=255)
                                                                              cATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCg                                                                      <  2:394033/150‑1 (MQ=255)
                                                                               aTTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCTGTGCGATACCTTGCCGGGTGACAGCgg                                                                     >  2:3693614/1‑150 (MQ=255)
                                                                                ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGt                                                                    >  2:1664515/1‑150 (MQ=255)
                                                                                 tAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGtt                                                                   <  2:1038724/150‑1 (MQ=255)
                                                                                  aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTc                                                                  >  1:4813324/1‑150 (MQ=255)
                                                                                  aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTc                                                                  <  2:4072104/150‑1 (MQ=255)
                                                                                   aaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCg                                                                 >  2:1765306/1‑150 (MQ=255)
                                                                                    aaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGc                                                                >  1:5167296/1‑150 (MQ=255)
                                                                                      gCGGCAGGTCGTAAAGTGACTCAGGCTGGCTACCGTGAAGATCCTGGCGCTCCGTTGGACAGTCCTCAACACTGTGATGTGACTGGCTGGGCGCAAACGTCGGTGCTGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCct                                                              <  1:1989551/150‑1 (MQ=255)
                                                                                        ggCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGTGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCctct                                                            <  2:381891/150‑1 (MQ=255)
                                                                                          cAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGa                                                          <  2:2767055/150‑1 (MQ=255)
                                                                                           aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGAt                                                         <  1:1869953/150‑1 (MQ=255)
                                                                                           aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGAt                                                         <  1:941312/150‑1 (MQ=255)
                                                                                           aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGAt                                                         <  1:2253707/150‑1 (MQ=255)
                                                                                              tCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGtt                                                      <  2:5183818/150‑1 (MQ=255)
                                                                                               cGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTg                                                     >  1:1753381/1‑150 (MQ=255)
                                                                                                   aaGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATa                                                 <  1:4457065/150‑1 (MQ=255)
                                                                                                         cTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATg                                           >  2:3487682/1‑150 (MQ=255)
                                                                                                          tCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGAGGTTCGCCTCTGATGTTGCATACCGATGa                                          >  2:2505360/1‑150 (MQ=255)
                                                                                                           caggcaggCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGAc                                         <  2:4663076/150‑1 (MQ=255)
                                                                                                             ggcaggcTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACgg                                       >  2:5221775/1‑150 (MQ=255)
                                                                                                              gcaggcTACCCAGAAGATCATCTCGATACGTTGTACAGTCATCAACACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACgcc                                      <  1:4533767/150‑3 (MQ=255)
                                                                                                               caggcTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCt                                     >  2:2286723/1‑150 (MQ=255)
                                                                                                               caggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCt                                     <  2:4561497/150‑1 (MQ=255)
                                                                                                                aggcTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTg                                    <  1:4718001/150‑1 (MQ=255)
                                                                                                                aggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTg                                    >  2:5157843/1‑150 (MQ=255)
                                                                                                                      cccTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaatt                              >  1:1050102/1‑150 (MQ=255)
                                                                                                                       ccTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaatta                             <  2:3385224/150‑1 (MQ=255)
                                                                                                                       ccTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaatta                             >  2:2413118/1‑150 (MQ=255)
                                                                                                                       ccTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaatta                             >  1:5312807/1‑150 (MQ=255)
                                                                                                                        cTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaa                            >  2:1613783/1‑150 (MQ=255)
                                                                                                                        cTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaa                            >  2:3321520/1‑150 (MQ=255)
                                                                                                                            aGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTgg                        <  2:2476182/150‑1 (MQ=255)
                                                                                                                             gATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCAGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTTGCAATTAATTggg                       >  1:2991128/1‑150 (MQ=255)
                                                                                                                              atcatcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTgggg                      >  2:819694/1‑150 (MQ=255)
                                                                                                                              atcatcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTgggg                      >  2:2843403/1‑150 (MQ=255)
                                                                                                                              atcatcTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTgggg                      <  2:2151221/150‑1 (MQ=255)
                                                                                                                               tcatcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGt                     >  1:682320/1‑150 (MQ=255)
                                                                                                                                catcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTg                    <  2:4300618/150‑1 (MQ=255)
                                                                                                                                catcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTg                    <  2:4813858/150‑1 (MQ=255)
                                                                                                                                catcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTg                    >  1:2641810/1‑150 (MQ=255)
                                                                                                                                catcTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTg                    >  1:985787/1‑150 (MQ=255)
                                                                                                                                 atcTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCTGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAAGTCATTGGGGTGc                   >  1:3407578/1‑150 (MQ=255)
                                                                                                                                     cGATACGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAg               <  1:812388/150‑1 (MQ=255)
                                                                                                                                       aTACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGtt             <  1:367920/150‑1 (MQ=255)
                                                                                                                                          cGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCgg          >  1:735646/1‑150 (MQ=255)
                                                                                                                                          cGTTGTACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCgg          >  1:1209903/1‑150 (MQ=255)
                                                                                                                                            ttGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCt        <  1:3216948/150‑1 (MQ=255)
                                                                                                                                               tACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCTGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCt     <  2:2964376/150‑1 (MQ=255)
                                                                                                                                               tACCGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCt     <  2:1386119/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCC  >  CP009273/1666659‑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: