Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,471 A→G 34.5% S27S (AGT→AGC)  tpn ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4710AG34.5% 125.3 / 104.3 130S27S (AGT→AGC) tpnIS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (48/37);  new base G (20/25);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 135 total aligned reads displayed.

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACC  >  CP000730/36331‑36607
                                                                                                                                            |                                                                                                                                         
cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAACTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACTCCACGTTTCTTTAATTTTTTAGATATATCACGGTaa                                                                                                                                           >  1:80076/1‑141 (MQ=11)
cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                           <  1:547506/141‑1 (MQ=18)
cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                           <  2:531112/141‑1 (MQ=21)
 cccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                          <  2:765457/140‑1 (MQ=40)
  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAAGTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  1:416904/1‑141 (MQ=11)
  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATTACGCTAGCt                                                                                                                                         >  2:573776/1‑141 (MQ=11)
  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  1:565915/1‑141 (MQ=11)
  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCGGATATTTTACGGTaatt                                                                                                                                         >  1:264366/1‑139 (MQ=11)
  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAcc                                                                                                                                         >  2:202201/1‑140 (MQ=21)
  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCCTGAACCCAACGGTAGATCGTTGCATGATGAACGTGTAGAGCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  2:716410/1‑141 (MQ=11)
   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATAATATTGGGGCTTATTCGTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACta                                                                                                                                        >  2:479688/1‑140 (MQ=11)
   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGGTAGCta                                                                                                                                        >  1:496066/1‑140 (MQ=11)
    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                       <  1:283539/141‑1 (MQ=11)
    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCCCTTTATATTTCAGATATATCACGATAGCTCa                                                                                                                                       <  2:437676/141‑1 (MQ=255)
    aTTTGTAATAAGCATATTTATGCTTATACTTCCAAAATTGATATAAAATTGGGGCATATTCTTGAACCAAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                       <  1:638641/141‑1 (MQ=11)
     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGCTATATCACGATAGCGCaa                                                                                                                                      >  2:269438/1‑141 (MQ=11)
     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  1:260230/141‑1 (MQ=11)
     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  1:257465/141‑1 (MQ=21)
     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:633660/141‑1 (MQ=21)
     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:540448/141‑1 (MQ=21)
     tttGTAATAAGCTTTCTTATGCTTTTTCTCCCAAATTTGATATATACCTGCGCCTCAGACTCGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  1:219437/141‑1 (MQ=11)
     tttCTAATAAGCTTTTTCTTGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTCGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:528253/141‑1 (MQ=11)
     gttGTAATAAGCTGTATTATGCTTGTTCTACCAAATTTGATACAAAATTGGGGCATATTCAGGAACCCAACGATAGACCGATGAATGATGAAGGTTTACACACCGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      >  2:353346/2‑141 (MQ=11)
      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                     >  1:307533/1‑141 (MQ=255)
      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                     <  1:648721/141‑1 (MQ=255)
      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                     >  1:289400/1‑141 (MQ=40)
      ttGTAATAAGCTTTTTTATGCTTTTTCTCCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                     <  2:55673/141‑1 (MQ=11)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGATATAAAATTGGGGCATATTTTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    >  2:167028/1‑141 (MQ=255)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGTACGTTTACACCACGTTCTTTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  1:595859/1‑141 (MQ=11)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTATACCACGTTCTCTTAATATTTCAGATATATTACGATAGCTCAATg                                                                                                                                    <  1:431441/141‑1 (MQ=11)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  1:786455/1‑141 (MQ=11)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTGTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  1:105735/1‑141 (MQ=11)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    >  1:144420/1‑141 (MQ=25)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    <  2:690831/141‑1 (MQ=18)
       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    <  1:30276/141‑1 (MQ=39)
        ttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                   <  1:233503/140‑2 (MQ=11)
        ttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                   <  2:225225/140‑2 (MQ=11)
        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                   <  1:312020/141‑2 (MQ=11)
        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                                                   >  2:642231/1‑140 (MQ=18)
         ttataaGCTTTTTTATGCTTTTTCTGTCAAATTTGATATAAAATGGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGCTAGCTCAATgca                                                                                                                                  <  2:214525/139‑3 (MQ=11)
         taataaGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATAATCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                  >  1:389432/1‑139 (MQ=255)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCACTgca                                                                                                                                  >  1:635242/1‑139 (MQ=11)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  <  1:544595/141‑3 (MQ=11)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  >  2:59428/1‑139 (MQ=11)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  >  2:423638/1‑139 (MQ=11)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  >  2:103200/1‑139 (MQ=255)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACtgcctgac                                                                                                                                  >  1:369917/1‑134 (MQ=255)
         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtt                                                                                                                                  >  2:535678/1‑140 (MQ=255)
          aataaGCTTTTTTTTTTTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACGATAACTCAAttcat                                                                                                                                 >  2:549609/1‑137 (MQ=11)
          aataaGCTTTTTTATTTTTTTTCTTCCAAATTTGATATAAAATTGTGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAtccat                                                                                                                                 >  2:635205/1‑137 (MQ=11)
          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACTCCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 <  2:129491/141‑4 (MQ=11)
          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 <  2:47209/141‑4 (MQ=11)
          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 <  2:147434/141‑4 (MQ=11)
          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGATCTCTTAATATTTCAGATATATCACGATAGCTCAATgcaa                                                                                                                                 >  2:685558/1‑138 (MQ=11)
          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGGACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgcat                                                                                                                                 >  1:495088/1‑138 (MQ=2)
          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCATTGtaa                                                                                                                                 >  2:268896/1‑140 (MQ=255)
           ataaGCTTTTTTATGCTTTATATTGCACATTTGATATAAACTTGGGGGAGATTAGAGAACCCGAAGATGGATCGGTGAGTGATGAACGTTTAAGAGACGTTATCTTAATATTTCAGATACATCACGATAGCGCaagatatt                                                                                                                                >  2:204848/1‑135 (MQ=255)
                     ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACAGCACGTTCACTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                      >  1:582588/1‑141 (MQ=25)
                     ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                      >  2:337639/1‑141 (MQ=35)
                      tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:96053/1‑141 (MQ=35)
                      tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:467483/1‑141 (MQ=35)
                      tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACAGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:512139/1‑141 (MQ=25)
                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                    >  2:332236/1‑141 (MQ=35)
                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                    >  1:94844/1‑141 (MQ=35)
                         gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                  >  1:564426/1‑141 (MQ=35)
                           tcttgcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCa                                                                                                                <  1:241332/136‑1 (MQ=25)
                                   ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTTTATCTTAGATAGTAGCCAACGGCTAc                                                                                                        >  1:599114/1‑141 (MQ=21)
                                    ttttttttATACAAAATTGGGGCATAATCTTGTACCCAACGGTAGACCGTTGAATGGTGAACGTTTACACCACTTTTTTTTATTATTTCAGATATATGATGATAACTCAATTTATATTTTAGATAGTAGCCAACgggtgaa                                                                                                       >  2:495916/5‑136 (MQ=11)
                                           gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGAAGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAGGGCTACAGTgctca                                                                                                >  1:521911/1‑137 (MQ=255)
                                                  aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATGCtt                                                                                        >  2:313297/1‑141 (MQ=17)
                                                    ttGGGGCATATTCTTGAACCCAACGGTATACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGt                                                                                      <  1:447623/141‑1 (MQ=17)
                                                    tgggggCATATTCTTGAACCCACCCGTAGACCGTTGAATGATGAACCTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGt                                                                                      <  1:162038/139‑1 (MQ=12)
                                                       gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTaa                                                                                   >  1:320265/1‑141 (MQ=32)
                                                       gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATTTCAGATATATCACGATAACTCAATTTAGATCTTAGATAGTAGCCAACGGCTACAGTG‑CTACGATCCTTGTTaa                                                                                   >  1:22184/1‑141 (MQ=11)
                                                        ggCCTATTCTTGAACCCAACGGTAGCCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATACCACGATAACTCAATGTATGTCTTAGATAGTAGCCACCGGCCACAGTG‑ATAACCTCCCTGTTaca                                                                                  >  2:92202/1‑139 (MQ=11)
                                                        agCATATTCTTGAACCCAACGGTAGACAGTTGAATGATGAACGTTTACACCACGATCCCTTCATATTTCAGAAATATCACGATAACTCACTGTAAATCTTAGATCATGGACAACGCCT‑CAATGAATAACATCCTTGTTcaa                                                                                  >  2:92214/2‑141 (MQ=11)
                                                                cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTCTATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGGTtat                                                                          >  2:2533/1‑141 (MQ=17)
                                                                    aaCCCATCCGTAGATTGTTTAATGATGAACGTTTACACGACGTCCTCCTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCt                                                                      <  2:537633/141‑1 (MQ=11)
                                                                      cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGa                                                                    >  2:394832/1‑141 (MQ=32)
                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                          >  1:303520/1‑141 (MQ=32)
                                                                                     ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                     >  2:773191/1‑141 (MQ=32)
                                                                                                aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttt                                          >  2:581495/1‑141 (MQ=17)
                                                                                                    ttACACCACGTTCCTTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                      <  2:610387/141‑1 (MQ=17)
                                                                                                    ttACACCACCTTCCCTTAACATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                      <  1:28769/141‑1 (MQ=17)
                                                                                                    tgacacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                      <  1:14448/139‑1 (MQ=17)
                                                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                    <  1:124833/141‑1 (MQ=32)
                                                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                    <  1:129676/141‑1 (MQ=32)
                                                                                                         ccacGTTCCCTTAATATTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                 <  1:23867/141‑1 (MQ=17)
                                                                                                            cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                              >  2:18857/1‑141 (MQ=32)
                                                                                                               tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                           >  2:497075/1‑141 (MQ=32)
                                                                                                               tCCCTTAACATATCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                           <  2:507431/141‑1 (MQ=17)
                                                                                                                     aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                     <  1:632673/141‑1 (MQ=35)
                                                                                                                           tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt               >  1:440501/1‑141 (MQ=35)
                                                                                                                            cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt              >  2:773679/1‑141 (MQ=35)
                                                                                                                              aatatatCACGATAACTCAATGGATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg            <  2:144553/140‑1 (MQ=17)
                                                                                                                                tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAATTTATAATATAAATTTAACTTTGCa          >  2:471689/1‑141 (MQ=255)
                                                                                                                                tatatCACGATAACTCAACGGATATCGTAGATAGAAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATGCAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCAACTTTGCa          <  1:202110/141‑1 (MQ=11)
                                                                                                                                taaatCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa          <  1:353346/138‑1 (MQ=21)
                                                                                                                                  tatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc        <  1:416965/141‑1 (MQ=255)
                                                                                                                                        gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAccc  >  2:193791/1‑141 (MQ=25)
                                                                                                                                            |                                                                                                                                         
CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACC  >  CP000730/36331‑36607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑