Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,090,918 T→A 7.7% 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,9180TA 7.7% 682.2 / 38.6 196intergenic (‑25/‑43)insE1/ymdEIS3 transposase A/pseudogene
Reads supporting (aligned to +/- strand):  ref base T (93/88);  new base A (10/5);  total (103/93)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 198 total aligned reads displayed.

CAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACA  >  CP009273/1090812‑1091065
                                                                                                          |                                                                                                                                                   
cAGCAAGCTTCAGGTCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGt                                                                                                           <  1:606177/149‑1 (MQ=11)
cAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGt                                                                                                           <  1:3229542/149‑1 (MQ=11)
   ctaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGTGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                        <  1:2769925/147‑1 (MQ=17)
   cAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGcc                                                                                                        <  1:1271097/149‑1 (MQ=32)
    aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCt                                                                                                       >  1:3613289/1‑149 (MQ=33)
    aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTAAGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTTAGAGTTTACTCACTTAGCCGCGTGTCCACTATTTCTGGGTAAGATCAGTGCCt                                                                                                       <  1:1666802/149‑1 (MQ=16)
     aGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                      >  2:2636776/1‑149 (MQ=34)
     aGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                      >  2:2902105/1‑149 (MQ=34)
     aGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                      <  2:1152413/149‑1 (MQ=34)
      gCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTgg                                                                                                     >  1:1953881/1‑149 (MQ=35)
      gCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTgg                                                                                                     <  1:2456564/149‑1 (MQ=35)
        ttCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGca                                                                                                   >  1:686575/1‑149 (MQ=37)
          cAGGGCGTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGcaca                                                                                                 <  2:3616591/149‑1 (MQ=38)
             ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                              <  1:647226/149‑1 (MQ=40)
             ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                              <  1:3825329/149‑1 (MQ=40)
              gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgcg                                                                                             <  2:3312333/149‑1 (MQ=255)
              gCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgcg                                                                                             >  1:4035342/1‑149 (MQ=255)
               cTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGt                                                                                            <  2:3861411/149‑1 (MQ=255)
                ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                           <  2:4790219/149‑1 (MQ=255)
                ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                           <  2:4933426/149‑1 (MQ=255)
                  cACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTg                                                                                         <  1:671367/149‑1 (MQ=255)
                     tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATg                                                                                      <  2:1050058/149‑1 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  <  1:3289649/149‑1 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  >  2:1191768/1‑149 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  <  2:1782602/149‑1 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  >  2:2606435/1‑149 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  <  2:3165394/149‑1 (MQ=255)
                         aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                  >  2:530187/1‑149 (MQ=255)
                          aaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGa                                                                                 <  2:1452910/149‑1 (MQ=255)
                          aaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGa                                                                                 <  2:2226002/149‑1 (MQ=255)
                          aaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGa                                                                                 <  2:645780/149‑1 (MQ=255)
                               aGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCg                                                                            >  2:533316/1‑149 (MQ=255)
                               aGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCg                                                                            >  2:786983/1‑149 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGt                                                                           <  2:4234469/149‑1 (MQ=255)
                                ggCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATAAAGTTTTTTTCATGTGATTCACCTATGACTGAGAGTTTACTCACATAGACTCGTGTCCACTATTTCTGGGTAAGATCAGTTCCTGGCACACTCGTTTGATGGACGACTGCGt                                                                           >  1:1785779/1‑149 (MQ=14)
                                  cGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTc                                                                         <  2:3588716/149‑1 (MQ=255)
                                   gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCa                                                                        >  1:2039785/1‑149 (MQ=255)
                                    aaTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAg                                                                       <  1:2676657/149‑1 (MQ=255)
                                      tGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGc                                                                     <  1:1683382/149‑1 (MQ=255)
                                       gCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCa                                                                    >  1:3069784/1‑149 (MQ=255)
                                         tGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                  >  2:1712307/1‑149 (MQ=255)
                                          gTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATg                                                                 <  2:4661719/149‑1 (MQ=255)
                                           tttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgc                                                                >  2:2250694/1‑149 (MQ=255)
                                             tACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgcgc                                                              >  1:4432411/1‑149 (MQ=255)
                                                ggggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAt                                                           >  1:391339/1‑149 (MQ=255)
                                                 gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAtt                                                          >  2:4092896/1‑149 (MQ=255)
                                                 gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAtt                                                          >  2:1475276/1‑149 (MQ=255)
                                                   gTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTAt                                                        <  2:2492561/149‑1 (MQ=255)
                                                   gTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTAt                                                        <  1:2939129/149‑1 (MQ=255)
                                                   gTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTAt                                                        >  2:4907475/1‑149 (MQ=255)
                                                   gTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCGGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTAt                                                        <  2:4864965/149‑1 (MQ=255)
                                                    ttttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                       >  1:855657/1‑149 (MQ=255)
                                                         tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGcc                                                  >  2:988050/1‑149 (MQ=255)
                                                         tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGcc                                                  >  2:1076492/1‑149 (MQ=255)
                                                           cTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCgg                                                >  2:2738643/1‑149 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                              <  1:4684350/149‑1 (MQ=255)
                                                             ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGAc                                              >  2:2899022/1‑149 (MQ=255)
                                                              gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                             >  2:3142934/1‑149 (MQ=255)
                                                              gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                             >  2:4424066/1‑149 (MQ=255)
                                                                tGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTgg                                           >  1:644727/1‑149 (MQ=255)
                                                                 gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGa                                          >  2:2885296/1‑149 (MQ=255)
                                                                  aTACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGat                                         >  1:3038712/1‑149 (MQ=255)
                                                                    aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                       >  1:2531621/1‑149 (MQ=255)
                                                                     cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATc                                      <  2:151389/149‑1 (MQ=255)
                                                                     cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATc                                      <  2:866988/149‑1 (MQ=255)
                                                                     cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATc                                      <  2:373806/149‑1 (MQ=255)
                                                                       gTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                    >  1:2155642/1‑149 (MQ=255)
                                                                        tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCt                                   <  2:3015002/149‑1 (MQ=255)
                                                                        tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCt                                   >  2:1963265/1‑149 (MQ=255)
                                                                             gTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGcc                              >  2:4434189/1‑149 (MQ=255)
                                                                              tCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCa                             >  2:442201/1‑149 (MQ=255)
                                                                               cATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCag                            >  1:761997/1‑149 (MQ=255)
                                                                                aTGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCaga                           <  2:3476191/149‑1 (MQ=255)
                                                                                 tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCagag                          >  1:4361292/1‑149 (MQ=255)
                                                                                 tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCagag                          >  1:4361611/1‑149 (MQ=255)
                                                                                 tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCagag                          >  2:4537753/1‑149 (MQ=255)
                                                                                 tgtgAGTCACCTCTGACTGAGAGTTAACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCagag                          >  2:2084678/1‑149 (MQ=255)
                                                                                   tgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCa                        <  1:3952680/149‑1 (MQ=255)
                                                                                   tgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCa                        >  1:526245/1‑149 (MQ=255)
                                                                                     aGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATc                      <  2:4006160/149‑1 (MQ=255)
                                                                                      gTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCa                     <  1:713468/149‑1 (MQ=255)
                                                                                        cACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                   <  2:4501838/149‑1 (MQ=255)
                                                                                           ctctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTc                <  2:1274762/149‑1 (MQ=255)
                                                                                           ctctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTc                <  1:3772345/149‑1 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt               <  1:3328815/149‑1 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt               >  1:4597996/1‑149 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt               >  1:1458356/1‑149 (MQ=255)
                                                                                            tctGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCt               >  2:1307810/1‑149 (MQ=255)
                                                                                                actgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg           <  1:492720/149‑1 (MQ=255)
                                                                                                 ctgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatga          >  1:1991548/1‑149 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa         >  2:2698611/1‑149 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa         <  1:1038914/149‑1 (MQ=255)
                                                                                                  tgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa         >  2:2864604/1‑149 (MQ=255)
                                                                                                    agagTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACg       >  1:4984000/1‑149 (MQ=255)
                                                                                                     gagTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGt      >  1:4147987/1‑149 (MQ=255)
                                                                                                      agTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTa     <  2:4644445/149‑1 (MQ=255)
                                                                                                      agTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTa     <  2:2433102/149‑1 (MQ=255)
                                                                                                       gTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTaa    <  1:2871100/149‑1 (MQ=255)
                                                                                                         ttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa  <  1:840779/149‑1 (MQ=255)
                                                                                                         ttACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa  <  1:5164670/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: