Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,161 T→C 32.5% I131V (ATT→GTT)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,1610TC32.5% 530.2 / inf 615I131V (ATT→GTT) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (225/186);  new base C (95/103);  total (322/293)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 630 total aligned reads displayed.

CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301
                                                                                                                                            |                                                                                                                                            
 tttcctACTTTAATATGGCGGTGATGTTGCTCAATGAGGTTTTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                                                             <  2:380129/141‑1 (MQ=11)
   tcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAc                                                                                                                                           >  1:782951/1‑141 (MQ=255)
     ttACTTTAATATGACGGTGATCTTGCTCAATGCGGTTATTCACATATTTCGATGTACAATGGCAGTCAGGTCTAAGGTTAAACGCGTTAATGACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                                                         <  1:401371/141‑1 (MQ=11)
     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                                                         >  2:136182/1‑141 (MQ=255)
     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCt                                                                                                                                         >  1:91474/1‑141 (MQ=255)
       aatgaaaTTTGACGGTGATCTTGCTCGATGAGGGTATTCAGATATTTCGGTGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCTATTGCTGCCTTAGGTGAAGGTGCCTGATATGTGACTACCttt                                                                                                                                       >  2:218027/6‑141 (MQ=11)
       aCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCttt                                                                                                                                       >  2:475081/1‑141 (MQ=255)
         tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                                                     >  2:551829/1‑141 (MQ=255)
         tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTCCCTGATCTATAATTACCtttta                                                                                                                                     >  1:173378/1‑140 (MQ=255)
         tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTg                                                                                                                                     <  2:796358/141‑1 (MQ=255)
          ttAATATGGCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACCATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGa                                                                                                                                    <  2:443785/141‑1 (MQ=11)
              tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAACTTTTGAGGtt                                                                                                                                >  1:385593/1‑141 (MQ=255)
                tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                              <  1:278790/141‑1 (MQ=255)
                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCGACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  2:535516/141‑1 (MQ=255)
                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                             <  1:158105/141‑1 (MQ=11)
                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                            >  2:432054/1‑141 (MQ=21)
                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                            <  2:612785/141‑1 (MQ=25)
                    ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaa                                                                                                                          >  2:68267/1‑141 (MQ=255)
                     gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                         <  1:227300/141‑1 (MQ=11)
                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCCTTGATTGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                        <  1:575538/141‑1 (MQ=11)
                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                     >  1:514919/1‑141 (MQ=17)
                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                     >  2:588910/1‑141 (MQ=12)
                          cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGCCAGTCAGATTTAAGTTTAAAGGCTTTAATGATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                    >  1:371182/1‑141 (MQ=11)
                          cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                    >  1:675770/1‑141 (MQ=11)
                           ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                                                   <  1:309359/141‑1 (MQ=17)
                            tGCTCAATGAGGTTATTCAGATATGTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGATTTAATTACTTTAACCATTGCTACCTTCGTTGGATGTGCCTGATCTGTAATTGCATTTTGGGGTGTACCAAATTTTTTa                                                                                                                  >  1:55351/1‑141 (MQ=11)
                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAATTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                 <  2:33639/141‑1 (MQ=12)
                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTATTTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                 >  1:342773/1‑141 (MQ=11)
                              cTCAATGAGTTTATTCAGAGATTGCTATGTACAATGGCAGTCAGGTTTAAGTTTAATAGCTTTAATTACTTTAACGATTGCTACCGTCGTTGAAGGTGCCTGATCTGGTATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                <  2:51421/141‑1 (MQ=11)
                               tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                               <  1:417273/141‑1 (MQ=17)
                               tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                               >  1:117431/1‑141 (MQ=17)
                               tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAGTGTTTAATg                                                                                                               <  2:183741/141‑1 (MQ=12)
                                     gggTTATTCAGATATTTCGATGTACAATGACAGTTAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGt                                                                                                         <  2:481614/140‑1 (MQ=11)
                                     aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                         >  2:85051/1‑141 (MQ=17)
                                      ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                        <  1:127452/141‑1 (MQ=12)
                                      ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCCTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                        <  2:318784/141‑1 (MQ=11)
                                      ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAACCAATGCTCCCTTCGTTCAAGGTGCCGAATCAGTAGCTGCCTTTTGAGGCTGACCAAATTGTTTAATGAGACGtt                                                                                                        <  2:579015/141‑1 (MQ=11)
                                       gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                       >  2:273123/1‑141 (MQ=12)
                                          aTTCAGATATTTCGATGTACAATGGCAGTCAGCTTTAAGTTTAAAAGCTTTAATTACTTTAACCCTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGCTTTACCAAATTGTTTAATGAGACGCTAGAt                                                                                                    >  2:500316/1‑141 (MQ=11)
                                          aTTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                    >  2:102404/1‑141 (MQ=255)
                                           ttCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                   <  2:235248/141‑1 (MQ=17)
                                           ttCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                   >  1:188480/1‑141 (MQ=255)
                                            tCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAACTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaa                                                                                                  >  1:458182/1‑141 (MQ=12)
                                             cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTGTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                 >  1:835408/1‑141 (MQ=12)
                                             cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                 >  1:296137/1‑141 (MQ=32)
                                             cAGATATTCGGATATTAAATGGCAGTCTTGTTTAAGTTTAACAGCTTTACTTGCTTTAAGCGTTACGACCTTTGTGTATGGTGCCTGATCTGTAATTACCTTTTGATGTTTTCCAAATTGTTACATGAGACGTTTGATaaa                                                                                                 <  2:114827/141‑1 (MQ=2)
                                                    ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                          >  2:761948/1‑141 (MQ=12)
                                                    ttCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                          >  2:391525/1‑141 (MQ=32)
                                                      cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  1:245677/141‑1 (MQ=12)
                                                       gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                       >  2:120935/1‑141 (MQ=32)
                                                         tGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                     <  1:746761/141‑1 (MQ=17)
                                                          gTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGaa                                                                                    <  1:629887/141‑1 (MQ=17)
                                                              aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTGAGCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                <  2:302236/141‑1 (MQ=11)
                                                                tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                              <  1:731950/141‑1 (MQ=11)
                                                                 gACATTCAGGTTTAAGTTTAAAAGCTTTAATGATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCGTGATCTGTAACTAGCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTAt                                                                             <  2:17431/141‑1 (MQ=11)
                                                                                                   tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           <  2:91555/141‑1 (MQ=2)
                                                                                                   tttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           <  2:668500/141‑1 (MQ=11)
                                                                                                      aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTCCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTACGTcc                                        <  1:591834/139‑1 (MQ=9)
                                                                                                       gCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                       <  1:810582/141‑1 (MQ=255)
                                                                                                       gCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                       >  2:819935/1‑141 (MQ=255)
                                                                                                       accATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                       >  1:710332/2‑141 (MQ=17)
                                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGGCCctct                                    >  1:817566/1‑141 (MQ=17)
                                                                                                           ttGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctctt                                   >  2:694183/1‑140 (MQ=255)
                                                                                                             gcgtccTTCGGTGAGGGTGCTTGATGTGCATTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTGGATAAACGCATATGCGGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                 <  1:101720/137‑1 (MQ=11)
                                                                                                             gCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgta                                 <  2:293221/141‑3 (MQ=255)
                                                                                                              ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTAt                                <  2:358404/138‑1 (MQ=255)
                                                                                                              cTACCTTCGTTGAAGTTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                >  1:21075/1‑141 (MQ=17)
                                                                                                              cTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTAt                                >  2:74814/1‑141 (MQ=17)
                                                                                                               tgccTTCGTTGAAGGTGCCTGATCTGTAACTCCCTTTTGATGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                               <  1:442294/139‑1 (MQ=11)
                                                                                                               tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAACTTGTTTAATGAGACGTGTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCGCTCTGTATc                               >  2:210373/1‑141 (MQ=255)
                                                                                                                aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTCTGCATCa                              >  1:135663/1‑141 (MQ=17)
                                                                                                                aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                              <  2:782951/141‑1 (MQ=17)
                                                                                                                aCCTTCCTTGAAGTTGCCTGTTTTGTAATTCACTTTTGAGGTTTACTAAACTGTTTAATTAGACGTTTGATAAACCCATATGCTGAATGTTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                              <  2:451749/141‑1 (MQ=11)
                                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                           <  2:567193/141‑1 (MQ=32)
                                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATCTTTGTCCCTCTGCATCACTg                           >  1:660256/1‑141 (MQ=12)
                                                                                                                     cGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGc                         >  2:124902/1‑141 (MQ=255)
                                                                                                                     cGTTGAACGTGCCTGATCTGTAATTAACGTTTGAGGTTTCCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATATAATGTATGTCCCTCTGCATCAATGGc                         <  1:450303/141‑1 (MQ=11)
                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCa                        >  1:472576/1‑141 (MQ=32)
                                                                                                                        tGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATTGCACg                      >  1:539062/1‑141 (MQ=17)
                                                                                                                        tGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTTTCAATGGCACg                      >  1:484162/1‑141 (MQ=255)
                                                                                                                         gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                     <  1:615099/141‑1 (MQ=32)
                                                                                                                            ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                  <  1:761948/141‑1 (MQ=17)
                                                                                                                            ggTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATCAGCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                  >  2:717147/1‑141 (MQ=11)
                                                                                                                              tGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCCCTGTATCAATGGCACGAAATaa                >  2:166091/1‑141 (MQ=11)
                                                                                                                               gCCTGATGTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATaaa               >  1:522872/1‑141 (MQ=17)
                                                                                                                               gCCTGATCTGGACTTTCCTTTTGAGGTGTACCATATTGTTCAATGAGACGGGTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATaaa               <  1:800931/141‑1 (MQ=11)
                                                                                                                                ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt              >  1:328144/1‑141 (MQ=17)
                                                                                                                                 cTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATa             <  2:408856/141‑1 (MQ=17)
                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCCACTATCTAATGAATGGCCCTCTACATCAATGGCACGCTATAAATac            >  2:418195/1‑140 (MQ=11)
                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATCAAAAg            >  1:641663/1‑141 (MQ=17)
                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAg            >  2:606301/1‑141 (MQ=32)
                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAg            >  1:689541/1‑141 (MQ=32)
                                                                                                                                    atttGTAATTACCTTTTGAGGTTTCCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt          <  2:328359/138‑1 (MQ=17)
                                                                                                                                    atttGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTATTGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt          <  2:629149/138‑1 (MQ=12)
                                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt          <  1:726076/141‑1 (MQ=17)
                                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt          <  2:736135/141‑1 (MQ=17)
                                                                                                                                     tCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc         <  2:651653/141‑1 (MQ=17)
                                                                                                                                         tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAACTAGCTCCCtt     >  2:248301/1‑141 (MQ=12)
                                                                                                                                           gttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAtttt   <  1:261066/140‑1 (MQ=17)
                                                                                                                                            ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc  <  1:741840/141‑1 (MQ=17)
                                                                                                                                            |                                                                                                                                            
CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301

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: