Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,800,129 C→T 38.5% W71* (TGG→TAG)  ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,800,1290CT38.5% 303.9 / inf 296W71* (TGG→TAG) ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (99/83);  new base T (72/42);  total (171/125)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 299 total aligned reads displayed.

CAGCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCAT  >  CP009273/1799982‑1800276
                                                                                                                                                   |                                                                                                                                                   
     tGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCa                                                                                                                                               >  2:2903947/1‑149 (MQ=255)
      gCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAt                                                                                                                                              >  1:698324/1‑149 (MQ=255)
       ccTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtg                                                                                                                                             >  1:2983631/1‑149 (MQ=255)
         tGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggt                                                                                                                                           >  1:4729489/1‑149 (MQ=255)
         tGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCAtggt                                                                                                                                           >  1:147137/1‑149 (MQ=255)
         tGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCAtggt                                                                                                                                           <  2:2407906/149‑1 (MQ=255)
         tGACAAAATAAATCGTCAAAATCATTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAATAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCAtggt                                                                                                                                           <  1:1727277/149‑1 (MQ=255)
            ccAAATAAATAGTCATAATCAGTTAAATTAAAGCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATggtggt                                                                                                                                        <  2:2054199/149‑1 (MQ=255)
            ccAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATggtggt                                                                                                                                        <  1:2338589/149‑1 (MQ=255)
                ataaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCa                                                                                                                                    <  1:2713762/149‑1 (MQ=255)
                ataaataGTCATAATCAGTTAAATAAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCa                                                                                                                                    <  1:519966/149‑1 (MQ=255)
                    ataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCg                                                                                                                                >  1:1026065/1‑149 (MQ=255)
                    ataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCg                                                                                                                                >  2:3222343/1‑149 (MQ=255)
                          aTAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTg                                                                                                                          >  1:3116851/1‑149 (MQ=255)
                           tAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTga                                                                                                                         >  1:1106661/1‑149 (MQ=255)
                            aaTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgag                                                                                                                        <  2:3034354/149‑1 (MQ=255)
                             aTCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                       >  1:1085471/1‑149 (MQ=255)
                              tCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTgagga                                                                                                                      >  1:1959124/1‑149 (MQ=255)
                                aGTTAAATTAACACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGc                                                                                                                    <  1:4169560/149‑1 (MQ=255)
                                 gTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCt                                                                                                                   <  1:2728642/149‑1 (MQ=255)
                                 gTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCt                                                                                                                   >  1:3534168/1‑149 (MQ=255)
                                  ttaaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtt                                                                                                                  <  2:3258340/149‑1 (MQ=255)
                                    aaattcacACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACCGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtttt                                                                                                                <  1:829665/149‑1 (MQ=255)
                                      attaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCg                                                                                                              >  1:2722203/1‑149 (MQ=255)
                                      attaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCg                                                                                                              >  2:4049248/1‑149 (MQ=255)
                                       ttaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                             <  1:4667005/149‑1 (MQ=255)
                                         aaaaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTAt                                                                                                           >  1:3883650/1‑149 (MQ=255)
                                           aaCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTg                                                                                                         >  1:1918177/1‑149 (MQ=255)
                                            aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                        >  1:2838976/1‑149 (MQ=255)
                                            aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                        >  1:4519387/1‑149 (MQ=255)
                                            aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                        <  2:3571780/149‑1 (MQ=255)
                                            aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGc                                                                                                        >  1:2023981/1‑149 (MQ=255)
                                              ggCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCt                                                                                                      >  1:122258/1‑149 (MQ=255)
                                              ggCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCt                                                                                                      >  2:4079959/1‑149 (MQ=255)
                                              ggCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCt                                                                                                      >  2:2622203/1‑149 (MQ=255)
                                                   cgcccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGAtgg                                                                                                 <  2:899530/149‑3 (MQ=255)
                                                     cccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATtgtg                                                                                               >  1:2808475/1‑149 (MQ=255)
                                                     cccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATtgtg                                                                                               <  2:1654466/149‑1 (MQ=255)
                                                        gcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATGACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGATGAGCTTTTCGTATTGCCTGTTTGTGCaa                                                                                            <  2:1976168/149‑1 (MQ=255)
                                                          ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGa                                                                                          <  1:2604932/149‑1 (MQ=255)
                                                                 atatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTg                                                                                   >  2:1846181/1‑149 (MQ=255)
                                                                  tatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGc                                                                                  >  1:3006064/1‑149 (MQ=255)
                                                                       ttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCa                                                                             >  2:4571265/1‑149 (MQ=255)
                                                                        tCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAg                                                                            >  2:1756321/1‑149 (MQ=255)
                                                                         cTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGa                                                                           >  2:3688893/1‑149 (MQ=255)
                                                                          tGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGAt                                                                          >  2:4131887/1‑149 (MQ=255)
                                                                            aagaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTa                                                                        >  1:1259180/1‑149 (MQ=255)
                                                                             agaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAc                                                                       <  1:916150/149‑1 (MQ=255)
                                                                             agaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAc                                                                       >  2:2330034/1‑149 (MQ=255)
                                                                             agaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAc                                                                       <  2:3953980/149‑1 (MQ=255)
                                                                              gaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAcc                                                                      >  2:1475767/1‑149 (MQ=255)
                                                                               aaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCt                                                                     >  1:2456521/1‑149 (MQ=255)
                                                                               aaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCt                                                                     <  2:1890994/149‑1 (MQ=255)
                                                                                   cTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGa                                                                 >  1:2188888/1‑149 (MQ=255)
                                                                                   cTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGa                                                                 >  1:4094911/1‑149 (MQ=255)
                                                                                         atTAGAATAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                           >  1:612367/1‑149 (MQ=255)
                                                                                          ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                          >  1:2078482/1‑149 (MQ=255)
                                                                                            agaagaAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                        <  1:3741562/149‑1 (MQ=255)
                                                                                             gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtttttt                                                       <  1:3333101/149‑1 (MQ=255)
                                                                                                gaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgc                                                    >  1:3430347/1‑149 (MQ=255)
                                                                                                   gTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCt                                                 >  2:4470318/1‑149 (MQ=255)
                                                                                                    tATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                                >  1:4016351/1‑149 (MQ=255)
                                                                                                    tATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATCTCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                                >  2:3305090/1‑149 (MQ=255)
                                                                                                    tATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                                >  1:321816/1‑149 (MQ=255)
                                                                                                      ttACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                              >  2:42969/1‑149 (MQ=255)
                                                                                                      ttACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGc                                              <  2:4370980/149‑1 (MQ=255)
                                                                                                       tACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCa                                             <  2:4678159/149‑1 (MQ=255)
                                                                                                         cTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                           >  2:2746582/1‑149 (MQ=255)
                                                                                                         cTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                           >  2:4153644/1‑149 (MQ=255)
                                                                                                             cGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatca                                       >  2:3269554/1‑149 (MQ=255)
                                                                                                               ttCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCTGGatcatc                                     >  2:12798/1‑149 (MQ=255)
                                                                                                               ttCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatcatc                                     >  1:2131872/1‑149 (MQ=255)
                                                                                                               ttCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatcatc                                     <  1:3958880/149‑1 (MQ=255)
                                                                                                               ttCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatcatc                                     >  1:4303529/1‑149 (MQ=255)
                                                                                                               ttCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGatcatc                                     >  2:1153488/1‑149 (MQ=255)
                                                                                                                tCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATcc                                    <  1:1650442/149‑1 (MQ=255)
                                                                                                                 cccccACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCa                                   >  1:2285817/1‑149 (MQ=255)
                                                                                                                 cccccACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCa                                   <  2:2148529/149‑1 (MQ=255)
                                                                                                                 cccccACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCa                                   <  1:2920000/149‑1 (MQ=255)
                                                                                                                    ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGa                                <  2:4786539/149‑1 (MQ=255)
                                                                                                                    ccACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGa                                <  1:1553406/149‑1 (MQ=255)
                                                                                                                     cacaGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGAc                               <  2:1171797/149‑1 (MQ=255)
                                                                                                                      acaGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                              >  1:426172/1‑149 (MQ=255)
                                                                                                                        aGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGt                            >  1:4081230/1‑149 (MQ=255)
                                                                                                                        aGCGACCAGGCAGTGGTGTGTCCATACTAGCTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGt                            <  2:850905/149‑1 (MQ=255)
                                                                                                                         gCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTa                           >  2:2241787/1‑149 (MQ=255)
                                                                                                                         gCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTa                           >  2:730687/1‑149 (MQ=255)
                                                                                                                           gACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaa                         >  2:4270934/1‑149 (MQ=255)
                                                                                                                               aGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACg                     >  1:542907/1‑149 (MQ=255)
                                                                                                                               aGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACg                     <  2:1431956/149‑1 (MQ=255)
                                                                                                                                   aGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAat                 <  1:1365009/149‑1 (MQ=255)
                                                                                                                                       gtgtgtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcatc             <  2:2295510/149‑1 (MQ=255)
                                                                                                                                       gtgtgtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTAGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcatc             >  2:4530717/1‑149 (MQ=255)
                                                                                                                                        tgtgtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCg            >  1:4561306/1‑149 (MQ=255)
                                                                                                                                         gtgtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCgg           >  2:2646079/1‑149 (MQ=255)
                                                                                                                                          tgtCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGc          <  1:252414/149‑1 (MQ=255)
                                                                                                                                             ccATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGc       <  1:973403/149‑1 (MQ=255)
                                                                                                                                             ccATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGc       >  2:4383257/1‑149 (MQ=255)
                                                                                                                                               aTACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTg     >  1:1595624/1‑149 (MQ=255)
                                                                                                                                                tACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGc    >  2:2818075/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCAT  >  CP009273/1799982‑1800276

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: