Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,463,906 C→A 19.9% Q433K (CAA→AAA)  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,9060CA19.9% 255.1 / 26.5 114Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (18/70);  new base A (17/9);  total (35/79)
Fisher's exact test for biased strand distribution p-value = 4.66e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  CP009273/2463758‑2464054
                                                                                                                                                    |                                                                                                                                                      
tGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTc                                                                                                                                                        <  1:1271539/149‑1 (MQ=255)
  tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaa                                                                                                                                                      <  2:1963392/149‑1 (MQ=255)
    gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTTTAAAAAAAATaaaaa                                                                                                                                                    >  1:904654/1‑148 (MQ=255)
    gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAAAAAATaaaat                                                                                                                                                    >  1:4245916/1‑148 (MQ=255)
             tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                           <  2:1388425/149‑1 (MQ=255)
             tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                           <  2:1388407/149‑1 (MQ=255)
              ttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTTTTTTtaaaaaaaataaaaaaaaaaaaaaa                                                                                                                                          >  1:875172/1‑142 (MQ=255)
                ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACaa                                                                                                                                        <  1:3162774/149‑1 (MQ=255)
                 gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                                                                                                       <  1:3911587/149‑1 (MQ=255)
                 gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAAAAATaaaaaaaaa                                                                                                                                       >  2:1605367/1‑143 (MQ=255)
                     atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                                                                                                   <  2:1505078/149‑1 (MQ=255)
                          aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                              <  2:217275/149‑1 (MQ=255)
                             tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATTATTAAAATAAATAAAGACAATAGAGATGttatt                                                                                                                           >  2:2719369/1‑145 (MQ=255)
                             tttATGTGTCATTTTTCGGGGAAGGAGATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                           <  2:3335100/149‑1 (MQ=255)
                                  gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAATAAAAAATCGAAGTAAAacacaa                                                                                                                      >  2:4227648/1‑147 (MQ=255)
                                   tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                     <  1:4231377/149‑1 (MQ=255)
                                   tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                     <  1:4231307/149‑1 (MQ=255)
                                   tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAAAAATAAAAAAAAAAAAAAAAAAAAACGATGTTAAGAACaaa                                                                                                                     >  1:530838/1‑149 (MQ=255)
                                       aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                <  2:3966790/149‑1 (MQ=255)
                                       aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                <  2:3965956/149‑1 (MQ=255)
                                          tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  1:2057125/149‑1 (MQ=255)
                                            tCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                                                                                           <  1:1919188/149‑1 (MQ=255)
                                             cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                                                                                                          <  1:1894695/149‑1 (MQ=255)
                                              ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                                                                                                         <  1:3468926/149‑1 (MQ=255)
                                               gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAAAAAACACAAAATCGATTTAAATCACAAAATTTaaaattag                                                                                                        >  2:2966177/1‑144 (MQ=255)
                                               gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAAAAAAAAATCGATGTAAATCAAAAAATTTAAAtttag                                                                                                        >  2:2870319/1‑146 (MQ=255)
                                                ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                                                                                                       <  1:1010683/149‑1 (MQ=255)
                                                   aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                                                                                                   <  2:1983375/149‑1 (MQ=255)
                                                       tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                                                                                               <  2:1838378/149‑1 (MQ=255)
                                                       tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                                                                                               <  2:3720552/149‑1 (MQ=255)
                                                            agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                          <  2:1726538/149‑1 (MQ=255)
                                                             gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg                                                                                         <  2:2033838/149‑1 (MQ=255)
                                                              aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                        <  2:3186789/149‑1 (MQ=255)
                                                                    gcaagcaTTTATTTGGAATTTATTTTTCGTTCGACCTTTGCTTATAAATAACTGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATAATGGCGCGCACaa                                                                                  <  1:1532253/149‑1 (MQ=255)
                                                                    gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTAAAAAAAAAATCAAACAAAACAACACAAAAACCATGTAAAGCACGAAGTTTAAATTTCGCGCCACGCAAGGCGCGcccac                                                                                  >  2:3181141/1‑148 (MQ=255)
                                                                    gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                  <  2:3505855/149‑1 (MQ=255)
                                                                    gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAAAAAGAAAAAATAGATGTAAATAACAAACTTTAAATTGAGCGACAAACATGTAGCGCACaa                                                                                  >  2:3427700/1‑149 (MQ=255)
                                                                      aagcaTTTATTTTGATTGTATTTTTCGTTAGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAACTAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGc                                                                                <  2:3332587/149‑1 (MQ=255)
                                                                        gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                              <  2:2664374/149‑1 (MQ=255)
                                                                           tttatttGGATTTTATTTTTCTTTCGACATTTGATTATAAATAACTGTCTTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGAAAATAGCGATTTAAACAACTATTTTTAATTTGCGCGCCAATCTGGGCGCGCACAAGCTataa                                                                           <  1:856294/149‑1 (MQ=255)
                                                                           tttatttGGATTGTATTTTTCGTTCGACCTTTTATTATACATAACTTTCTTTTTTTTAATTTATAAAATATTTCAAATAAATAAAGAAAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                           <  1:856310/149‑1 (MQ=255)
                                                                            ttatttGGAATTTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                                                          <  2:89508/149‑1 (MQ=255)
                                                                            ttatttGGAATGTTTTTTTCTTTCCACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCGATTTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataat                                                                          <  2:782506/149‑1 (MQ=255)
                                                                               tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAAATAAAAAAAATAAAAAAAAAAGAGATGTAAAGCACAAAGTTTAAATTTAGCGCCCAAAAATGGGCGCACAAGCAaaaataac                                                                       >  2:474281/1‑147 (MQ=255)
                                                                                 tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAATTAAAATAAATAAAGAAAAAAGAGATGTTAATAACTAAGTTTAAATTGAGCGCCAATCATTGCGCGCACAAGCTagaaaaaaaa                                                                     >  2:3159055/1‑143 (MQ=255)
                                                                                    aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAATAAAATATCGATTTTAAACACTAATTTTAAAATTAGCGCCCATCAATGGGCGCACAAGCAataataaaaaacg                                                                  >  2:2110713/1‑142 (MQ=255)
                                                                                    aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGAAAATATAGAATTTAATAACTAATTTTAAATTTGGCGCCCATCATTGCGCGCAAAAGATATAAAACAAAACt                                                                  >  2:2110471/1‑149 (MQ=255)
                                                                                    aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAAAAATCAAGAAAAAACCCCTGTTAAGCAAAAAGTTTTAAATTGCGGCCAAAAATGGGGCGCCCCAGCGGAAATAcccaaca                                                                  >  2:1878765/1‑144 (MQ=255)
                                                                                        tATTTTTCGTTCTACCTTTGATTATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATAGCGTTGTAAAGCAATAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                              <  1:2329141/149‑1 (MQ=255)
                                                                                        tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAAGAAAACAATAAAAAATCGATGTAAACCACAAAGTTTAAATTGCGCGCCCAACAATGCGCGCACAACCAATAAAAACAAACAAAtt                                                              >  2:3716356/1‑149 (MQ=255)
                                                                                        tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCTCACAAGCTATAATACCAACCTAAtt                                                              <  1:3461611/149‑1 (MQ=255)
                                                                                        tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTAAAAAAAAAATCAAATAAAAAAAGAAAAAATAGATTTAAAGCACAAAGTTTAAAATTAGCGGCAAACATGGCGCCCACAAGCAATAAAACAAACAAAAtt                                                              >  2:2876658/1‑149 (MQ=255)
                                                                                          tttttCGTTCGACCTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                            <  2:1462714/149‑1 (MQ=255)
                                                                                                tttctACCTTTGATTATACATAACTTTTTTTTTTTTAATTTATAAAATAATTAAAATAAATCAACAAAATAGCTATTTAAAGCACTAATTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                      <  1:2571915/148‑1 (MQ=255)
                                                                                                gttcgACCTTTGATTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATAAAAATAAATAAATAAAATCGCGTTGTAAAGCACTATTTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                      <  1:2571902/149‑1 (MQ=255)
                                                                                                   cgACGTTTGATTATACATAACTTTCTTTTTTTTAATTTATAAAATAATACAAATAAATCAAGACAATAGCGATGTAAAGCAATATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                   <  1:839274/149‑1 (MQ=255)
                                                                                                   cgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                   <  2:2809575/149‑1 (MQ=255)
                                                                                                   cgACCTTTGATTATAAATAACAGTAGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                   <  2:3890347/149‑1 (MQ=255)
                                                                                                    gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                  <  1:3133168/149‑1 (MQ=255)
                                                                                                     aCCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAATTCAAAGAAATCAAAACAAAATCGATGTAAAGCACAAAATTTAAATTTCGCGCCAAACATGGCGCGCACAAGCAATAAAACCAACCAAAATTATCCAACTCTTa                                                 >  2:3875232/1‑149 (MQ=255)
                                                                                                       cTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACATTAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  1:1012698/149‑1 (MQ=255)
                                                                                                       atttGCTTATACATAACCTTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCAATAATTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  1:3389451/148‑1 (MQ=255)
                                                                                                                taaaTAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtata                                      <  1:727452/146‑1 (MQ=255)
                                                                                                                 aaaTAACATTCTTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                     <  1:2215956/147‑1 (MQ=255)
                                                                                                                  aaTAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtc                                    <  2:718859/148‑1 (MQ=255)
                                                                                                                     aaCAGTCTTTTTTTTAATTTAAAAAATAATTCAAAGAAATCAAGACAATAGCGATGAAAACCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                 <  1:848150/149‑1 (MQ=255)
                                                                                                                     aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATTTAAAGCAATATGTTTAAAGTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                 <  1:848121/149‑1 (MQ=255)
                                                                                                                      aCAGTCGTTTTTTTAATTTATCAAATAATACAAAGAAATCAAGACAATAGCGATGAAAAGCACTAAGTTAAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                <  1:852071/149‑1 (MQ=255)
                                                                                                                        attgttttttttAATTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAACCACTATGTTTAAATTGCGCGCCAATCATGGCGCCCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                              <  1:441814/144‑1 (MQ=255)
                                                                                                                            gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                          >  2:88791/1‑149 (MQ=255)
                                                                                                                            gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                          <  1:2870319/149‑1 (MQ=255)
                                                                                                                             tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                         >  1:2591210/1‑149 (MQ=255)
                                                                                                                              ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                        <  2:2926445/149‑1 (MQ=255)
                                                                                                                                  taaaTTTATAAAATAATTAAAATAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    <  2:3252834/147‑1 (MQ=255)
                                                                                                                                   tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                   <  1:3159507/149‑1 (MQ=255)
                                                                                                                                    aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                  <  2:4047843/149‑1 (MQ=255)
                                                                                                                                    aaTTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                  <  1:1881455/149‑1 (MQ=255)
                                                                                                                                      tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                <  1:2979569/149‑1 (MQ=255)
                                                                                                                                      tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                <  2:3638317/149‑1 (MQ=255)
                                                                                                                                      tttATAAAATAATAAAAATAAATAAAGAAAATAGCGATTAAAACCAATAATTTAAAATTGCGCCCCAAACATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                <  1:2949556/149‑1 (MQ=255)
                                                                                                                                       ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt               <  2:2083107/149‑1 (MQ=255)
                                                                                                                                       ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt               <  2:1337322/149‑1 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             >  2:2488599/1‑149 (MQ=255)
                                                                                                                                           aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt           <  1:1060380/149‑1 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa          <  1:3159055/149‑1 (MQ=255)
                                                                                                                                            aaaTAATACAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa          <  1:2956258/149‑1 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt         <  1:3503211/149‑1 (MQ=255)
                                                                                                                                              ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        >  1:3979494/1‑149 (MQ=255)
                                                                                                                                              ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        >  2:2044855/1‑149 (MQ=255)
                                                                                                                                              ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        >  1:433435/1‑149 (MQ=255)
                                                                                                                                              ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTACGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc        >  2:3146601/1‑149 (MQ=255)
                                                                                                                                               taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa       >  1:1797813/1‑149 (MQ=255)
                                                                                                                                               taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa       >  2:2545982/1‑149 (MQ=255)
                                                                                                                                               taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa       >  2:2545983/1‑149 (MQ=255)
                                                                                                                                                aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa      >  2:3565709/1‑149 (MQ=255)
                                                                                                                                                 atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa     <  1:4227648/149‑1 (MQ=255)
                                                                                                                                                 atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa     >  2:2439493/1‑149 (MQ=255)
                                                                                                                                                 atTCAAAGAAATCAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCTGGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa     <  1:4227328/149‑1 (MQ=255)
                                                                                                                                                  ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt    >  2:638721/1‑149 (MQ=255)
                                                                                                                                                   tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg   <  1:3128209/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  CP009273/2463758‑2464054

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