Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,090,918 T→A 11.4% intergenic (‑25/‑43) insE1 ← / → ymdE IS3 transposase A/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,090,9180TA11.4% 562.9 / 54.7 176intergenic (‑25/‑43)insE1/ymdEIS3 transposase A/pseudogene
Reads supporting (aligned to +/- strand):  ref base T (97/59);  new base A (11/9);  total (108/68)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 182 total aligned reads displayed.

CAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACA  >  CP009273/1090812‑1091065
                                                                                                          |                                                                                                                                                   
   cAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                        <  2:3061632/149‑1 (MQ=32)
   cAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                        >  1:4183782/1‑149 (MQ=32)
   cAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                        <  2:3757555/149‑1 (MQ=32)
    aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCt                                                                                                       >  2:892888/1‑149 (MQ=33)
     aGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                      >  2:2346610/1‑149 (MQ=34)
     aGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGACGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                      >  2:2397980/1‑149 (MQ=21)
      gCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTgg                                                                                                     >  2:2536913/1‑149 (MQ=35)
        ttCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGca                                                                                                   >  1:1152247/1‑149 (MQ=37)
          cAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGcaca                                                                                                 <  1:1640942/149‑1 (MQ=38)
          cAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGcaca                                                                                                 <  1:2101007/149‑1 (MQ=38)
           aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCacac                                                                                                >  2:728524/1‑149 (MQ=39)
            gggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcg                                                                                               >  2:1371545/1‑149 (MQ=40)
             ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTATGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                              <  1:2847244/149‑1 (MQ=255)
              gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgcg                                                                                             >  1:4558158/1‑149 (MQ=255)
              gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgcg                                                                                             >  1:4558192/1‑149 (MQ=255)
              gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgcg                                                                                             <  2:31625/149‑1 (MQ=255)
               cTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGt                                                                                            <  1:4540785/149‑1 (MQ=255)
                ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                           <  1:1720473/149‑1 (MQ=255)
                 tCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGttt                                                                                          >  2:3398477/1‑149 (MQ=255)
                 tCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGttt                                                                                          >  2:833242/1‑149 (MQ=255)
                    cTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGAt                                                                                       <  1:3454932/149‑1 (MQ=255)
                     tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATg                                                                                      >  1:1718305/1‑149 (MQ=255)
                        gAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTTATGGaa                                                                                   >  2:4419837/1‑149 (MQ=255)
                        gAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGaa                                                                                   >  2:455329/1‑149 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  >  2:400547/1‑149 (MQ=255)
                           aTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGAc                                                                                <  1:1498662/149‑1 (MQ=255)
                            ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACg                                                                               >  1:2533483/1‑149 (MQ=255)
                            ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACg                                                                               <  2:3514662/149‑1 (MQ=255)
                             tCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTCAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACgg                                                                              >  1:3527079/1‑149 (MQ=255)
                             tCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACgg                                                                              >  1:4437942/1‑149 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                           <  2:979559/149‑1 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                           >  2:3230712/1‑149 (MQ=255)
                                   gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCa                                                                        <  1:2185541/149‑1 (MQ=255)
                                   gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCa                                                                        >  1:2273197/1‑149 (MQ=255)
                                   gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCa                                                                        <  1:4300613/149‑1 (MQ=255)
                                    aaTGCTGTTTACGGGGTTTTTTACTGGTTGATACAGTTTTTGTCATCTGAGTAACCTCTGACTGAGAGTTTACTCACTAAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGACTGGCACACGCGATTGATGGACGACGGCGGTCAg                                                                       >  2:2865561/1‑149 (MQ=17)
                                     aTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAgg                                                                      >  2:2695775/1‑149 (MQ=255)
                                      tGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGc                                                                     >  1:2778983/1‑149 (MQ=255)
                                         tGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                  >  2:2490789/1‑149 (MQ=255)
                                          gTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATg                                                                 >  1:4062005/1‑149 (MQ=255)
                                           tttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgc                                                                >  2:4093099/1‑149 (MQ=255)
                                           tttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgc                                                                >  2:827949/1‑149 (MQ=255)
                                           tttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgc                                                                <  2:3124630/149‑1 (MQ=255)
                                              aCGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGcgcg                                                             <  2:1963558/149‑1 (MQ=255)
                                              aCGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGcgcg                                                             >  1:380150/1‑149 (MQ=255)
                                                ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAt                                                           >  2:2518126/1‑149 (MQ=255)
                                                ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAt                                                           >  2:1680304/1‑149 (MQ=255)
                                                ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAt                                                           >  2:1245903/1‑149 (MQ=255)
                                                ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAt                                                           >  2:1245995/1‑149 (MQ=255)
                                                 gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAtt                                                          >  1:4044054/1‑149 (MQ=255)
                                                 gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAtt                                                          <  2:1531217/149‑1 (MQ=255)
                                                  ggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTa                                                         >  2:2181759/1‑149 (MQ=255)
                                                    ttttttACTGGTTGATACAGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                       <  1:3339251/149‑1 (MQ=255)
                                                      ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCtt                                                     >  2:261980/1‑149 (MQ=255)
                                                       tttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTg                                                    <  1:1411695/149‑1 (MQ=255)
                                                          aCTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCg                                                 >  1:1223460/1‑149 (MQ=255)
                                                            tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGa                                               <  2:1226724/149‑1 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                              <  1:3560247/149‑1 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                              <  1:3207729/149‑1 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                              <  1:3689353/149‑1 (MQ=255)
                                                              gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTTACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                             >  2:4571022/1‑149 (MQ=255)
                                                              gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                             <  2:1464717/149‑1 (MQ=255)
                                                                 gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGa                                          >  1:2319128/1‑149 (MQ=255)
                                                                 gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGa                                          >  2:574294/1‑149 (MQ=255)
                                                                  aTACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGat                                         >  1:1455500/1‑149 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                       >  1:700275/1‑149 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                       >  1:4033862/1‑149 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                       <  1:2010593/149‑1 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                       >  2:4178423/1‑149 (MQ=255)
                                                                       gttgttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    <  2:2379749/149‑1 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    >  1:10619/1‑149 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    >  1:10723/1‑149 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    >  1:2568767/1‑149 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    >  1:3780500/1‑149 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    <  2:1643361/149‑1 (MQ=255)
                                                                       gTTTGTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    <  2:4456948/149‑1 (MQ=255)
                                                                             gTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGcc                              >  1:1746775/1‑149 (MQ=255)
                                                                               cATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCag                            >  2:3826983/1‑149 (MQ=255)
                                                                                        cACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                   >  2:3101940/1‑149 (MQ=255)
                                                                                         aCCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGt                  <  2:3307612/149‑1 (MQ=255)
                                                                                          ccTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                 <  1:4067315/149‑1 (MQ=255)
                                                                                          ccTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                 >  1:109274/1‑149 (MQ=255)
                                                                                           ctctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTc                >  1:2981975/1‑149 (MQ=255)
                                                                                               gactgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgat            >  1:2490061/1‑149 (MQ=255)
                                                                                               gactgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgat            >  2:4021566/1‑149 (MQ=255)
                                                                                               gactgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgat            >  1:1365221/1‑149 (MQ=255)
                                                                                                actgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg           >  2:721144/1‑149 (MQ=255)
                                                                                                actgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg           >  1:3848035/1‑149 (MQ=255)
                                                                                                 ctgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatga          >  1:1828895/1‑149 (MQ=255)
                                                                                                 ctgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatga          >  2:1752190/1‑149 (MQ=255)
                                                                                                   gagaGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAAc        >  2:3464004/1‑149 (MQ=255)
                                                                                                   gagaGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAAc        >  2:1757584/1‑149 (MQ=255)
                                                                                                      agTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTa     >  1:977151/1‑149 (MQ=255)
                                                                                                       gTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTaa    <  1:872353/149‑1 (MQ=255)
                                                                                                       gTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTaa    <  1:3072464/149‑1 (MQ=255)
                                                                                                        tttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAAc   <  1:3016981/149‑1 (MQ=255)
                                                                                                         ttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa  <  2:4368438/149‑1 (MQ=255)
                                                                                                         ttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa  <  2:2336574/149‑1 (MQ=255)
                                                                                                         ttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa  <  2:4131226/149‑1 (MQ=255)
                                                                                                         ttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa  <  2:3107547/149‑1 (MQ=255)
                                                                                                          |                                                                                                                                                   
CAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACA  >  CP009273/1090812‑1091065

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