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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,1610TC35.0% 418.6 / inf 589I131V (ATT→GTT) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (179/198);  new base C (108/96);  total (293/296)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 600 total aligned reads displayed.

CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301
                                                                                                                                            |                                                                                                                                            
 tttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                                                             <  1:753632/141‑1 (MQ=255)
   tcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGAACTGTAATTAc                                                                                                                                           >  1:629934/1‑141 (MQ=255)
     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCt                                                                                                                                         <  2:308794/141‑1 (MQ=255)
       tcTTTAATAGGACGGTTATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTCCGTTGAAGGTGCCTGATCTGTAATTACCttt                                                                                                                                       <  1:600733/140‑1 (MQ=11)
       aCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCttt                                                                                                                                       >  2:460768/1‑141 (MQ=255)
        ttttAATAGGATGGGTATCTTGCTCAATGAGGTTATTCAGATATTTCGGTGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtttt                                                                                                                                      <  2:287414/140‑1 (MQ=11)
        cTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtttt                                                                                                                                      <  1:508412/141‑1 (MQ=255)
        cTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCGAGTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCtttt                                                                                                                                      <  1:729410/141‑1 (MQ=11)
        cTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCtttt                                                                                                                                      >  1:223041/1‑141 (MQ=255)
        cTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtttt                                                                                                                                      <  2:418726/141‑1 (MQ=35)
         tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTCCCTTTTg                                                                                                                                     >  2:558030/1‑141 (MQ=255)
           tAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAg                                                                                                                                   >  2:588121/1‑141 (MQ=25)
           tAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAg                                                                                                                                   >  1:660032/1‑141 (MQ=255)
            aaTATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAgg                                                                                                                                  <  1:328522/141‑1 (MQ=255)
             atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTTTAACTACCTTTTGAggc                                                                                                                                 >  2:57359/1‑140 (MQ=255)
              tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGCTGTACAATGACAGTCCGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAggct                                                                                                                                >  1:49006/1‑139 (MQ=11)
              tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGtt                                                                                                                                >  1:357052/1‑141 (MQ=25)
               atGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGctt                                                                                                                               >  1:535651/1‑141 (MQ=255)
                tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                              >  2:562829/1‑141 (MQ=255)
                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGGCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                            >  2:483134/1‑141 (MQ=11)
                      tGATCTTGCTCAATGAGTTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAt                                                                                                                        <  2:168222/141‑1 (MQ=11)
                        aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAttt                                                                                                                      >  2:618476/1‑140 (MQ=12)
                        aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                                                      <  2:258789/141‑1 (MQ=17)
                        aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTg                                                                                                                      >  1:645018/1‑141 (MQ=255)
                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                     >  2:369861/1‑141 (MQ=12)
                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTGTGAGGTTTTCCCAGTTGt                                                                                                                     >  2:724397/1‑141 (MQ=11)
                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGt                                                                                                                     >  2:667415/1‑141 (MQ=255)
                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTAATCTGTAACTACCTTTTGAGGATTACCCAATTGt                                                                                                                     >  2:688648/1‑141 (MQ=11)
                          cTTGCTCAATGAGGTTTTTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATATGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                    >  1:384264/1‑141 (MQ=17)
                          cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                    >  2:222265/1‑141 (MQ=17)
                          cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                    >  1:392852/1‑141 (MQ=255)
                           ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                                                                   >  2:359990/1‑141 (MQ=12)
                            tGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTa                                                                                                                  >  1:165780/1‑141 (MQ=255)
                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                 >  2:745962/1‑141 (MQ=12)
                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTTCTGCCTTAGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTat                                                                                                                 >  1:368225/1‑140 (MQ=11)
                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTTCCAAATTGTTTat                                                                                                                 >  1:253246/1‑140 (MQ=17)
                               tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                               <  1:301081/141‑1 (MQ=12)
                                cAATGAGGTTATTCCGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaaata                                                                                                              <  2:710864/141‑4 (MQ=11)
                                   tGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                                                                                                           >  1:686423/1‑141 (MQ=17)
                                     gggTTATTTAGATTTTTTGATGTACATTGGCTGGCAGTTTTATGTTTAAAAGGTTTAATTGCTTTATCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                         <  2:408221/140‑1 (MQ=11)
                                     aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                         >  1:669310/1‑141 (MQ=17)
                                         tATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                     <  1:151919/141‑1 (MQ=255)
                                          aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAGCAAATTGTTTAATGAGACGTTTGAt                                                                                                    >  1:617508/1‑141 (MQ=12)
                                          aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATAACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGTTCTGTAATTAACTTTTGCGGTTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                    >  1:517027/1‑141 (MQ=11)
                                           ttCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                   <  2:82691/141‑1 (MQ=255)
                                             cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGGTaaa                                                                                                 >  2:104717/1‑141 (MQ=255)
                                             cAGATAGTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                 <  2:471972/141‑1 (MQ=11)
                                                 tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                             <  1:770659/141‑1 (MQ=11)
                                                 tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                             <  1:548766/141‑1 (MQ=32)
                                                    ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCctat                                                                                          >  1:89789/1‑141 (MQ=12)
                                                     tCGATGTAAAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                         <  2:646366/141‑1 (MQ=11)
                                                      cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  1:238339/141‑1 (MQ=17)
                                                      cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTGTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAAGCATATGc                                                                                        >  2:244957/1‑141 (MQ=255)
                                                      cGATGTACAATGGCAGTCAGGATTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        >  1:401648/1‑141 (MQ=255)
                                                      cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  1:582766/141‑1 (MQ=11)
                                                         tGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                     >  2:297480/1‑141 (MQ=11)
                                                          gTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGaa                                                                                    <  1:415137/141‑1 (MQ=17)
                                                            aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                  >  1:677318/1‑141 (MQ=12)
                                                            aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                  <  1:563740/141‑1 (MQ=17)
                                                             cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                 <  1:47646/141‑1 (MQ=32)
                                                             cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                 <  1:696369/141‑1 (MQ=32)
                                                              aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                <  1:104700/141‑1 (MQ=14)
                                                              aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                >  2:193249/1‑141 (MQ=11)
                                                               aTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTGGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGCCTTACCAACTTGTTTAATGAGACGTTTGATAAACGCGTATGCTGAATGAtt                                                                               >  2:451907/1‑141 (MQ=11)
                                                                                       agatttACTTATTTCACCCATTGCTGCCTGCGATGAAGGTGCCTCATCTGCAACTACCTTTTGAGGCTTACCAAATTGTTTAAGGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAat                                                       <  2:531642/138‑1 (MQ=1)
                                                                                                   tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           >  1:18233/1‑141 (MQ=2)
                                                                                                   tttAACCAGTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           <  2:763168/141‑1 (MQ=11)
                                                                                                     ttgacATGGCTGCCTTCCGTGAAGGTGCCTAAACTGTAACTACCTTTTGAGGCTTAGCAAATTGTTTAACGAGACGTTAGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  2:1307/137‑1 (MQ=11)
                                                                                                     taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  1:139739/138‑1 (MQ=17)
                                                                                                     taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  1:732947/138‑1 (MQ=17)
                                                                                                     taaccATTGCAACCTTGTTTGAGGGTGCCTGATCGTCAATACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        <  1:346474/138‑1 (MQ=11)
                                                                                                        ccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAGTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCct                                      <  2:509129/141‑1 (MQ=17)
                                                                                                         cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  2:755074/141‑1 (MQ=32)
                                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTACTTACCGTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTAACGCAACCAAATATCTAACGGATGTcccccc                                    >  2:118440/1‑138 (MQ=11)
                                                                                                           ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCCCTg                                   >  2:607409/1‑141 (MQ=17)
                                                                                                           ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCACTg                                   >  2:628133/1‑141 (MQ=17)
                                                                                                             cgtACCTTCGTTGACGGTTCCTCATCTGTAATTCCCTTTTGAGGTTTACCAACTTGTTTAAGGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                 <  1:275860/139‑1 (MQ=11)
                                                                                                              cTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAACTATCTAATGTATGTCCCTTTGCAt                                >  1:292007/1‑141 (MQ=17)
                                                                                                               tgccTTCGTTGAAGGTGCCTGAGCTGTAACTAGCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                               <  2:273011/139‑1 (MQ=11)
                                                                                                                gccTTCGTTGAAGGTGCCTGATCTGTAACTACATTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                              <  1:196072/140‑1 (MQ=255)
                                                                                                                aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATCATCACTGAAACGCAATGATTATCTCGTTGCTTACGCAACCAAATATATAATGTATGTCCCTCTGCATCa                              >  1:761582/1‑141 (MQ=11)
                                                                                                                 ccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCaa                             <  2:679993/141‑1 (MQ=32)
                                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGGAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                           >  1:428992/1‑141 (MQ=17)
                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCACTGGCa                        >  2:400619/1‑141 (MQ=17)
                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATTGCa                        >  2:176829/1‑141 (MQ=17)
                                                                                                                       ttGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                       <  2:638971/141‑1 (MQ=17)
                                                                                                                        tGAAGGTGTCGGATCTGTAATTTCCTTTTGAGGTTTACCGAATTGTTTCATGAGACGTTTGCTAAGCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACg                      <  1:113942/141‑1 (MQ=11)
                                                                                                                        tGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGGTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATATAATGTATGTCCGTCTGTATCAATGGCAca                      >  2:386207/1‑140 (MQ=11)
                                                                                                                          aaGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGat                    >  1:706705/1‑141 (MQ=32)
                                                                                                                            ggTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                  >  2:543458/1‑141 (MQ=12)
                                                                                                                                ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTACGCAACCAAATCTCTAATGTATGTCCCTCTGTGTCAATGGCACGTTATAAAt              >  1:213547/1‑141 (MQ=11)
                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGTACGATATAAATAg            >  2:658502/1‑141 (MQ=17)
                                                                                                                                  tGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATACAATGTATGTCCCTCTGTATCAATGGCACGAAATAAATAg            >  2:393317/1‑141 (MQ=11)
                                                                                                                                   gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATTTATGTCCCTCTGTATCAATGGCACGAGATAAATAgg           >  2:144938/1‑140 (MQ=12)
                                                                                                                                   gATCCGTAACTACCTTTTTAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc           <  2:331681/141‑1 (MQ=11)
                                                                                                                                      cTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTTTATCAATGGCACGGTATAAATAGCTcc        >  2:429601/1‑141 (MQ=12)
                                                                                                                                       cgTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa       <  1:342485/140‑1 (MQ=17)
                                                                                                                                        gTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAt      <  2:344243/141‑1 (MQ=17)
                                                                                                                                            ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc  <  1:462753/141‑1 (MQ=17)
                                                                                                                                            ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTc  >  1:261466/1‑141 (MQ=32)
                                                                                                                                            |                                                                                                                                            
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: