Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 49,029 A→C 10.0% S421R (AGC→CGC)  kefC → potassium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927349,0290AC10.0% 728.4 / 79.6 279S421R (AGC→CGC) kefCpotassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (112/139);  new base C (7/21);  total (119/160)
Fisher's exact test for biased strand distribution p-value = 6.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 285 total aligned reads displayed.

CGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCC  >  CP009273/48880‑49178
                                                                                                                                                     |                                                                                                                                                     
cGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCa                                                                                                                                                       <  1:3841167/150‑1 (MQ=255)
      ttCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCgggg                                                                                                                                                 >  2:1145311/1‑150 (MQ=255)
      ttCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCgggg                                                                                                                                                 >  1:2918431/1‑150 (MQ=255)
       tCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGt                                                                                                                                                <  2:2409319/150‑1 (MQ=255)
         tgctgGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGa                                                                                                                                              >  2:2067267/1‑150 (MQ=255)
                gATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAAtggt                                                                                                                                       >  1:2261828/1‑150 (MQ=255)
                  tCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAAtggtgg                                                                                                                                     <  1:4228705/150‑1 (MQ=255)
                  tCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAAtggtgg                                                                                                                                     <  1:3173531/150‑1 (MQ=255)
                  tCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCACGTGAAGCCGATTAGATCTACTAAGAACAGCCGCGCGTGATTATCGCCGTATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGATCTCCCGCTGGGTGAAAAtgttgg                                                                                                                                     >  2:323985/1‑150 (MQ=255)
                  tCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCACGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCCGCGGGGTGAAAAtggtgg                                                                                                                                     >  2:4370248/1‑150 (MQ=255)
                   ccTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATggtggt                                                                                                                                    <  1:4600049/150‑1 (MQ=255)
                     tCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTAc                                                                                                                                  <  2:3195038/150‑1 (MQ=255)
                        aTCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCg                                                                                                                               <  2:2452618/150‑1 (MQ=255)
                         tCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGa                                                                                                                              >  2:2561488/1‑150 (MQ=255)
                           gCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCCGCGGGGTGAAAATGGTGGTACTCGATc                                                                                                                            <  1:2022671/150‑1 (MQ=255)
                             cTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTTCTCTCCAGCGGGGTGAAAATGGTTGTACTCGATCAc                                                                                                                          >  2:447383/1‑150 (MQ=255)
                             cTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCAc                                                                                                                          <  1:9892/150‑1 (MQ=255)
                              ttGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACg                                                                                                                         <  2:1917050/150‑1 (MQ=255)
                                  tcAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGTCGAAGAACAGCCGCGCGTGATTATCGCCTGATTCGGTCGTTTTGGGCAGATTACCGTACGTTTACTGCTCTCCCTCGGGGTGAAAATGGTGGTACTCGATCACGATcc                                                                                                                     <  2:594934/149‑1 (MQ=255)
                                  gCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATcc                                                                                                                     <  1:342367/150‑1 (MQ=255)
                                    aaTCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCgg                                                                                                                   <  1:703914/150‑1 (MQ=255)
                                    aaTCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCgg                                                                                                                   <  2:706754/150‑1 (MQ=255)
                                      tcttctACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGAc                                                                                                                 <  1:2449865/150‑1 (MQ=255)
                                      tcttctACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGAc                                                                                                                 <  2:3598093/150‑1 (MQ=255)
                                       cttctACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGAcc                                                                                                                <  2:4580944/150‑1 (MQ=255)
                                          ctACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCata                                                                                                             >  1:1728492/1‑150 (MQ=255)
                                          ctACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACAGCTCTCCCGCGGGGTGAAAATGGTGGTACGCGATCACGATCCGGACCata                                                                                                             >  1:2744882/1‑150 (MQ=255)
                                           tACTGATGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCatat                                                                                                            <  1:1657747/150‑1 (MQ=255)
                                            aCTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATc                                                                                                           <  2:3369857/150‑1 (MQ=255)
                                            aCTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATc                                                                                                           <  1:1866899/150‑1 (MQ=255)
                                             cTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCg                                                                                                          <  2:1277268/150‑1 (MQ=255)
                                              tGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGa                                                                                                         >  2:1812181/1‑150 (MQ=255)
                                                aGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGaaa                                                                                                       >  1:2824771/1‑150 (MQ=255)
                                                 ggAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAAc                                                                                                      >  1:3267397/1‑150 (MQ=255)
                                                       gcgTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCg                                                                                                <  2:1303997/150‑1 (MQ=255)
                                                        cgTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGt                                                                                               <  2:1843994/150‑1 (MQ=255)
                                                           gAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTaaa                                                                                            <  2:3000642/150‑1 (MQ=255)
                                                            aaGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAAt                                                                                           <  2:3156824/150‑1 (MQ=255)
                                                            aaGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAAt                                                                                           <  2:11861/150‑1 (MQ=255)
                                                              gCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAAttt                                                                                         <  1:1491549/150‑1 (MQ=255)
                                                              gCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAAttt                                                                                         <  2:2867977/150‑1 (MQ=255)
                                                                cGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTGTTGGGCAGATTACCGGACGTTTACTTCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTgg                                                                                       >  1:2546464/1‑150 (MQ=255)
                                                                  aTGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTa                                                                                     >  1:1023924/1‑150 (MQ=255)
                                                                     agaTCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGa                                                                                  <  1:92728/150‑1 (MQ=255)
                                                                     agaTCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGa                                                                                  >  2:2166695/1‑150 (MQ=255)
                                                                     agaTCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGa                                                                                  <  1:3156507/150‑1 (MQ=255)
                                                                      gaTCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTAATCGATCACGATCCGGACCATATCGAAACCTTGCGGAAATTTTGTATGaa                                                                                 >  1:1164392/1‑150 (MQ=255)
                                                                       aTCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGaaa                                                                                <  2:2693797/150‑1 (MQ=255)
                                                                          gacgaAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATTGTGGTACTCGATCACGATCCGGACCATATCGAAACCATGCGTAAATTTGGTATGAAAgtg                                                                             >  1:401582/1‑150 (MQ=255)
                                                                          gacgaAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAgtg                                                                             >  1:2665362/1‑150 (MQ=255)
                                                                           acgaAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAgtgt                                                                            >  2:3347263/1‑150 (MQ=255)
                                                                           acgaAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAgtgt                                                                            >  2:4101873/1‑150 (MQ=255)
                                                                            cgaAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGtt                                                                           <  2:834838/150‑1 (MQ=255)
                                                                             gaagaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGttt                                                                          <  2:113198/150‑1 (MQ=255)
                                                                             gaagaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGttt                                                                          <  1:2965641/150‑1 (MQ=255)
                                                                             gaagaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGATCTCCAGCGGGGTGAAAATGGTGGTACTCGATAACGATCCGGACCATATCGAAACCATGCGTAAATTTGGTATGAAAGTtttt                                                                          >  1:4402374/1‑150 (MQ=255)
                                                                               agaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTa                                                                        >  1:680269/1‑150 (MQ=255)
                                                                               agaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTa                                                                        >  2:3506768/1‑150 (MQ=255)
                                                                               agaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTa                                                                        <  2:3538552/150‑1 (MQ=255)
                                                                                gaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTAt                                                                       >  2:2727531/1‑150 (MQ=255)
                                                                                gaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTAt                                                                       >  1:3460150/1‑150 (MQ=255)
                                                                                gaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTAt                                                                       <  2:3762206/150‑1 (MQ=255)
                                                                                 aaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATg                                                                      <  1:2312443/150‑1 (MQ=255)
                                                                                  aCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCCGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATgg                                                                     <  1:3826148/150‑1 (MQ=255)
                                                                                  aCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATgg                                                                     <  1:1957888/150‑1 (MQ=255)
                                                                                  aCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATgg                                                                     <  2:211117/150‑1 (MQ=255)
                                                                                       cgcgcgTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATg                                                                <  2:3482908/150‑1 (MQ=255)
                                                                                        gcgcgTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGc                                                               <  2:4001705/150‑1 (MQ=255)
                                                                                           cgTGATTCTCGCCGGATTCGGTCGGTTTGGGCAGATTACCGGACCTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACACGATCACGATCCGGACCATATCGAAACCTTGCGTAAAGTTGGTATGAAAGTGTTTTATGTCgcatccac                                                            >  2:1542058/1‑143 (MQ=255)
                                                                                                  aTCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCggatg                                                     >  2:760059/1‑150 (MQ=255)
                                                                                                   tCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCggatgg                                                    >  2:4635643/1‑150 (MQ=255)
                                                                                                    cGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCCGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCggatgga                                                   <  1:3994363/150‑1 (MQ=255)
                                                                                                    cGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCCGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCggatgga                                                   <  1:412438/150‑1 (MQ=255)
                                                                                                       cGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGAttt                                                <  2:1338065/150‑1 (MQ=255)
                                                                                                        ggATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTa                                               >  2:4240163/1‑150 (MQ=255)
                                                                                                          aTTCGGTCGTTTTGGGCAGATTACCGGACGTTTAATTCTCTCCAGCGGGGTGAAAAAGGTGGTACACGATAACCATCCGGACCAAATCGAAACCATTCGTAAATTTTGTATTAAATTGTTTTATTGCGATTCCACGCGGATGGATTTaat                                             >  1:3917420/1‑148 (MQ=255)
                                                                                                           ttagggcGTTTTGGGGACTTTAACGGTAGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTg                                            <  1:1516124/144‑1 (MQ=255)
                                                                                                             cGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTgg                                           >  2:957353/1‑149 (MQ=255)
                                                                                                             cGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGa                                          <  1:3198126/150‑1 (MQ=255)
                                                                                                             cGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCATGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGa                                          >  2:3219563/1‑150 (MQ=255)
                                                                                                              ggTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGaa                                         >  1:3030911/1‑150 (MQ=255)
                                                                                                                    tttGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGcc                                   <  1:2511745/150‑1 (MQ=255)
                                                                                                                      tGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCgg                                 >  1:3129479/1‑150 (MQ=255)
                                                                                                                       gggCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGa                                <  2:680269/150‑1 (MQ=255)
                                                                                                                         gCAGATTACCGGACGTTTACTGCTCTCCCGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAgc                              <  2:1637478/150‑1 (MQ=255)
                                                                                                                          cAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAgcg                             >  2:391825/1‑150 (MQ=255)
                                                                                                                           aGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGACTCTGCCGGAgcgg                            >  1:42851/1‑150 (MQ=255)
                                                                                                                             aTTACCGGACGTTTACTTCTATCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACATTGCGTAAATTTGGTATGAAAGTGTTTTAATGCGATGCCACGCGGATGGATTTACTGGAATATGCCGGAgcggcg                          >  1:971449/1‑150 (MQ=255)
                                                                                                                              ttACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGa                         <  1:3677253/150‑1 (MQ=255)
                                                                                                                               tACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGaa                        >  2:628796/1‑150 (MQ=255)
                                                                                                                                      cGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGaa                 <  2:2777353/150‑1 (MQ=255)
                                                                                                                                      cGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGaa                 <  2:4156588/150‑1 (MQ=255)
                                                                                                                                        tttACTGCTCTCCCGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGt               <  2:309223/150‑1 (MQ=255)
                                                                                                                                            cTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTg           >  1:3968331/1‑150 (MQ=255)
                                                                                                                                            cTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTg           <  1:1260507/150‑1 (MQ=255)
                                                                                                                                            cTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTg           <  1:1099570/150‑1 (MQ=255)
                                                                                                                                             tGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGa          >  2:1998820/1‑150 (MQ=255)
                                                                                                                                             tGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGa          >  1:1718255/1‑150 (MQ=255)
                                                                                                                                               ctctCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGAtt        <  2:2942312/150‑1 (MQ=255)
                                                                                                                                               ctctCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGAtt        >  1:1362404/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGGGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCC  >  CP009273/48880‑49178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: