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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,1740TC30.3% 569.1 / inf 625K126K (AAA→AAGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (233/195);  new base C (85/102);  total (322/303)
Fisher's exact test for biased strand distribution p-value = 4.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-02

Only 100 of 642 total aligned reads displayed.

TATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36035‑36314
                                                                                                                                            |                                                                                                                                             
tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGtt                                                                                                                                              >  2:119578/1‑141 (MQ=35)
 atGACGGTGATCTTGCTCAAGGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGttt                                                                                                                                             <  2:813996/141‑1 (MQ=21)
  tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                                            <  1:117907/141‑1 (MQ=25)
  tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                                            >  2:835126/1‑141 (MQ=255)
    gcGGTGATCTTGCTCAATGAGTTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                                          <  1:189159/140‑1 (MQ=11)
    aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                          <  2:695532/141‑1 (MQ=21)
    aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                                                                          >  1:265534/1‑141 (MQ=255)
    aCGGTGATCTTGCTAAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGCGGTTTAcc                                                                                                                                          >  1:457263/1‑141 (MQ=255)
     cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCa                                                                                                                                         <  2:109052/141‑1 (MQ=11)
          aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAATTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAttt                                                                                                                                    >  2:504416/1‑140 (MQ=12)
           tcacgCTCAATGAGGTAATAACAATGTTGCGATGTACAACGGCAGTCAGGTTTAAGTTTCAAAGGTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGAGCCTGATCTGTAATTACCTTTGGAGGATTACCAAATTGt                                                                                                                                   >  2:385761/5‑141 (MQ=11)
           tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                   >  1:514919/1‑141 (MQ=17)
           tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                   >  2:742396/1‑141 (MQ=12)
           tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCGGGTTTAAGTTTAAGAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAATGTGCCTGATCTGTGACTACCTTTTGAGGCTTAGCAAATTGt                                                                                                                                   >  2:408084/1‑141 (MQ=11)
            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCATGGTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                  >  1:788821/1‑141 (MQ=255)
            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:305584/1‑141 (MQ=255)
            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:675770/1‑141 (MQ=11)
            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAACCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:618625/1‑141 (MQ=11)
            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAgttgtt                                                                                                                                  >  2:609937/1‑141 (MQ=17)
               gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGGTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                               >  1:73176/1‑141 (MQ=255)
               gCTCAATGAGGTTATTCAGATAGTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                               <  1:9171/141‑1 (MQ=12)
                  cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                            <  1:520892/141‑1 (MQ=12)
                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                                            >  1:293221/1‑141 (MQ=11)
                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                                            <  2:792591/141‑1 (MQ=11)
                   aaTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                                                                           <  2:324384/141‑1 (MQ=32)
                      gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGGGACg                                                                                                                       >  2:5758/1‑141 (MQ=12)
                       atgTTATGCATATATTTCGATGTACAATGACAGTCATGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTTAAGTTGCCTATTCTGTAACTTCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  2:602823/139‑1 (MQ=11)
                        ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAGGGTGGCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                                     <  2:341094/141‑1 (MQ=11)
                         gTTATTCAGATATTTCGATGTACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  2:457499/141‑1 (MQ=12)
                         gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    >  1:348580/1‑141 (MQ=17)
                          ttagtgagagATTTCGATGTACAATGACAGTGAGATTTAAGTTTAAAAGCTTTATTTATTTTCGCCATTGCTGCCTTTGTTGAAGGTGCCTGATCTGTAACTAGCTTTGGAGGCTTACGAAATTGTAAAATCAGACGTTTg                                                                                                                   <  2:183703/135‑1 (MQ=2)
                          ttATTCGAATATTTCGATTTACAAGGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  2:444276/141‑1 (MQ=11)
                          ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  2:362691/141‑1 (MQ=17)
                               cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTTACCATTTCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTGGATatg                                                                                                              >  1:81684/1‑139 (MQ=11)
                                    atTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         >  2:160577/1‑141 (MQ=12)
                                     ttgcGATGTAAAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAACTACTTTAACCATTGCTACCTTCGGTGAAGGTGCCGGATCTGGAACTACCTTTTGAGGGTTACCAAATTGTTTAATGAGACGTTTGATAAACTCata                                                                                                        <  2:215974/138‑1 (MQ=11)
                                     tcacGACGTACCATGACAGTCAAGTTTAAGTTTAAAAGCTTTAAATACTTTAGCCATTGCTAACTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGACACGTTTGATAAACGCata                                                                                                        <  2:752406/138‑1 (MQ=11)
                                      ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTGTGACCATTGCTACCTTCGTTGAAGGTGACTGATCTGTAATTACCTTTTGATGTTTCCCTAATTGTTTAAAGAGTACGttattaaatggatat                                                                                                        >  2:23086/1‑128 (MQ=11)
                                        cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACGTTTTGAGGTTTACCAAATTGTTTAATGAAGCGTTTGATAAACGCATAAGc                                                                                                     >  2:548021/1‑141 (MQ=11)
                                        cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  2:778947/141‑1 (MQ=17)
                                         gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                                    >  1:481828/1‑141 (MQ=255)
                                          aTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                                   >  1:2328/1‑141 (MQ=12)
                                            gTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGaa                                                                                                 <  1:419923/141‑1 (MQ=32)
                                                aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             <  1:717763/141‑1 (MQ=14)
                                                  tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           >  1:443034/1‑141 (MQ=255)
                                                   gACAGGCCGGGTTAAGTTTAAAAGCTGTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACTTTTTGAGGGTTACCAAATTGTTTAATGAGAGGTTTGATAAACGCATCTGCtaccggattat                                                                                          >  2:260001/1‑131 (MQ=11)
                                                                                     tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        <  2:91474/141‑1 (MQ=2)
                                                                                     tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        <  1:383292/141‑1 (MQ=2)
                                                                                       taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                                      <  1:804006/138‑1 (MQ=17)
                                                                                        aaccATTGCTATCTTCGTTGAAGGCGCCTGATCTGTAATTACCTTTTAAGGTTTACCAAATTGTTTAATGAGACGTTTGAGCAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  2:300272/139‑1 (MQ=11)
                                                                                        aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTCCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTACGTcc                                                     <  1:591834/139‑1 (MQ=9)
                                                                                        aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     >  2:332934/3‑141 (MQ=17)
                                                                                        aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTACGCAACCAAATATTTAaagtatgtac                                                     >  2:43232/1‑139 (MQ=11)
                                                                                           cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAACTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  >  1:741243/1‑141 (MQ=17)
                                                                                            aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTTTCTCGTTGCTTACGCAAACAAATATCTAATGTATGTCCctct                                                 >  1:171141/1‑141 (MQ=17)
                                                                                            aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAAGCAAATATCTAATGTATGTCCctct                                                 >  2:523133/1‑141 (MQ=17)
                                                                                             ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                                <  1:579765/141‑1 (MQ=32)
                                                                                             ttGCTACCTTCGTTGAAAGTGCCTCACCTCAGGTGACGGTTTGAGGGTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTCACTTACGCAACAAAATATCTAATCTAAGTCCctcta                                                <  2:173588/141‑2 (MQ=11)
                                                                                              tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               <  2:136123/141‑1 (MQ=32)
                                                                                              tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               >  2:277428/1‑141 (MQ=32)
                                                                                              tGCTACCTTCGTTGAAGGTGCCTGGCGATAGAGGACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               >  1:830313/1‑141 (MQ=11)
                                                                                               tgtggcTTCGTTGAAGGTGCCTGATCCGTAACTACCTTTCGAGGCTTACCAAATTGTTTATTGAGTAGTTTGATAAACGCATATGCTGAATTATTATCTCGTTGCTTACGCAGCCAAATATCTAATGTATGACCCTCTgta                                              <  2:67873/136‑3 (MQ=11)
                                                                                               gcgtccTTCGGTGAGGGTGCTTGATGTGCATTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTGGATAAACGCATATGCGGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                              <  1:101720/137‑1 (MQ=11)
                                                                                                ttACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             <  2:582052/140‑4 (MQ=17)
                                                                                                cTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                             <  1:252591/141‑1 (MQ=32)
                                                                                                 tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            >  1:690943/1‑141 (MQ=32)
                                                                                                  gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTCTCTCGTTGCTTACGCAACCAGATATCTAATGTATGTCCCTCTGTATCa                                           >  2:668396/2‑141 (MQ=11)
                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTCTGCATCa                                           >  1:135663/1‑141 (MQ=17)
                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  2:115843/141‑1 (MQ=17)
                                                                                                   ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  1:631987/141‑1 (MQ=12)
                                                                                                     ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCACATTTTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  1:373940/141‑1 (MQ=12)
                                                                                                     ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                                        <  1:551829/141‑1 (MQ=32)
                                                                                                       cGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGc                                      >  2:124834/1‑141 (MQ=12)
                                                                                                         ttGAAGGTGGCTGCTCTGTAATTAGCTTTTGAGTTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:316545/141‑1 (MQ=255)
                                                                                                         ttGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                                    >  1:754477/1‑141 (MQ=17)
                                                                                                          tGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACg                                   <  2:558472/141‑1 (MQ=12)
                                                                                                             aGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGata                                >  2:648720/1‑141 (MQ=17)
                                                                                                              ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAGTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  2:623405/1‑141 (MQ=17)
                                                                                                              ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  1:665852/1‑141 (MQ=17)
                                                                                                              ggTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  1:772190/1‑141 (MQ=12)
                                                                                                                  ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAGATAAAt                           >  2:310603/1‑141 (MQ=17)
                                                                                                                  ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCACTCTGCATCAATGGCACGATATAAAt                           >  2:686619/1‑141 (MQ=17)
                                                                                                                    tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCCACTATCTAATGAATGGCCCTCTACATCAATGGCACGCTATAAATac                         >  2:418195/1‑140 (MQ=11)
                                                                                                                      aTCTGTAAATACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  1:746443/141‑1 (MQ=17)
                                                                                                                              actACCTGTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTATTTTATGCTCCTCGGTATCAATGGCACGATATAAATAGCTCCAtttt                <  1:538199/139‑1 (MQ=11)
                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                >  1:331475/1‑141 (MQ=17)
                                                                                                                                tACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              >  1:651498/1‑141 (MQ=17)
                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTAACGCTGCATAAATGTCACGATATAAATAGATCCATTTTGCt             >  2:759318/1‑141 (MQ=11)
                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:730044/141‑1 (MQ=17)
                                                                                                                                  ccTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  1:127401/141‑1 (MQ=17)
                                                                                                                                   ctgttcaGGTCTTCCAAATTGTTTAATGAGACGTTTGATAAACGAATATGCTGAATTATTATCTCGTTGCTTACGCAACCAAATCTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  2:474933/135‑1 (MQ=11)
                                                                                                                                   cTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  1:177128/141‑1 (MQ=17)
                                                                                                                                   cTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           >  1:220048/1‑141 (MQ=12)
                                                                                                                                     tttGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTa         <  2:481828/141‑1 (MQ=12)
                                                                                                                                      ttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAt        >  1:464637/1‑141 (MQ=17)
                                                                                                                                       tGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTCCGCAACCAAATTTCTGTTGGAGGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtt       <  1:401333/141‑1 (MQ=11)
                                                                                                                                        gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAATTAAATAGCTCCATTTTCCTTTTAttt      >  1:141752/1‑141 (MQ=12)
                                                                                                                                          ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTg    <  2:246837/141‑1 (MQ=32)
                                                                                                                                          ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCCCCATTTTCCTTTTATTTTg    >  2:201900/1‑141 (MQ=17)
                                                                                                                                          ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTACCTCTGCATCAATGGCACGATATAAGTAGGTCCATTTTCCTTTTATTTTg    >  1:708341/1‑141 (MQ=12)
                                                                                                                                            |                                                                                                                                             
TATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36035‑36314

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