Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,174 T→C 30.2% 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.2% 613.2 / inf 693K126K (AAA→AAGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (236/242);  new base C (86/123);  total (327/366)
Fisher's exact test for biased strand distribution p-value = 5.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 711 total aligned reads displayed.

ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑36314
                                                                                                                                            |                                                                                                                                            
atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTGTAGGt                                                                                                                                              >  2:663500/1‑141 (MQ=255)
 tatGACGGTGATCTTGCTCAATGAGGTTATTCATATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGtt                                                                                                                                             <  2:553407/141‑1 (MQ=255)
    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                                          <  1:485102/141‑1 (MQ=25)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  1:46860/141‑1 (MQ=25)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                                         <  2:682733/141‑1 (MQ=11)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                                         >  1:268577/1‑141 (MQ=11)
     aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                                         <  1:137031/141‑1 (MQ=11)
      cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCa                                                                                                                                        <  2:329505/141‑1 (MQ=21)
        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                                      <  1:805402/141‑1 (MQ=17)
         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     <  1:211558/141‑1 (MQ=12)
         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     <  2:648897/141‑1 (MQ=17)
         tGATCTTGCTCAATGAGGTTATTCAGATATTGCGATGTACAATGGCCGTCAGTTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                     <  2:593443/141‑1 (MQ=11)
          gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACGAAAtt                                                                                                                                    >  2:571700/1‑141 (MQ=12)
          gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAtt                                                                                                                                    <  2:839214/141‑1 (MQ=12)
            tCTTGCTCAATGAGGTTCTTCAGATATTTCGATGTACAATGACAGTCAGCTTTAAGTTTAACAACTTTCATTATTGTAGCCTTTGCTTCCTTCGTTGAAGGTGCCTGATCTGTCACTACATTTTGAGGCTTACCAAATTGt                                                                                                                                  >  2:100644/1‑141 (MQ=11)
            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  1:887413/1‑141 (MQ=12)
            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  2:81839/1‑141 (MQ=12)
             cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                 >  2:526952/1‑141 (MQ=11)
             cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTACGCAAATTGtt                                                                                                                                 >  2:36061/1‑141 (MQ=11)
               tGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTTCCAAATTGtttt                                                                                                                               >  2:182482/1‑140 (MQ=12)
                gcccAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              <  2:470693/138‑1 (MQ=11)
                gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGTCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                              <  2:212322/141‑1 (MQ=12)
                gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                              >  1:91127/1‑141 (MQ=12)
                gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              <  1:305700/141‑1 (MQ=11)
                  tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                                            >  1:869672/1‑141 (MQ=17)
                       gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTATTGAGACg                                                                                                                       >  1:234659/1‑141 (MQ=12)
                       gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACg                                                                                                                       <  2:240580/141‑1 (MQ=17)
                        aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  1:443446/141‑1 (MQ=17)
                         ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTGACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                     >  1:572785/1‑141 (MQ=12)
                          gTTATTCAGATATTTCGATGTACAAGGGCAGTCAGGTTTCAGTTTAAAAGCTTTAATGACTTTAACCATTGCTACCGTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  2:788660/141‑1 (MQ=11)
                           ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  1:123277/1‑141 (MQ=255)
                             aTTCAGATATTGCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                                 <  2:761119/141‑1 (MQ=11)
                              ttCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                                >  2:299263/1‑141 (MQ=17)
                              ttCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                                >  1:836404/1‑141 (MQ=255)
                                cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              <  1:750000/141‑1 (MQ=17)
                                  gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATACACg                                                                                                            >  2:913621/1‑141 (MQ=255)
                                  gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  1:859113/141‑1 (MQ=17)
                                    tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          >  2:165441/1‑141 (MQ=11)
                                    tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  1:775671/141‑1 (MQ=11)
                                     atTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         >  2:673532/1‑141 (MQ=32)
                                      tttCGATGTACCATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGTGGATTACCAAATTGTTAAATTAGACGTTTGATAAACGCata                                                                                                        <  1:583817/141‑1 (MQ=11)
                                      tttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCata                                                                                                        <  2:603802/141‑1 (MQ=17)
                                      tgtCGAGGTAGAATGACAGTCAGGTTTACGTGTAAAAGCTTTAATGATTTTAGCCATTGCTGCCCTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCata                                                                                                        <  1:518513/139‑1 (MQ=11)
                                       ttCGATGTGCCATGACAGTCAGTTTTAAGTTTAAAATCTGTAATTATTTTACCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAGTGTTTAATGAGACGTTTGATAAACGCatat                                                                                                       <  2:190544/141‑1 (MQ=11)
                                        tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      >  2:502498/1‑141 (MQ=12)
                                          gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTACTGAGACGTTTGATAAACGCATATGCt                                                                                                    <  1:330966/141‑1 (MQ=17)
                                           aTGTACAATGGCAGTCAGGCTTAAGCTTAAAAGCTTTAATTACTTTAACCATAGCTACCTTCGTTGAAGGTGCCTGATCTGTAGTTACCTTGTGAAGTTTACCAAATTGTTTAATGAGACGTTTGATAAACACATAGGCTg                                                                                                   >  2:383531/1‑141 (MQ=11)
                                                cAATGTCAGTCAGGTTTAAGTTTAAAAGCTTTAATTCCTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                              <  1:194423/141‑1 (MQ=12)
                                                 aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             <  2:271388/141‑1 (MQ=11)
                                                    gACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTAt                                                                                          <  2:696357/141‑1 (MQ=11)
                                                                                   tattttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTCAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgt                                                           <  1:649944/138‑1 (MQ=1)
                                                                                      tttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        <  1:823722/141‑1 (MQ=11)
                                                                                       ttaaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTACGCAACCAAATATCTAatgtatgt                                                       >  2:492686/5‑141 (MQ=14)
                                                                                       ttaaccATTGCTACCCTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       >  1:677482/5‑141 (MQ=14)
                                                                                         aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     >  1:379448/3‑141 (MQ=17)
                                                                                          accATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                                    >  2:736987/2‑141 (MQ=17)
                                                                                          accATTGCTACCGTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATCTCATGTATGTcac                                                    >  2:440250/2‑139 (MQ=255)
                                                                                             aTTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATCAACGCATATGCTGAATGGTTATCTCGTTGCTTACGCAACCAACCATCTAATGTATGTCCctct                                                 >  1:539028/1‑141 (MQ=11)
                                                                                             aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                                 <  2:634596/141‑1 (MQ=32)
                                                                                              ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                                >  1:256305/1‑141 (MQ=32)
                                                                                               tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAAATATCTAATGTATGTCCCTCTgt                                               >  1:69509/1‑140 (MQ=17)
                                                                                                gatgccTTTGTGGAAGGTGCCTGACTTTTAACTATCTTTTGAGGCTTACTAAATTGTTTAATGTGACGTTTGATAAACGTATTTGCTGAATTATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgta                                              <  1:44878/137‑3 (MQ=11)
                                                                                                gCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCTACCATAGATCTAATGTATGTCCCTCTgta                                              <  2:416902/141‑3 (MQ=11)
                                                                                                 ctgccTTCGTTGAAGGTGCCTGTTCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             <  1:724618/138‑4 (MQ=11)
                                                                                                  tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTAta                                            >  2:801113/1‑140 (MQ=17)
                                                                                                    ccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAAGTTTGAAAAACGCATATGCTGAATGATTCTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          >  2:712048/1‑141 (MQ=255)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATTAGACGTGTGATAAACGCAGATGCTCAATGATTATCGCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         >  2:620/1‑141 (MQ=11)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         >  2:574783/1‑141 (MQ=17)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCagt                                         >  2:40838/1‑139 (MQ=17)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATACTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                         >  1:363141/1‑141 (MQ=17)
                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAACTACCTGTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         <  2:855879/141‑1 (MQ=255)
                                                                                                      ttCGTTGAAGGCGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  2:875408/141‑1 (MQ=255)
                                                                                                       tCGTTGAATGTGCCTGATTTGTAATTACCTGTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATgg                                       <  2:904747/141‑1 (MQ=12)
                                                                                                        cGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                                      >  2:342127/1‑141 (MQ=32)
                                                                                                            gAAGGTGCCTGATCTGTAATCACCTTTTGAGGTTTACCAAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACACAACCCAATATCTAATGTATGTCCCTCCGCATCAATGGCACGa                                  >  1:641095/1‑141 (MQ=11)
                                                                                                             aaGGTGCCTGATCTGTAACTACCTTTTGCGGATTCCCCAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTGTGTCCCTCTGTATCCATGGCTCGat                                 >  1:34977/1‑141 (MQ=11)
                                                                                                              aGGTGCCTGAGCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGata                                <  2:138387/141‑1 (MQ=17)
                                                                                                               ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  1:10389/1‑141 (MQ=17)
                                                                                                                       ttCTGTAACTGCCGTTTGAGGCCTACCAAATTGTTTAATGAGACGTTTGATAAACGCTTATGCTGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGTATGTCCCGCTGTATCAATGGCACGATATAAATAGCt                       <  2:312166/140‑1 (MQ=11)
                                                                                                                       aTCTTTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt                       <  1:415054/141‑1 (MQ=17)
                                                                                                                         ttGTAACTACCTTTTGAGGTTTACCAAATTGTTAAATGAGACGTATGCGAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATATCTcc                     <  1:618049/140‑1 (MQ=11)
                                                                                                                         cTGTAACTACCTTTTGAGGCTTACTAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGGATGTCCCTCTGTATCAATTGCTCGATATAAATAGCTcc                     <  1:508482/141‑1 (MQ=11)
                                                                                                                         cTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATATACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTcc                     <  1:171654/141‑1 (MQ=255)
                                                                                                                          tGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                    <  2:123373/141‑1 (MQ=17)
                                                                                                                          tGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                    <  1:741035/141‑1 (MQ=17)
                                                                                                                           gTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAACTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAt                   <  1:602760/141‑1 (MQ=17)
                                                                                                                            tAACTACCTTTTGAGGCTTACCAAATAGTTTAATGAGACGTTTGATAAACGCATATGCTGAATCATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGAACGCTATAAATAGATCCact                  >  2:545969/1‑139 (MQ=11)
                                                                                                                              actACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTAAGCAAGAAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                >  1:661838/3‑141 (MQ=11)
                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:319143/141‑1 (MQ=17)
                                                                                                                                  ccTTTTGTGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:255997/141‑1 (MQ=17)
                                                                                                                                  ccTTTTGAGTTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:454091/141‑1 (MQ=17)
                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTGACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAACGGCACGATATAAATAGCTCCATTTTCCtt            <  1:158133/141‑1 (MQ=255)
                                                                                                                                   cTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           >  2:723676/1‑141 (MQ=17)
                                                                                                                                   cTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCAATTTTCCttt           >  2:430897/1‑141 (MQ=12)
                                                                                                                                    ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtttt          <  1:595191/141‑1 (MQ=17)
                                                                                                                                    ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGACCCTCTGCATCAATGGCACGATATAAATGGCTCCTTTTTCCttgt          >  2:351672/1‑139 (MQ=12)
                                                                                                                                       tGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtt       >  1:290905/1‑141 (MQ=32)
                                                                                                                                          ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATTTAATTTATGTCCCTCTGTATCAGTGGCACGATATAATTAGCTCCATTTTTCttttgttttg    >  1:784249/1‑141 (MQ=11)
                                                                                                                                           gTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGa   >  1:397318/1‑141 (MQ=32)
                                                                                                                                            tttACCAAATTGTTTAATGAAACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt  <  2:894592/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: