Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,997 T→C 31.3% T185T (ACA→ACG USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,9970TC31.3% 1153.8 / inf 1009T185T (ACA→ACGUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (347/342);  new base C (160/154);  total (513/496)
Fisher's exact test for biased strand distribution p-value = 8.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1012 total aligned reads displayed.

CAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTT  >  CP000730/35867‑36137
                                                                                                                                  |                                                                                                                                             
cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAGGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                                                                                     <  2:346327/141‑1 (MQ=2)
cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                                                                                     <  1:357173/141‑1 (MQ=2)
 aGTGTTCGTTTAACTTGTTAGGATGATGCTAATTTCGTGGCATGGTGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt                                                                                                                                    <  2:217599/141‑1 (MQ=11)
 aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                                                                                                                    <  1:336484/141‑1 (MQ=11)
 aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                                                                                                                    <  1:521340/141‑1 (MQ=255)
  gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                                                   >  1:430628/1‑141 (MQ=17)
       cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                                                                                                              >  1:499051/4‑141 (MQ=255)
       cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                                                                                                              >  1:139743/4‑141 (MQ=255)
        gtttAGCTGGCTAGCACGATGCTAATGTCGTGGCATGCGGAAAATCCGTAGATCTTAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATat                                                                                                                             <  1:507247/139‑2 (MQ=11)
         tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                                                                            <  1:571000/140‑1 (MQ=17)
            aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                                                                                                         <  2:539906/141‑1 (MQ=17)
               tggCTAGCATGATGCTAATTTCCTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGTTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                                                                                      <  2:422744/139‑1 (MQ=12)
                tGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctt                                                                                                                     <  1:148359/141‑1 (MQ=11)
                tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                     <  1:101021/141‑1 (MQ=17)
                tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                     <  1:365024/141‑1 (MQ=17)
                tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                     >  1:418161/1‑141 (MQ=17)
                 gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                                                                                                                    >  2:646660/1‑141 (MQ=255)
                 gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                                                                                                                    >  2:126053/1‑141 (MQ=255)
                 gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                                                                                    >  1:240026/1‑141 (MQ=17)
                  cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttt                                                                                                                   >  2:392307/1‑140 (MQ=17)
                  cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                                                   <  2:471382/141‑1 (MQ=17)
                   tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                                                  <  2:28014/141‑1 (MQ=32)
                   tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGGATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                                                  >  1:527407/1‑141 (MQ=17)
                       atgatgATAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTTCCGTATTGATACTTTGATATCTTGTCTTTCTTACt                                                                                                              <  2:608502/141‑1 (MQ=11)
                        tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                                                                                                             >  1:439158/1‑141 (MQ=11)
                         gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                                                                                                            <  2:291389/141‑1 (MQ=25)
                          atgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTa                                                                                                           <  1:111775/141‑1 (MQ=11)
                             cTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTCTATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATAGCGTGTCTTTCTTGCTTTAata                                                                                                        >  1:448042/1‑141 (MQ=11)
                              tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                                                                                                       >  2:511774/1‑141 (MQ=25)
                                aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                                                                                                     >  1:250744/1‑141 (MQ=35)
                                     ggggCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTACATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                                <  1:218834/139‑1 (MQ=14)
                                     gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                                <  1:78893/141‑1 (MQ=255)
                                     gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                                <  2:635871/141‑1 (MQ=255)
                                        gCATGTCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc                                                                                            <  2:595035/141‑1 (MQ=255)
                                          aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCtt                                                                                          <  1:564041/141‑1 (MQ=255)
                                           tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTTCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                         <  2:58702/141‑1 (MQ=255)
                                           tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                         <  2:496685/141‑1 (MQ=38)
                                           tGGCGAAAATCCGTAGATATGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATATTg                                                                                         <  2:174069/141‑1 (MQ=255)
                                              cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc                                                                                      >  1:275268/1‑141 (MQ=38)
                                                aaaaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTCAGTCCGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCaa                                                                                    <  1:288898/141‑1 (MQ=14)
                                                 aaaTCCGTAGATCTAAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGACTTTCTTACTTTAATATGACGGTGATCTTGCTCAAt                                                                                   >  1:369572/1‑141 (MQ=255)
                                                 aaaTCCGTAGAACTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGCTACTTTGCTATCTTGTCTTTCTTCCTTTAATATGACGGTGATTTTGCTCAAt                                                                                   >  2:376348/1‑141 (MQ=11)
                                                  aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                                                  <  2:3164/141‑1 (MQ=255)
                                                   aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGa                                                                                 >  1:577370/1‑141 (MQ=37)
                                                   aTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGa                                                                                 >  2:148379/1‑141 (MQ=255)
                                                      cGTAGATCTGAAGAGACCTGCGGTTCATTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTCATACTTTGATACCTTGTCTTTCTTACTTTAATATGACTGTGATCTTGCTCAATGAGGt                                                                              <  1:403134/141‑1 (MQ=255)
                                                      cGTAGAACTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                                                                              >  2:368680/1‑141 (MQ=255)
                                                       gTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGtt                                                                             <  1:659088/141‑1 (MQ=255)
                                                        tAGATCTGAAGAGACCTACGGTTCTTGTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTa                                                                            <  2:158951/141‑1 (MQ=255)
                                                            tCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAttaa                                                                        >  2:246451/1‑139 (MQ=37)
                                                            tCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTAATACTTTGATCCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCa                                                                        >  1:192628/1‑141 (MQ=255)
                                                              tGAAAAGACCTGCGGTTCTTTTTATATTGTGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGTTACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa                                                                      <  2:567056/141‑1 (MQ=255)
                                                               gAAGAGACCTACGGTTCTTTTTATATTTAGTTTAAATACATTCAATACCTTTTAAAGTATTGTTTGCCGTCTTGATACTTTGTTATCTTTTCTTTCTTCCTTTAATATGACGTGAGACTTGCTCAATATTGTTGTTCAGag                                                                     >  2:409436/1‑140 (MQ=2)
                                                                aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCatagt                                                                    <  2:79771/141‑5 (MQ=255)
                                                                aaGAGACCTACGATACTGTTTATATAGAGCGTAAATACATTCAATGCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATGTTACGGTGGTCTTGCTCAATGAGAGTATTCAGata                                                                    >  2:226368/1‑141 (MQ=11)
                                                                   agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAttt                                                                 >  2:485117/1‑141 (MQ=255)
                                                                    gaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGTTCTTTCTCAATGAGGTTATTCATATAtttt                                                                >  1:317812/1‑140 (MQ=255)
                                                                      ccTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGa                                                              <  2:62316/141‑1 (MQ=38)
                                                                      ccTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTGGCTCAATGAGGTTATTCAGATATTTCGa                                                              >  1:123105/1‑141 (MQ=255)
                                                                        tacGGTTCTTTTTATATAGAGCGTAAATACAATCAACACCTTCTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATg                                                            <  1:16470/139‑1 (MQ=11)
                                                                            gTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACa                                                        <  2:269120/141‑1 (MQ=38)
                                                                              tCTTTTTATATAGAGCGTAAATAACATCAATACCTTTTAAATTATTCTTTTCTTTATTGATACTTTGATACCTTGTCTTTCTTACTTGAATATGACGGTTATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                      <  2:628657/141‑1 (MQ=11)
                                                                                 ttttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGAc                                                   >  1:626129/1‑141 (MQ=255)
                                                                                  tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGGTATTTCGATGTACCATGGCa                                                  <  2:522653/141‑1 (MQ=255)
                                                                                    tatataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGt                                                <  1:380637/141‑1 (MQ=255)
                                                                                     atataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTc                                               >  1:90334/1‑141 (MQ=255)
                                                                                        aagagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAgg                                            <  1:603517/140‑1 (MQ=255)
                                                                                         agagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTTATACTTTGCTACCTTGTCTTTCTTACTTTAATATGACGGTGAGCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgt                                           >  2:434277/1‑141 (MQ=255)
                                                                                           ggcttAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttt                                         <  2:657031/137‑1 (MQ=255)
                                                                                           agCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttt                                         >  1:390023/1‑141 (MQ=255)
                                                                                            gCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttta                                        <  2:360809/141‑1 (MQ=255)
                                                                                             cGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCCTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaa                                       <  1:459740/141‑1 (MQ=255)
                                                                                               tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCCGgtttaagt                                     >  2:101851/1‑141 (MQ=255)
                                                                                               tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagt                                     >  2:446674/1‑141 (MQ=255)
                                                                                                 aatacacgccagaCCTTTTAAAGTATGCTGGGCTGTATTGGTACTTTGATCCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTAGAATGGTAGTCAGgtttaagttt                                   >  2:576185/13‑141 (MQ=14)
                                                                                                 aaTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTGTGATATGACGGTGAGCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttt                                   <  2:554976/141‑1 (MQ=255)
                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTGTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                   >  2:189656/1‑141 (MQ=255)
                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                   >  1:228990/1‑141 (MQ=255)
                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                   >  1:614731/1‑141 (MQ=255)
                                                                                                  aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta                                  >  1:64149/1‑141 (MQ=255)
                                                                                                   tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATATTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCTGgtttaggtttaa                                 >  2:332770/1‑141 (MQ=11)
                                                                                                   tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtttaa                                 >  1:175630/1‑141 (MQ=255)
                                                                                                   tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                 <  1:593294/141‑1 (MQ=255)
                                                                                                     cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATCCTTTGCTAAATTGCCTTTATTACCTTAATATGATGGGGATCTGGCGCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCGGGTTTAAGTTTaaaa                               >  1:132733/4‑141 (MQ=14)
                                                                                                      actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATTCCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCTATGTACAATGGCAGTCAGgttgtagtttacatg                              >  1:207421/3‑137 (MQ=14)
                                                                                                      actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGCTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCCGGTTTACGTTTAAAAg                              >  1:295053/3‑141 (MQ=255)
                                                                                                          aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCttt                          >  1:490856/1‑141 (MQ=255)
                                                                                                            ttccTTTTAAATTATTCTTTGCTGAATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTaa                        <  2:667934/139‑1 (MQ=255)
                                                                                                             aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                       >  1:536084/1‑141 (MQ=255)
                                                                                                             aCCTTTTAAAGTATTCTTTGCAGTTTTAAGGCTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                       >  2:305697/1‑141 (MQ=14)
                                                                                                                ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAc                    <  2:205563/141‑1 (MQ=255)
                                                                                                                  ttAAAGAATTCATTGCCGTATTGATCCTTTGATATCTTATCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTAttt                  >  2:255487/1‑141 (MQ=14)
                                                                                                                    aaaGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTa                >  2:45613/1‑141 (MQ=255)
                                                                                                                       gTATTCTTTGCTGTATTGCTGCTTTGATACCTTGTCTTTCTTAGTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTGAGATTTTTCGATGTACAATGGGAGTCAGGTTTAAGTTTGAACGCTTTAATTACTTTaatc             >  2:91227/1‑138 (MQ=11)
                                                                                                                        tATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAGAGCTTTAATTACTTTAAcct            >  1:547455/1‑140 (MQ=255)
                                                                                                                         aTTCTTTGCCGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAt           <  1:262590/141‑1 (MQ=255)
                                                                                                                          ttCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAtt          >  2:112129/1‑141 (MQ=255)
                                                                                                                            cTTTGCTGTATTGATACTTTGATACCTTGGCTTTCTTACTTTAATATGACGGTGATCTTGCTTAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGc        <  1:155151/141‑1 (MQ=255)
                                                                                                                               tGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTCGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAc     >  2:382368/1‑141 (MQ=255)
                                                                                                                                 cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCGTGCTCAATGAGGTTATTCAGATATTGCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCt   <  1:477974/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                             
CAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTT  >  CP000730/35867‑36137

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