Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,174 T→C 26.7% K126K (AAA→AAG USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,1740TC26.7% 574.0 / inf 524K126K (AAA→AAGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (179/201);  new base C (58/83);  total (240/284)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 539 total aligned reads displayed.

ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑36314
                                                                                                                                            |                                                                                                                                            
  atGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGttt                                                                                                                                            <  2:48266/141‑1 (MQ=255)
  aggACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTTCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGttt                                                                                                                                            <  1:32059/139‑1 (MQ=255)
   tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                                           >  1:274197/1‑141 (MQ=255)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTGAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                                          <  1:152381/141‑1 (MQ=255)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACCATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGTTCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                                         >  1:42054/1‑141 (MQ=11)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  2:136889/1‑141 (MQ=255)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  1:394436/1‑141 (MQ=255)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  2:393895/141‑1 (MQ=19)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  1:54609/1‑141 (MQ=21)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGAACTGTAATTACCTTTTGCGGTTTAcc                                                                                                                                         >  1:81557/1‑141 (MQ=255)
       ggTGATCTTTCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaa                                                                                                                                       <  1:75483/141‑1 (MQ=12)
        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGCCTTACGaaa                                                                                                                                      >  2:33712/1‑141 (MQ=11)
         tGATCTTGCTCAATTACGGTATTCAGATATTTCGATGTACAATGGCAGGCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTATCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     <  1:179589/141‑1 (MQ=11)
         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTTCCAAAt                                                                                                                                     >  2:247469/1‑141 (MQ=12)
         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     >  1:383931/1‑141 (MQ=12)
         tGATCTTGCTCAATGAGGTTAGTCATATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     >  2:281298/1‑141 (MQ=255)
          gATCTTGCTCAATGAGGTTATTCACATATTTCTATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAtt                                                                                                                                    <  2:117799/141‑1 (MQ=11)
           agcTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTg                                                                                                                                   >  2:200002/3‑141 (MQ=11)
            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  <  2:54068/141‑1 (MQ=17)
            tCTTGCCCAATGGGGTTATCCCGATATTTCGATGGACAATGGCAGTCAGGTTTAAGTTGAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  <  1:144792/141‑1 (MQ=11)
              ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                                                                >  1:101446/1‑141 (MQ=12)
              ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                                                                >  2:209934/1‑141 (MQ=17)
              ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGttt                                                                                                                                <  2:162610/141‑1 (MQ=11)
                gCTCAATGAGGTTATTCAGATCTTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              <  1:205968/141‑1 (MQ=11)
                gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              >  2:40599/1‑141 (MQ=255)
                 cTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAAt                                                                                                                             <  1:362424/141‑1 (MQ=255)
                 cTCAAGGAGGTTATTCAAATATTTCGATGGACAATGGCAGTCAGGTTTAATTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                             <  2:175987/141‑1 (MQ=11)
                   cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           <  2:61063/141‑1 (MQ=17)
                      tGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                                                                                                                        >  1:307188/1‑141 (MQ=17)
                      tGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                                                                                                                        >  1:115304/1‑141 (MQ=17)
                        aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTCCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      >  1:89725/1‑141 (MQ=12)
                         ggTTATTCAGATATTTCGATGTACGATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCTTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                                     >  1:288224/1‑141 (MQ=11)
                         ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                     >  1:62373/1‑141 (MQ=12)
                          gTTATTCATATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  2:362747/141‑1 (MQ=12)
                           ttATCCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  1:130043/141‑1 (MQ=12)
                           ttATACAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTCTAGCCATTGCTGCCTTCGTTGAAGGTGCCTTATCTGTAACTGCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtttt                                                                                                                   <  2:226788/141‑2 (MQ=11)
                               tCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATaa                                                                                                               <  1:350118/141‑1 (MQ=11)
                                 aGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                                                                             >  1:315601/1‑141 (MQ=12)
                                   atatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                                           >  1:324245/1‑141 (MQ=12)
                                   atatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAACTTGCTTCATGCGACGTTTGCTCAACGc                                                                                                           >  2:370418/1‑141 (MQ=11)
                                    tatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTATCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  2:286412/141‑1 (MQ=12)
                                    tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAACTTGTTTAATGAGACGTTTGATGAACGCa                                                                                                          >  1:157947/1‑141 (MQ=11)
                                    tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  1:35496/141‑1 (MQ=255)
                                        tCGATGTACCATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      <  1:238377/141‑1 (MQ=12)
                                         cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  2:68341/141‑1 (MQ=255)
                                          gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCaacgaaa                                                                                                    >  2:18419/1‑135 (MQ=255)
                                          gATGTACAATGAGAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCGTACCAAATTGTTTAATGAGAGGTTTGATAGACGCATATGCt                                                                                                    >  2:298928/1‑141 (MQ=11)
                                          gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                                    <  2:141450/141‑1 (MQ=11)
                                           aTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGTTAAACGCATATGCTg                                                                                                   >  1:370367/1‑141 (MQ=12)
                                            tGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                                  >  1:298276/1‑141 (MQ=17)
                                               aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                               <  1:281298/141‑1 (MQ=17)
                                                   tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           <  1:389990/141‑1 (MQ=255)
                                                                                       ttAGCCATTGCTGCCTTCGTTGAAGGTGCTTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       <  2:387323/141‑1 (MQ=255)
                                                                                       ttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       <  1:230079/141‑1 (MQ=9)
                                                                                       ttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTGACGCAACCAAATATCTAatgtatgt                                                       <  2:381612/141‑1 (MQ=14)
                                                                                        taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGGTAAACGCACAGGCTGAATGATTATCTCGTGGCTTATGCAGCCAAATATATAATGTGTGTc                                                      >  2:219489/4‑141 (MQ=11)
                                                                                          accATTGCTACCTTCGTTGACGGTGCCTCATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAGAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                                    >  1:256264/2‑141 (MQ=12)
                                                                                           ccATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATACCTAATGTATGTCCct                                                   <  2:322680/141‑1 (MQ=255)
                                                                                             tttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                                 <  2:263434/140‑1 (MQ=17)
                                                                                             aTTGCGTCCGTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                                 <  1:340769/141‑1 (MQ=12)
                                                                                               tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTCACCAAATTGTTTAATGAGACGGTTGATAAACGCATATGCTGACTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                               >  2:258032/1‑140 (MQ=12)
                                                                                                gctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgta                                              >  1:242043/5‑139 (MQ=255)
                                                                                                gCTATCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                              <  2:109520/141‑1 (MQ=17)
                                                                                                 ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTAt                                             >  1:251027/4‑141 (MQ=255)
                                                                                                  tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            <  2:73470/141‑1 (MQ=32)
                                                                                                    ccTTCGTTGCAGGTGACTTATCTGCAATTACCTTTTGAAGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATATCGTTGCTTACGCAACCAAATAACTCATGTATGTCCCGCTGCATTaa                                          >  2:124268/1‑141 (MQ=11)
                                                                                                     ctacGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                         <  1:21375/138‑1 (MQ=17)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAAGCATATCTTGAAGGAGTATCTCGCTGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                         <  2:83189/141‑1 (MQ=11)
                                                                                                       tCGTTGAAGGTGCCTGATCCTGAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATgg                                       <  1:302056/141‑1 (MQ=12)
                                                                                                        cGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                                      <  2:16523/141‑1 (MQ=32)
                                                                                                          ttGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:110900/141‑1 (MQ=32)
                                                                                                          tggAAGGTGCCTGCTCTGTTATTACCTTTTTAGGTTTACTAAATTTTTTAATCAGTCGTTTGATAAACGCATATGCTGAATTATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:140894/139‑1 (MQ=11)
                                                                                                            gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGa                                  <  1:98125/141‑1 (MQ=17)
                                                                                                            gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                                  >  1:6710/1‑141 (MQ=32)
                                                                                                             aaGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGGGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGat                                 <  1:164223/141‑1 (MQ=17)
                                                                                                              aGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCGAATATTTAATGTATGTCCCTCTGTATCAATGGCACGatt                                >  1:219090/1‑140 (MQ=11)
                                                                                                                  gCCTGATATGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCGACCAAATATCTAATGTGTGTCCCTCTGCAGCAATGGCACGATATaaa                            >  1:137222/1‑141 (MQ=255)
                                                                                                                     tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAg                         >  2:311075/1‑141 (MQ=32)
                                                                                                                      gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  1:200076/141‑1 (MQ=12)
                                                                                                                       aTCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  2:268732/141‑1 (MQ=12)
                                                                                                                        tCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc                      <  1:333511/141‑1 (MQ=12)
                                                                                                                         cTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTTTCTCTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTcc                     >  1:183751/1‑141 (MQ=11)
                                                                                                                            taactACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtt                  >  2:347861/5‑141 (MQ=12)
                                                                                                                              actACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                                      >  1:331627/3‑103 (MQ=12)
                                                                                                                              actACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                                      <  2:331627/101‑1 (MQ=12)
                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAtttt                <  2:192694/141‑1 (MQ=32)
                                                                                                                               ttACCTTTTGAGGTTTACCAAATTGTGTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc               >  2:237573/1‑141 (MQ=17)
                                                                                                                                tgccTTTTGGCGTTGACCAAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              <  2:324785/139‑1 (MQ=11)
                                                                                                                                tACCTTTTGAGGTTTACCAAATTTTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTcc              <  2:157061/141‑1 (MQ=17)
                                                                                                                                tACCTTTTGAGGCTTACCCAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              <  1:123537/141‑1 (MQ=12)
                                                                                                                                tACCGTTCGAGGCTTACCAAATTGTTAAATTATACGTTTGATACACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              <  1:19215/141‑1 (MQ=11)
                                                                                                                                 aCCGTTTCAGGCTTACCAAATTGTTTAATGAGACGTGTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  2:275071/141‑1 (MQ=11)
                                                                                                                                  cgttttGAGGTGTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:90071/139‑1 (MQ=17)
                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:148613/141‑1 (MQ=17)
                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCtt            >  1:376735/1‑141 (MQ=32)
                                                                                                                                   cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACCCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTGCttt           <  1:268488/141‑1 (MQ=12)
                                                                                                                                    ttttGCGGTTTACCAAGTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAAGATCTAATGTATGTCCCTCTGCATCAAGGGCACGATATAAATAGCTCCATTTTCCtttt          <  2:296561/141‑1 (MQ=255)
                                                                                                                                        gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttt      >  1:72140/1‑141 (MQ=17)
                                                                                                                                         aGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtttt     <  2:354204/141‑1 (MQ=17)
                                                                                                                                            tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt  <  2:376015/141‑1 (MQ=17)
                                                                                                                                            |                                                                                                                                            
ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑36314

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: