Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,463,489 G→T 6.5% E294* (GAA→TAA)  gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,463,4890GT 6.5% 682.0 / 16.0 184E294* (GAA→TAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base G (98/74);  new base T (7/5);  total (105/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 185 total aligned reads displayed.

CTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAATCAAACAT  >  CP009273/2463342‑2463637
                                                                                                                                                   |                                                                                                                                                    
cTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGa                                                                                                                                                     <  2:583952/149‑1 (MQ=255)
 tttATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGaa                                                                                                                                                    >  1:3730401/1‑149 (MQ=255)
 tttATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGaa                                                                                                                                                    <  1:1955098/149‑1 (MQ=255)
  ttATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGaaa                                                                                                                                                   <  2:3556326/149‑1 (MQ=255)
  ttATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGaaa                                                                                                                                                   <  1:1798317/149‑1 (MQ=255)
   tataTGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAAc                                                                                                                                                  >  1:4049929/1‑149 (MQ=255)
    atatGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTTTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACg                                                                                                                                                 <  2:501340/149‑1 (MQ=255)
    atatGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACg                                                                                                                                                 >  1:1055806/1‑149 (MQ=255)
    atatGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACg                                                                                                                                                 >  1:1771812/1‑149 (MQ=255)
      atGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGtt                                                                                                                                               >  1:2151618/1‑149 (MQ=255)
       tGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGttt                                                                                                                                              <  2:2617013/149‑1 (MQ=255)
        gAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGtttt                                                                                                                                             <  2:4114267/149‑1 (MQ=255)
         aaaaaCAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTc                                                                                                                                            >  1:575285/1‑149 (MQ=255)
           aaaCAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGa                                                                                                                                          >  2:612935/1‑149 (MQ=255)
             acaacGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAAt                                                                                                                                        >  1:659192/1‑149 (MQ=255)
             acaacGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAAt                                                                                                                                        >  2:187494/1‑149 (MQ=255)
               aacGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGc                                                                                                                                      >  2:3542171/1‑149 (MQ=255)
                    gTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCActct                                                                                                                                 >  2:2730558/1‑149 (MQ=255)
                       taAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCtttt                                                                                                                              >  1:790768/1‑149 (MQ=255)
                         aaaTGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTa                                                                                                                            >  2:2673177/1‑149 (MQ=255)
                         aaaTGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTa                                                                                                                            >  2:533583/1‑149 (MQ=255)
                          aaTGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATTAAACGTTTTCGAATGCACTCTTTTTAt                                                                                                                           <  1:1349428/149‑1 (MQ=255)
                          aaTGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTAt                                                                                                                           >  1:2641436/1‑149 (MQ=255)
                          aaTGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTAt                                                                                                                           <  1:1983410/149‑1 (MQ=255)
                           aTGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTAtt                                                                                                                          >  1:494804/1‑149 (MQ=255)
                              ccATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGt                                                                                                                       <  2:3591632/149‑1 (MQ=255)
                               cATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCCTAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGtt                                                                                                                      <  1:3810199/149‑1 (MQ=255)
                               cATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGtt                                                                                                                      >  2:2881112/1‑149 (MQ=255)
                                   gTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCa                                                                                                                  >  1:2265247/1‑149 (MQ=255)
                                   gTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCa                                                                                                                  >  2:1742223/1‑149 (MQ=255)
                                   gTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCa                                                                                                                  >  2:3135930/1‑149 (MQ=255)
                                           atgatAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGc                                                                                                          >  2:1591894/1‑149 (MQ=255)
                                           atgatAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGc                                                                                                          >  1:2240521/1‑149 (MQ=255)
                                              atAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGACCGGAATGTTTCGAATCTCATAAAGATGACACGTGTTCGAATGCACTCTTTGTATTGGTTAGCAAACCAAGCAcc                                                                                                       >  2:2542417/1‑149 (MQ=255)
                                              atAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCAcc                                                                                                       >  2:2542206/1‑149 (MQ=255)
                                              atAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCAcc                                                                                                       >  2:3818424/1‑149 (MQ=255)
                                               tAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCa                                                                                                      >  2:3344601/1‑149 (MQ=255)
                                                aattgtgttgtgtTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCAt                                                                                                     <  2:3092246/146‑1 (MQ=255)
                                                 aGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATc                                                                                                    >  1:1556342/1‑149 (MQ=255)
                                                  gtgtgtTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCt                                                                                                   <  1:1091547/149‑1 (MQ=255)
                                                       ttGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCaaa                                                                                              <  2:3212300/149‑1 (MQ=255)
                                                        tGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAAc                                                                                             <  2:1543768/149‑1 (MQ=255)
                                                          ggTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACtt                                                                                           >  2:3772200/1‑149 (MQ=255)
                                                            ttAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTcc                                                                                         <  1:1127730/149‑1 (MQ=255)
                                                            ttAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTcc                                                                                         >  2:1043412/1‑149 (MQ=255)
                                                            ttAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTcc                                                                                         <  1:2393931/149‑1 (MQ=255)
                                                               gATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGAAGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCAtt                                                                                      >  2:737521/1‑149 (MQ=255)
                                                                 tGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTtg                                                                                    <  2:2710775/149‑1 (MQ=255)
                                                                  gCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTtga                                                                                   <  2:518639/149‑1 (MQ=255)
                                                                  gCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTtga                                                                                   <  1:1422573/149‑1 (MQ=255)
                                                                   ccTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTtgat                                                                                  >  1:3520383/1‑149 (MQ=255)
                                                                         gTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGgt                                                                            >  2:1525049/1‑149 (MQ=255)
                                                                         gTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGgt                                                                            >  2:1283836/1‑149 (MQ=255)
                                                                          taataaATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGgtg                                                                           >  2:1079739/1‑149 (MQ=255)
                                                                             taaATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGa                                                                        >  2:749569/1‑149 (MQ=255)
                                                                              aaaTTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGAt                                                                       >  1:3197580/1‑149 (MQ=255)
                                                                              aaaTTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGAt                                                                       >  1:3196651/1‑149 (MQ=255)
                                                                               aaTTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATg                                                                      >  2:955471/1‑149 (MQ=255)
                                                                                aTTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGt                                                                     <  2:1055806/149‑1 (MQ=255)
                                                                                   cGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGtctc                                                                  <  2:2462963/149‑1 (MQ=255)
                                                                                    gACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCa                                                                 <  1:444444/149‑1 (MQ=255)
                                                                                      cATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGt                                                               <  1:1980025/149‑1 (MQ=255)
                                                                                       atatCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGta                                                              >  2:1878352/1‑149 (MQ=255)
                                                                                           cATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATaaa                                                          >  2:89979/1‑149 (MQ=255)
                                                                                            aTTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAg                                                         >  2:3797785/1‑149 (MQ=255)
                                                                                             tttGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGa                                                        >  2:1488322/1‑149 (MQ=255)
                                                                                               tGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGaaa                                                      <  1:2700185/149‑1 (MQ=255)
                                                                                                 gcgcAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAat                                                    >  2:893526/1‑149 (MQ=255)
                                                                                                 gcgcAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAat                                                    >  2:893551/1‑149 (MQ=255)
                                                                                                   gcAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAatta                                                  >  1:1436520/1‑149 (MQ=255)
                                                                                                     aaCCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAattatt                                                <  2:1029795/149‑1 (MQ=255)
                                                                                                       ccccAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTc                                              <  2:346460/149‑1 (MQ=255)
                                                                                                       ccccAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTc                                              <  2:1975643/149‑1 (MQ=255)
                                                                                                         ccAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCa                                            >  1:2260219/1‑149 (MQ=255)
                                                                                                            aCAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGt                                         >  2:2149050/1‑149 (MQ=255)
                                                                                                             cAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGta                                        <  2:2106976/149‑1 (MQ=255)
                                                                                                              aGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGtat                                       <  1:3438207/149‑1 (MQ=255)
                                                                                                               gAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGtata                                      >  2:165985/1‑149 (MQ=255)
                                                                                                                   ttGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATaaa                                  >  1:1003370/1‑149 (MQ=255)
                                                                                                                      gAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATaaaaaa                               <  1:1023982/149‑1 (MQ=255)
                                                                                                                         cGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTa                            >  1:2962075/1‑149 (MQ=255)
                                                                                                                          ggAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTac                           >  2:851762/1‑149 (MQ=255)
                                                                                                                          ggAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTac                           <  1:3818424/149‑1 (MQ=255)
                                                                                                                          ggAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTac                           >  1:3364741/1‑149 (MQ=255)
                                                                                                                           gAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTaca                          >  2:362250/1‑149 (MQ=255)
                                                                                                                           gAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTaca                          <  2:1901543/149‑1 (MQ=255)
                                                                                                                            aaTGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTacac                         <  2:451740/149‑1 (MQ=255)
                                                                                                                            aaTGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTacac                         >  2:718200/1‑149 (MQ=255)
                                                                                                                            aaTGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTacac                         <  2:4050190/149‑1 (MQ=255)
                                                                                                                                tttCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTaa                     <  2:2641436/149‑1 (MQ=255)
                                                                                                                                 ttCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAat                    >  2:3010128/1‑149 (MQ=255)
                                                                                                                                  tCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAata                   >  2:3986295/1‑149 (MQ=255)
                                                                                                                                   cGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAatat                  <  1:2730686/149‑1 (MQ=255)
                                                                                                                                    gAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAatata                 <  1:1742223/149‑1 (MQ=255)
                                                                                                                                     aaTCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAatatat                >  1:22052/1‑149 (MQ=255)
                                                                                                                                      aTCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATg               >  1:2903808/1‑149 (MQ=255)
                                                                                                                                            tAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAAt         >  1:2054040/1‑149 (MQ=255)
                                                                                                                                            tAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAAt         >  1:2053175/1‑149 (MQ=255)
                                                                                                                                            tAAAGATTAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAAt         >  2:805101/1‑149 (MQ=255)
                                                                                                                                            tAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAAt         <  2:2558033/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTTATATGAAAAACAACGGGTATAAAATGCCATCGTATGTTGATGATAAGTGTGTTGGGTTAGATGCCTGGGGTAATAAATTCGACATATCATTTGGCGCAACCCCAACAGAAGTTGGAACGGAATGTTTCGAATCTCATAAAGATGAAACGTTTTCGAATGCACTCTTTTTATTGGTTAGCAAACCAAGCACCATCTTCAAACTTCCATTTGATGATGGTGTGATGTCTCAGTATAAAGAAAATTATTTCCATGTATATAAAAAACTACACGTAATATATGGAGAATCAAACAT  >  CP009273/2463342‑2463637

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: