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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,1740TC30.3% 850.6 / inf 942K126K (AAA→AAGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (347/300);  new base C (114/170);  total (466/476)
Fisher's exact test for biased strand distribution p-value = 1.60e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-02

Only 100 of 967 total aligned reads displayed.

ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑36314
                                                                                                                                            |                                                                                                                                            
atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGt                                                                                                                                              <  1:467997/141‑1 (MQ=25)
   tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTa                                                                                                                                           >  2:1100601/1‑141 (MQ=255)
   tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                                           <  2:115247/141‑1 (MQ=35)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGTAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTATCTTTTGAGGTTTAc                                                                                                                                          <  2:756845/141‑1 (MQ=255)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                                          <  2:294325/141‑1 (MQ=21)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                                          >  2:134110/1‑141 (MQ=255)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTACGTTTAAAAGATTTAATTACTGTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGCGGTTTAc                                                                                                                                          >  2:946279/1‑141 (MQ=255)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                                          >  1:491226/1‑141 (MQ=255)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                                          <  1:47853/141‑1 (MQ=11)
       ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaa                                                                                                                                       <  1:600406/141‑1 (MQ=11)
       ggTGATCTTGCTCAATGAGGTTATACAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATGACCTTTTGAGGTTTACCaa                                                                                                                                       >  1:705207/1‑141 (MQ=18)
        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      >  2:1088096/1‑141 (MQ=11)
        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      <  2:705911/141‑1 (MQ=11)
         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAt                                                                                                                                     >  1:43003/1‑141 (MQ=11)
           aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTTCCTTCGTTGTAGGTGCCTGATCTGTAATTAGCTTTTGTGGTTTACCAAATTg                                                                                                                                   >  1:53478/1‑141 (MQ=11)
           aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                                                                   >  1:762306/1‑141 (MQ=17)
           aTCTTGCTCAATGAGGTTATTAGCATTTTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                                                                   >  1:1098461/1‑141 (MQ=11)
              ttGCTCAATGAGGTTTTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                                                                <  1:999476/141‑1 (MQ=12)
              ttGCTCAATGAGGTTATTCAGATATTCCGATGTACAATGGCAGTCAGGTGTAAGTTTAAAAGCTTTAATTACATTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                                                                >  2:683325/1‑141 (MQ=11)
                 cTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCATATTGTTTAAt                                                                                                                             >  1:828793/1‑141 (MQ=255)
                 cTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                             >  1:339823/1‑141 (MQ=17)
                 cTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAAt                                                                                                                             <  1:936981/141‑1 (MQ=11)
                   gacgacgcattatCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAATGATAAGAAGACTGTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           <  1:804590/129‑1 (MQ=11)
                   cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           <  2:7915/141‑1 (MQ=17)
                    aaaggggTTAGTGAGATATTTCGATGTAAAATGCCAGAAAGGTTTAGGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAGGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                                                                          <  2:229110/136‑1 (MQ=11)
                    aaTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                                                                          >  2:1056576/1‑141 (MQ=12)
                     aTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATgaga                                                                                                                         >  1:626796/1‑141 (MQ=11)
                      tGAGGTTATTCAGATCTTTCGATGTACAATGGCAGCCAGGTCTACGTTTAAAACCTTTAATTACTTTAACCATTGCTACCTTCGTTGACGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                                                                                                                        >  2:96575/1‑141 (MQ=11)
                      tGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGAc                                                                                                                        >  1:1054284/1‑141 (MQ=255)
                      ggAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTGTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                                                                                                                        <  2:416490/140‑1 (MQ=12)
                        gggTTATTCAGATATTTCGGTTTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTCAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  2:275033/140‑1 (MQ=11)
                        gggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  2:190872/140‑1 (MQ=17)
                        aGGTTATTCAGATATTTCGATGTACAAGGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  2:211290/141‑1 (MQ=12)
                         ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTCACCATTGCTACCTTAGTTGAAGGTGCCTGATCTGTCCTTCCCTTTCGAGGTATACCAGCTTGTGTAATTAGACGtt                                                                                                                     >  2:577419/1‑141 (MQ=11)
                         ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                     >  1:954352/1‑141 (MQ=17)
                          gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  2:165090/141‑1 (MQ=17)
                          gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  1:982372/141‑1 (MQ=17)
                           ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCCTTGAAGGTGCCTGATCTTTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtttt                                                                                                                   >  2:472611/1‑140 (MQ=255)
                           ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  2:256838/141‑1 (MQ=255)
                           ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  2:1076519/141‑1 (MQ=11)
                            tATTCAGATATTTCGATGTACAATGGCAGTCAGATTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  2:617482/141‑1 (MQ=12)
                             tttCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                                 <  2:393608/140‑1 (MQ=17)
                                cAGATATTTCGATATACAATGGCAGCCAGGTGTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCAGGAGATATAAGGACATTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              >  2:335259/1‑141 (MQ=11)
                                  gATATTTCGATGTACAATGGCAGTCAGGTTTACGTTTATAAGATTTAATTACTTTAATCATTGCTACCTTCTTTTAAGGTGCCTTATCTTTAATTACCTTTTGTGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:369730/141‑1 (MQ=11)
                                  gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            >  2:880356/1‑141 (MQ=17)
                                  gATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  1:879001/141‑1 (MQ=11)
                                     atTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGAAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         <  2:596412/141‑1 (MQ=11)
                                     agttCGATGTACAACGGCAGGTAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         <  2:928736/139‑1 (MQ=11)
                                      tttCGATGTACAATGGCAGTCAGGTTTAGGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCata                                                                                                        >  1:445700/1‑141 (MQ=12)
                                        tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      <  1:215276/141‑1 (MQ=11)
                                         cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     >  2:637085/1‑141 (MQ=12)
                                              tACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAt                                                                                                >  2:82613/1‑141 (MQ=12)
                                              tACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGCGACGTTTGATAAGCGCATATGCTGAAt                                                                                                >  2:659213/1‑141 (MQ=11)
                                               aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGTTGAATg                                                                                               >  1:196793/1‑141 (MQ=12)
                                                 aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             <  1:14120/141‑1 (MQ=11)
                                                 aaTGACAGTTAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             <  2:15539/141‑1 (MQ=11)
                                                   ggACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           <  2:738763/140‑1 (MQ=255)
                                                    gACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTAt                                                                                          >  2:777545/1‑141 (MQ=11)
                                                                    ttaaaaaGCTTTAATTACCTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCaa                                                                          <  1:376450/138‑1 (MQ=0)
                                                                     ttcaaaGCTTTATTTACTGTAGCCATTGCTACCTTCGTTGATGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAAc                                                                         <  2:495216/138‑1 (MQ=0)
                                                                                      tttaccgagtGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACTTTTCGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        <  1:15352/132‑1 (MQ=11)
                                                                                         aGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTTCCAAATTGTTTAATGAGTCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  2:436843/141‑1 (MQ=17)
                                                                                              ttGCTACCTTCGTTGAAGGTGCCTGGTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                                >  1:554651/1‑141 (MQ=17)
                                                                                              ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                                <  2:2455/141‑1 (MQ=32)
                                                                                              ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                                <  2:984068/141‑1 (MQ=32)
                                                                                               tgctgcCGTCGTTGAAGGTGCTTGATCTGTTACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                               <  1:1027077/141‑2 (MQ=11)
                                                                                               tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               >  2:323035/1‑141 (MQ=32)
                                                                                                gCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgta                                              <  1:1073946/141‑3 (MQ=255)
                                                                                                 cTACCTTCGTTGAAGGTGCTTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                             <  1:247985/141‑1 (MQ=17)
                                                                                                  tACCTTCGTTGAAGGTGCCTGATCTGTAATTTCCTTTTGAGTTTTACCAAATTGTTTAATGAGACGTTTGATAACCGCATACGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            <  1:902651/141‑1 (MQ=255)
                                                                                                  tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            >  1:401856/1‑141 (MQ=32)
                                                                                                   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGcttta                                           >  1:682326/1‑137 (MQ=17)
                                                                                                    ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:552386/141‑1 (MQ=255)
                                                                                                      ttagtcataaatGCCTCATTTGTTATTATCCGTGGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  1:1081135/130‑1 (MQ=11)
                                                                                                      ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  1:173107/141‑1 (MQ=17)
                                                                                                      ttCGTTGAAGGGCCCTGATCGCTAACTACTTTTTGAGGCTTACCATATTGTTTGCTTGTACCTCTGATACACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  1:31794/141‑1 (MQ=11)
                                                                                                       tCGTTGAAGGTGCCTTATCTGTAATTAACTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATgg                                       <  1:650737/141‑1 (MQ=17)
                                                                                                        caacgAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                                      >  2:355279/5‑141 (MQ=17)
                                                                                                         gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATTCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATTGCa                                     >  1:232236/1‑141 (MQ=12)
                                                                                                          acgAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    >  1:844571/3‑141 (MQ=17)
                                                                                                           tGAAGGTGCCTGATTTGTAATTACCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACg                                   <  1:1002839/141‑1 (MQ=17)
                                                                                                             aaGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTTTCAATGGCACGat                                 >  1:725309/1‑141 (MQ=17)
                                                                                                               ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGTTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatat                               >  1:293196/1‑141 (MQ=17)
                                                                                                                gTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAtata                              <  2:844465/141‑1 (MQ=12)
                                                                                                                 tGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATaa                             >  2:144128/1‑141 (MQ=32)
                                                                                                                    cttaTCTGTAATTACCTTTTGATGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATa                          <  1:124744/138‑1 (MQ=17)
                                                                                                                    cTGATCTGTAACTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACTCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATa                          <  1:779473/141‑1 (MQ=17)
                                                                                                                      gATGTGTCATTATCCTGTGAGGTTTACCAAATTGTTTAATGAGCCGTTTGATAAACGCAGATGCTGAATGATTATCTCGTTACTTACGCCACCACATATCTAATGTATGTCCCTCTGCATCAATGGCGCGGTATCAATAGc                        >  2:511907/1‑141 (MQ=11)
                                                                                                                      gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  1:962597/141‑1 (MQ=17)
                                                                                                                      gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  2:1027207/141‑1 (MQ=12)
                                                                                                                       gtctgaaaccaCCTTTTGAGGCGTACCAAATTGTTTAATGGGAAGTTTGATAAACGCATATGCTGAATGATTATCGCGTTGCTTACGCAACTAAATATTTAATGTATCTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  2:727123/131‑1 (MQ=1)
                                                                                                                       aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACAAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt                       <  2:446534/141‑1 (MQ=17)
                                                                                                                         cTGTAATCACCTTTTGAGTTTTACCCAATTCTTACCTGACATGTTCGTTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                     <  1:748728/141‑1 (MQ=11)
                                                                                                                           gTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAGTATCTAATGTATGTCCCTCTGCATCCATGGCAAGATATACATAGCTCCAt                   <  2:813637/141‑1 (MQ=255)
                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCCACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            >  2:531716/1‑141 (MQ=17)
                                                                                                                                   cTTTTGATGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  1:515744/141‑1 (MQ=12)
                                                                                                                                   cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  2:911369/141‑1 (MQ=32)
                                                                                                                                        gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACTCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttt      >  1:579703/1‑141 (MQ=17)
                                                                                                                                            tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAACTATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGag  >  1:917603/1‑140 (MQ=17)
                                                                                                                                            tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAATTTGAt  >  1:866754/1‑141 (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: