Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,090,918 T→A 11.9% 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.9% 760.4 / 81.8 236intergenic (‑25/‑43)insE1/ymdEIS3 transposase A/pseudogene
Reads supporting (aligned to +/- strand):  ref base T (118/90);  new base A (21/7);  total (139/97)
Fisher's exact test for biased strand distribution p-value = 6.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 239 total aligned reads displayed.

CAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAAC  >  CP009273/1090812‑1091064
                                                                                                          |                                                                                                                                                  
cAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGt                                                                                                          >  1:851249/1‑149 (MQ=11)
  gCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGc                                                                                                        >  2:5033557/1‑149 (MQ=32)
   cAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                       >  2:2389893/1‑149 (MQ=32)
   cAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                       <  1:3715171/149‑1 (MQ=32)
    aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCt                                                                                                      <  1:1091903/149‑1 (MQ=33)
    aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCt                                                                                                      >  2:1511521/1‑149 (MQ=33)
     aGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                     >  1:3920767/1‑149 (MQ=34)
      gCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTgg                                                                                                    >  2:2327369/1‑149 (MQ=35)
       attCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGATTTTACTCACTTATCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGc                                                                                                   >  1:50338/2‑149 (MQ=255)
            gggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcg                                                                                              >  2:2453000/1‑149 (MQ=40)
             ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                             <  1:1690222/149‑1 (MQ=40)
             ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                             <  1:4348493/149‑1 (MQ=40)
              gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgcg                                                                                            >  1:1285200/1‑149 (MQ=255)
                ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                          <  1:4284906/149‑1 (MQ=255)
                ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                          <  2:4283625/149‑1 (MQ=255)
                  cACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTg                                                                                        >  1:1220296/1‑149 (MQ=255)
                   aCTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGa                                                                                       >  1:4975450/1‑149 (MQ=255)
                    cTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGAt                                                                                      >  1:1571058/1‑149 (MQ=255)
                     tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATg                                                                                     >  1:2063733/1‑149 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                 <  1:4152878/149‑1 (MQ=255)
                          aaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGa                                                                                <  2:3205151/149‑1 (MQ=255)
                            ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACg                                                                              >  1:226398/1‑149 (MQ=255)
                            ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACg                                                                              >  1:4603019/1‑149 (MQ=255)
                             tCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACgg                                                                             >  2:3327956/1‑149 (MQ=255)
                              cAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGc                                                                            <  2:466242/149‑1 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                          <  1:767239/149‑1 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                          >  1:4220021/1‑149 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                          >  2:4307017/1‑149 (MQ=255)
                                 gCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGtt                                                                         <  1:920695/149‑1 (MQ=255)
                                    aaTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAg                                                                      <  1:1229532/149‑1 (MQ=255)
                                     aTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAgg                                                                     >  1:3599694/1‑149 (MQ=255)
                                        cTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCaa                                                                  >  1:2084747/1‑149 (MQ=255)
                                        cTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCaa                                                                  <  1:1762705/149‑1 (MQ=255)
                                         tGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                 >  2:2173993/1‑149 (MQ=255)
                                         tGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                 <  2:3634567/149‑1 (MQ=255)
                                         tGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                 <  2:726528/149‑1 (MQ=255)
                                          gTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATg                                                                >  1:4829145/1‑149 (MQ=255)
                                                ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAt                                                          >  1:2786521/1‑149 (MQ=255)
                                                  ggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTa                                                        >  2:1756964/1‑149 (MQ=255)
                                                    ttttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                      >  1:1167166/1‑149 (MQ=255)
                                                    ttttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTAAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                      >  2:1093837/1‑149 (MQ=255)
                                                      ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCtt                                                    >  2:2804728/1‑149 (MQ=255)
                                                        ttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGc                                                  >  1:295562/1‑149 (MQ=255)
                                                            tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGa                                              >  2:1400830/1‑149 (MQ=255)
                                                            tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGa                                              >  1:2166079/1‑149 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                             <  2:668323/149‑1 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                             <  2:2063733/149‑1 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                             <  1:5173689/149‑1 (MQ=255)
                                                                 gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGa                                         <  2:763254/149‑1 (MQ=255)
                                                                 gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGa                                         >  1:3509361/1‑149 (MQ=255)
                                                                  aTACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGat                                        >  1:668644/1‑149 (MQ=255)
                                                                  aTACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGat                                        >  2:1849162/1‑149 (MQ=255)
                                                                   tACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGata                                       <  1:2645406/149‑1 (MQ=255)
                                                                   tACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGata                                       <  1:3311969/149‑1 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                      >  2:2935500/1‑149 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                      <  2:1494806/149‑1 (MQ=255)
                                                                     cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATc                                     <  2:844587/149‑1 (MQ=255)
                                                                     cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATc                                     <  2:4216818/149‑1 (MQ=255)
                                                                     cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATc                                     >  2:1538872/1‑149 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                   >  2:2532364/1‑149 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                   >  2:4146730/1‑149 (MQ=255)
                                                                        tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCt                                  <  2:318193/149‑1 (MQ=255)
                                                                        tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCt                                  >  1:1613174/1‑149 (MQ=255)
                                                                        tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCt                                  >  2:281609/1‑149 (MQ=255)
                                                                         ttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCtt                                 >  2:1147406/1‑149 (MQ=255)
                                                                          tttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTc                                >  2:150313/1‑149 (MQ=255)
                                                                          tttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTc                                >  2:2519940/1‑149 (MQ=255)
                                                                               cATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCag                           <  1:2389893/149‑1 (MQ=255)
                                                                                aTGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCaga                          >  1:3157441/1‑149 (MQ=255)
                                                                                 tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCagag                         >  1:4908148/1‑149 (MQ=255)
                                                                                 tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCagag                         >  1:2175847/1‑149 (MQ=255)
                                                                                    gAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCTATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCAt                      <  1:153453/149‑1 (MQ=255)
                                                                                      gacACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCa                    >  1:2497241/3‑149 (MQ=255)
                                                                                      gTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCa                    <  2:128444/149‑1 (MQ=255)
                                                                                      gTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCa                    >  1:2497210/1‑149 (MQ=255)
                                                                                       tCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAg                   <  1:1056473/149‑1 (MQ=255)
                                                                                       tCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAg                   >  2:1387165/1‑149 (MQ=255)
                                                                                       tCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAg                   >  2:3557595/1‑149 (MQ=255)
                                                                                        cACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                  <  1:3579062/149‑1 (MQ=255)
                                                                                         aCCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGt                 >  2:487416/1‑149 (MQ=255)
                                                                                         aCCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGt                 >  2:3831071/1‑149 (MQ=255)
                                                                                          ccTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                <  1:2071981/149‑1 (MQ=255)
                                                                                          ccTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                >  1:4458298/1‑149 (MQ=255)
                                                                                           ctctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTc               <  2:3493078/149‑1 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCACTGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt              >  2:1183806/1‑149 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt              >  1:551912/1‑149 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt              >  1:4434938/1‑149 (MQ=255)
                                                                                             ctgactgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtg             >  1:2525288/1‑149 (MQ=255)
                                                                                              tgactgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtga            <  1:4661304/149‑1 (MQ=255)
                                                                                                actgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg          <  1:2424881/149‑1 (MQ=255)
                                                                                                actgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg          <  1:3423127/149‑1 (MQ=255)
                                                                                                 ctgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatga         >  1:4557575/1‑149 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa        <  2:4027360/149‑1 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa        <  2:10179/149‑1 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa        >  2:2206189/1‑149 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa        <  1:1322626/149‑1 (MQ=255)
                                                                                                   gagaGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAAc       >  2:1402685/1‑149 (MQ=255)
                                                                                                   gagaGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAAc       <  1:1451249/149‑1 (MQ=255)
                                                                                                     gagTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGt     <  2:4982044/149‑1 (MQ=255)
                                                                                                        tttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAAc  <  2:1674211/149‑1 (MQ=255)
                                                                                                          |                                                                                                                                                  
CAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAAC  >  CP009273/1090812‑1091064

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: