Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,174 T→C 28.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,1740TC28.7% 545.6 / inf 568K126K (AAA→AAGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (205/193);  new base C (72/90);  total (284/284)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 582 total aligned reads displayed.

ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑36314
                                                                                                                                            |                                                                                                                                            
  atGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTAGAATGGCAGTCAGGTGTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGATGCCTGATCTGTAATTACCTTTTGAGGatt                                                                                                                                            >  2:489526/1‑141 (MQ=11)
   tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                                           >  1:489824/1‑141 (MQ=255)
   tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTACCTACCTTTTGCGGCTTa                                                                                                                                           >  2:294596/1‑141 (MQ=11)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTGAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                                          <  1:295541/141‑1 (MQ=11)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAATTTTAAGAGCTTTAATTACCTTAATCATTGCTATGTTCGTTGTAGGTGCCTGTTCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  2:349657/1‑141 (MQ=11)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  1:737126/141‑1 (MQ=255)
     aCGGTGATCTTACTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTCATTACCTTTTGAGGTTTTcc                                                                                                                                         >  2:754082/1‑141 (MQ=11)
      cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCa                                                                                                                                        <  1:265298/141‑1 (MQ=18)
        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                                      <  1:483391/141‑1 (MQ=12)
        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      >  2:500566/1‑141 (MQ=255)
         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     >  1:507094/1‑141 (MQ=12)
           ttCTTGCTCAATGAGCTTATTCAGATCTTTCGATCTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACGTTTTGAGGTTTACCAAATTg                                                                                                                                   <  1:465696/140‑1 (MQ=11)
            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTGCTAACTTGt                                                                                                                                  >  2:456340/1‑141 (MQ=11)
            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  2:662342/1‑141 (MQ=12)
            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTGAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  2:125594/1‑141 (MQ=12)
             cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  1:533228/1‑141 (MQ=12)
             cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                 >  1:343142/1‑141 (MQ=11)
              ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTCCCAAATTGttt                                                                                                                                >  1:143435/1‑141 (MQ=255)
                gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              >  1:383206/1‑141 (MQ=255)
                gCTCAATGAGGTTATTCAGATATTGCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              <  1:194525/141‑1 (MQ=11)
                 ctctATGAGTTTATTCAGATATTTCGATGTACAATGACAGGCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAAt                                                                                                                             <  2:305487/141‑1 (MQ=11)
                       gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGTTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACg                                                                                                                       >  2:330529/1‑141 (MQ=255)
                       gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACg                                                                                                                       >  2:145892/1‑141 (MQ=12)
                       gAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACg                                                                                                                       <  1:100419/141‑1 (MQ=255)
                       gAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGCCg                                                                                                                       >  2:679871/1‑141 (MQ=17)
                        aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  2:618529/141‑1 (MQ=17)
                         ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                     >  2:58691/1‑141 (MQ=12)
                         ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                                     <  1:622166/141‑1 (MQ=11)
                          ggtATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  1:107158/139‑1 (MQ=12)
                          gTTATTCAGATAATTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAGCTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGttt                                                                                                                    >  1:542390/1‑141 (MQ=11)
                           ttattaAGCTATTACGATGTACCATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATCACCTTAACCAGTGCTACCGTCGTTTAAGGCGCCTGATCTGTAGTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  2:241539/1‑141 (MQ=11)
                           ttATTCAGATACTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  1:35741/141‑1 (MQ=12)
                             aTTCAGATATTTCGATGTACAATGACAGTCCGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGGTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGCCGTTTCAt                                                                                                                 >  1:309916/1‑141 (MQ=255)
                               tCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTCTGATaa                                                                                                               >  2:562478/1‑141 (MQ=255)
                                cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              <  2:601944/141‑1 (MQ=11)
                                  gATATTTCGATTTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:408400/141‑1 (MQ=12)
                                  gATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:252966/141‑1 (MQ=255)
                                    tatTTCGATGGACAATGGCAGTTAGGTTTAAGTTTAAATGCTTTAATTACTTTAACCATCGCTACCTTCGTTGACGGCGCCTGAGCTGTAATTAACTGTTGAGGTTTAACAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  1:31870/141‑1 (MQ=11)
                                     ttttCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         <  2:207482/140‑1 (MQ=17)
                                       ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                                       <  1:659305/141‑1 (MQ=17)
                                       ttCGATGTACAATGACAGTCAGGTTTAAGTTTACAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                                       <  1:706512/141‑1 (MQ=11)
                                        tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      >  1:498232/1‑141 (MQ=17)
                                         tgATGTACAATGGCAGTCAGGTTTTAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCCGTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  1:421631/140‑1 (MQ=11)
                                         cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  1:279422/141‑1 (MQ=11)
                                          gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAACAATTGCTGCCTTCGTTGAAGTTGCCTGATCTGTAACTAACTTTTTAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                                    >  1:453784/1‑141 (MQ=11)
                                          gATGTACAATGACAGTCAGGTTTAAGGTTAAAAGCTTTAATTATTTTAGACATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCAAATgat                                                                                                    >  2:155933/1‑139 (MQ=11)
                                               aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCTTGATCTGTCATTACCTTGTGAAAGTTCCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAta                                                                                               >  1:333652/1‑140 (MQ=11)
                                               aCAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGTCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                               <  2:362764/141‑1 (MQ=11)
                                                 aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTATCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             <  2:548663/141‑1 (MQ=17)
                                                                                       ttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       <  1:722230/141‑1 (MQ=9)
                                                                                       ttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       <  2:519956/141‑1 (MQ=14)
                                                                                          accATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                                    <  2:77726/140‑1 (MQ=17)
                                                                                             aTTGCTACCTTCGGTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                                 <  1:232621/141‑1 (MQ=17)
                                                                                               tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                               >  1:616611/1‑140 (MQ=17)
                                                                                                gCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGTATGTCCCTCTGCa                                              <  2:571806/141‑1 (MQ=17)
                                                                                                 cTACCTTCGTTGAAGGTGCCTGATGTGTGATTAGCTTTTGGTGTTTCGGACTGTCTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             <  1:554224/141‑4 (MQ=11)
                                                                                                  tgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                                            >  2:632124/3‑141 (MQ=255)
                                                                                                  tgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                                            >  1:199088/3‑141 (MQ=255)
                                                                                                  tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            >  1:161876/1‑141 (MQ=32)
                                                                                                  tACCTTCGGTGAAGGTGTCTTATCTGTAATTACCTTTTGATGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            <  2:220143/141‑1 (MQ=255)
                                                                                                   gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCCAATTGTTTAATGAGTCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  1:112696/140‑1 (MQ=11)
                                                                                                   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCCTCa                                           >  1:124380/1‑141 (MQ=17)
                                                                                                   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                           <  1:562842/141‑1 (MQ=32)
                                                                                                    ccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCaa                                          <  1:81557/141‑1 (MQ=32)
                                                                                                    ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          >  2:320179/1‑141 (MQ=12)
                                                                                                    ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTGTGAGGCTTACTAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:167046/141‑1 (MQ=11)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                         <  2:479400/141‑1 (MQ=32)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                         >  1:692492/1‑141 (MQ=32)
                                                                                                      ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  2:13876/141‑1 (MQ=17)
                                                                                                         aacgAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAGTGGcc                                     >  1:246349/4‑140 (MQ=17)
                                                                                                              aGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatt                                >  1:145924/1‑140 (MQ=17)
                                                                                                                   ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTTGTTGTCTATTCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAAt                           >  1:594307/1‑141 (MQ=11)
                                                                                                                    cTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATa                          <  2:235813/141‑1 (MQ=32)
                                                                                                                      gacatGTAACTACCGTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTTATAAACGCATAGGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  2:714450/137‑1 (MQ=11)
                                                                                                                      gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATAACTAATCTATGTCCCTCCGCATCCATGGCACGATATaaaaacc                        >  2:113447/1‑137 (MQ=11)
                                                                                                                      gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTAGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGc                        >  1:518259/1‑141 (MQ=17)
                                                                                                                      gATCTGTAACTACCTTTTGAGGCTTACCAAAGTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  2:464104/141‑1 (MQ=255)
                                                                                                                       aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAATTTGCt                       >  1:584272/1‑141 (MQ=17)
                                                                                                                        tCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc                      >  1:606290/1‑141 (MQ=17)
                                                                                                                         ttGTAATGACCTTTGGAGGCTTACCAAATTGTTTAATGAGACTTTTGATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                     <  2:537073/140‑1 (MQ=11)
                                                                                                                         cTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGTCTGTCCCTATGTATCAATTTCACGATATAAATAGCTcc                     <  2:69540/141‑1 (MQ=11)
                                                                                                                          tGTCATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTCGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATCTCTAATGTATGTCCCTCTGCATCACTGGCACGATATAAATAGCTCCa                    <  2:183821/141‑1 (MQ=255)
                                                                                                                          tGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                    <  2:383297/141‑1 (MQ=12)
                                                                                                                          tGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                    >  1:671121/1‑141 (MQ=255)
                                                                                                                           gTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGACTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCCCGATCTCAATAGGTCCAt                   >  1:470087/1‑141 (MQ=11)
                                                                                                                            tAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATGTAAATAGCTCCAtt                  >  1:440356/1‑141 (MQ=255)
                                                                                                                             aagtACCTTCTGAGGTTTACCCAATTGCTTAATGCGACGTTCGATTAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAttt                 <  1:18060/138‑1 (MQ=11)
                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa                    <  2:392234/137‑1 (MQ=33)
                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa                    >  1:392234/1‑137 (MQ=33)
                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCCAATGTATGCCTCTCTGCATCAATGGCACGATACAAATAGCTCCAattc                >  2:8606/1‑140 (MQ=11)
                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             >  1:98132/1‑141 (MQ=17)
                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCACTTTCCt             <  2:424005/141‑1 (MQ=17)
                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCCTATGCTGATTGATTAGCGCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:272693/141‑1 (MQ=11)
                                                                                                                                  ccTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            >  2:433502/1‑141 (MQ=17)
                                                                                                                                   cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  2:651477/141‑1 (MQ=17)
                                                                                                                                    ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCtttt          >  2:726876/1‑141 (MQ=32)
                                                                                                                                    ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCtttt          >  1:710516/1‑141 (MQ=32)
                                                                                                                                    ttttGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtttt          <  1:370238/141‑1 (MQ=17)
                                                                                                                                         gggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtttt     <  1:215518/140‑1 (MQ=17)
                                                                                                                                          ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCAATTTTCCTTTTAtttag    >  2:445450/1‑139 (MQ=12)
                                                                                                                                            |                                                                                                                                            
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: