Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,174 T→C 33.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,1740TC33.3% 760.4 / inf 984K126K (AAA→AAGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (323/323);  new base C (116/208);  total (445/539)
Fisher's exact test for biased strand distribution p-value = 2.95e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 1005 total aligned reads displayed.

AATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36033‑36314
                                                                                                                                             |                                                                                                                                            
 atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTACTTACCTTTTGAGGt                                                                                                                                              >  2:438598/1‑141 (MQ=255)
 atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAACTTTAATTACTTTAACCATTGCTACCTTCGTTGAGGGTGCCAGATCTGTAATTACATTTTGAGGt                                                                                                                                              >  2:321736/1‑141 (MQ=11)
  tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGtt                                                                                                                                             >  2:625794/1‑141 (MQ=21)
   atGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGttt                                                                                                                                            <  2:403809/141‑1 (MQ=21)
      tcGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  1:377645/140‑1 (MQ=18)
      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAAGCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  2:695391/1‑141 (MQ=255)
      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  1:106163/141‑1 (MQ=19)
      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCCTTTGAGGTTTAcc                                                                                                                                         >  2:53840/1‑141 (MQ=255)
      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  2:167480/141‑1 (MQ=34)
        ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTATGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTCATTACCTTTTGAGGTTTACCaa                                                                                                                                       >  1:708609/1‑141 (MQ=255)
        ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaa                                                                                                                                       <  1:541971/141‑1 (MQ=11)
         gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGACCTGTAATCACCTTTTGAGGTTTACCaaa                                                                                                                                      <  1:505449/141‑1 (MQ=255)
         gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTGCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      >  1:682082/1‑141 (MQ=11)
         gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGGACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      <  1:582302/141‑1 (MQ=11)
           gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAtt                                                                                                                                    >  1:406723/1‑141 (MQ=17)
           gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAtt                                                                                                                                    <  1:289232/141‑1 (MQ=11)
             tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  1:520826/1‑141 (MQ=17)
             tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  2:740833/1‑141 (MQ=12)
             tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  1:636652/1‑141 (MQ=17)
             tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  1:166249/1‑141 (MQ=32)
              gttGGTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGCTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAGCTTCGTTGAAGGTGCCTGATCTCTAATTACCTTTTGAGGTTTACCAAAGTGtt                                                                                                                                 <  2:574275/140‑1 (MQ=11)
              cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  1:744842/1‑141 (MQ=17)
              cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  2:481569/1‑141 (MQ=32)
              cTTGCACCATGAGGCTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTGTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTCAGGATAACCAAATtttt                                                                                                                                 <  2:111009/141‑1 (MQ=11)
                 gCTCAATGTGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                              <  1:477998/141‑1 (MQ=255)
                  cTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                             <  2:883022/141‑1 (MQ=17)
                    cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           >  1:813856/1‑141 (MQ=17)
                       tGAGGGTATTCAGATATGTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGAACAATGACTTGCTCGTCCCTTTTGAGGCTTACCAAATTGTTTAATGAGAc                                                                                                                        >  2:195502/1‑141 (MQ=11)
                          ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                     <  1:733975/141‑1 (MQ=32)
                            ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  1:943341/1‑141 (MQ=12)
                            ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  1:793025/141‑1 (MQ=12)
                            ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGTAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAACTTGTTTAATGAGACGTTTg                                                                                                                   <  1:567065/141‑1 (MQ=11)
                            ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  2:639514/1‑141 (MQ=255)
                            ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  1:30909/1‑141 (MQ=255)
                             tATTCAGATATTTCGATGTACAATGGCAGTCATGTTTAAGTTTAAAAGCTTTAATTCCATTATCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  1:448796/141‑1 (MQ=11)
                             tATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  1:655667/141‑1 (MQ=17)
                              aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGaa                                                                                                                 >  1:832948/1‑140 (MQ=12)
                              aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                                 >  1:458669/1‑141 (MQ=17)
                               ttCAGATATTTCGATGGACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                                >  1:468857/1‑141 (MQ=12)
                                 cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTATAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              <  1:668093/141‑1 (MQ=12)
                                   gATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:936112/141‑1 (MQ=255)
                                   aatatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:280023/140‑1 (MQ=17)
                                     tatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  2:758916/141‑1 (MQ=12)
                                         tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAAGTTTGATAAACGCATATg                                                                                                      >  2:800686/1‑141 (MQ=12)
                                         tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      >  2:361359/1‑141 (MQ=255)
                                         tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      <  1:359604/141‑1 (MQ=255)
                                           gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGTGACGTTTGATAAAAGCATATGCt                                                                                                    >  2:507597/1‑141 (MQ=255)
                                            gtGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                                   <  2:296153/140‑1 (MQ=11)
                                                aCAATGTCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGATGAAtc                                                                                               <  2:244951/141‑2 (MQ=12)
                                                    tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           <  2:411506/141‑1 (MQ=255)
                                                                                    tattttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgt                                                           <  2:778686/138‑1 (MQ=0)
                                                                                    tactgtaaacATTGCTACCTTCGTTGAAGGTCCCTGATCTGTAACTACCTTTTGAGGTTTACCCAATTGTTTAATGCGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgt                                                           <  1:393291/132‑1 (MQ=1)
                                                                                       tttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCCTTGCTTATGCAACCAAATATCTAatgtatg                                                        >  1:875549/1‑141 (MQ=2)
                                                                                        ttAACCATTGCTACCTTCGTTGAATGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       <  1:215119/141‑1 (MQ=14)
                                                                                         taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAACATCTAATGTATGTc                                                      <  1:493428/138‑1 (MQ=17)
                                                                                          aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  2:869232/141‑1 (MQ=12)
                                                                                          aGCCATTGCTGCCTTCGTTAGAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACACAACCAAATATCTAATGTATGTcc                                                     <  1:600493/141‑1 (MQ=11)
                                                                                            ccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGATTGATTATTTCGTTGCTTATGCAACGAAGTATCTAATGTATGTctct                                                   >  1:267343/1‑141 (MQ=11)
                                                                                             cattgtgatcTTGGTTGAAGGTGCCTGATCTTTAATTACCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGATAACCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  1:34912/132‑1 (MQ=11)
                                                                                             cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTCTGTccggc                                                  <  2:816210/141‑4 (MQ=17)
                                                                                             cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  2:5973/141‑1 (MQ=32)
                                                                                             cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  2:727352/141‑1 (MQ=32)
                                                                                               ttGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTg                                                >  2:228065/1‑141 (MQ=11)
                                                                                                tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                               >  1:864622/1‑140 (MQ=17)
                                                                                                 tcTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                              <  2:389014/140‑1 (MQ=17)
                                                                                                  ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGTGGCTGCAGGTCGTGGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             >  1:923128/4‑138 (MQ=11)
                                                                                                   tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCAGATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAAGGTATGTCCCTCTGTATc                                            >  1:608365/1‑141 (MQ=12)
                                                                                                    aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAACTATCTAATGTATGTCCCTCTGTATCa                                           <  2:523849/141‑1 (MQ=17)
                                                                                                    aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           >  1:624478/1‑141 (MQ=17)
                                                                                                     ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:178206/141‑1 (MQ=255)
                                                                                                     ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTATTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          >  2:668002/1‑141 (MQ=255)
                                                                                                     ccTTCGTTGAAGGTGCCTGATCTGGAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:208183/141‑1 (MQ=255)
                                                                                                      ctgcggtGAAGGTGCCTTATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         <  2:460441/135‑1 (MQ=255)
                                                                                                      cTTCGGTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTATGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         <  2:683227/141‑1 (MQ=255)
                                                                                                       ttCGTGGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                                        <  1:281026/141‑1 (MQ=17)
                                                                                                          gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGcc                                     >  2:838263/1‑140 (MQ=17)
                                                                                                           ttGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    >  1:244375/1‑141 (MQ=32)
                                                                                                           tcgAAGTTGCATGATTTGTAATTACCTTTTGACGTTTACTACATTGTTTAATGCGACGTTTGATACACGCATATGCGTAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:321736/139‑1 (MQ=11)
                                                                                                            tGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATTAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAAGATCTAATGTATGTCCCTCTGCATCCATGGCACg                                   >  1:129055/1‑141 (MQ=12)
                                                                                                            tGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACg                                   >  1:910536/1‑141 (MQ=17)
                                                                                                            tGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCGCTCTGTATCAATGGCACg                                   >  1:189830/1‑141 (MQ=255)
                                                                                                                   gCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATaaa                            >  1:413758/1‑141 (MQ=32)
                                                                                                                    ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATCTGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                           >  1:436979/1‑141 (MQ=255)
                                                                                                                    ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                           >  1:853246/1‑141 (MQ=12)
                                                                                                                       gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  2:30909/141‑1 (MQ=12)
                                                                                                                        acttGTAACTACCTTTTGAGGCTTACCAAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  1:230531/138‑1 (MQ=11)
                                                                                                                        aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt                       <  2:894246/141‑1 (MQ=32)
                                                                                                                          ttGTAATTACCTTTTTATGTTTACCAAATTGTTTAATGAGACGTTTGATAACCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                     <  1:780956/140‑1 (MQ=12)
                                                                                                                              aaTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAttt                 <  2:787996/141‑1 (MQ=32)
                                                                                                                               aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                >  2:709670/1‑141 (MQ=17)
                                                                                                                               aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTACGCAACCAAATATCTAATGTATGTCCCTCTGGATCAATGGGACGATATAAATAGCGCCActtt                >  2:913769/1‑141 (MQ=11)
                                                                                                                                  tccTTGTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:150237/140‑1 (MQ=255)
                                                                                                                                  aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  2:143432/141‑1 (MQ=17)
                                                                                                                                  aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAACGGCACCACATAAATAGcccaatctccct             >  2:56186/1‑130 (MQ=11)
                                                                                                                                  aCCTTTTGAGGCTGACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACAAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  2:848576/141‑1 (MQ=255)
                                                                                                                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCCTATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATCTAATGTATGTCCCTCTGTATCCATGCCCCGATATAAAACGCCCCCTTGTCCttt           >  2:128111/1‑141 (MQ=11)
                                                                                                                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt           >  1:893515/1‑141 (MQ=32)
                                                                                                                                      tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttttt         >  2:350328/1‑140 (MQ=17)
                                                                                                                                         gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttt      >  1:198646/1‑141 (MQ=17)
                                                                                                                                          aGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtttt     >  1:75830/1‑141 (MQ=17)
                                                                                                                                             |                                                                                                                                            
AATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36033‑36314

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