Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,463,915 C→A 28.4% Q436K (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,9150CA28.4% 145.7 / 30.1 70Q436K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base A (18/3);  total (44/26)
Fisher's exact test for biased strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATT  >  CP009273/2463775‑2464061
                                                                                                                                            |                                                                                                                                                   
gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                                                                                                            <  1:2129993/150‑1 (MQ=255)
 aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAAAAAGaaaaaaa                                                                                                                                           >  2:5039475/1‑147 (MQ=255)
   tatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTATTTTATTAAATAAATAAAATAAAAGAAAAAAATAgag                                                                                                                                         >  2:4289059/1‑148 (MQ=255)
    atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                        <  2:3662786/150‑1 (MQ=255)
    atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                        <  2:1038467/150‑1 (MQ=255)
     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                       <  1:1709138/150‑1 (MQ=255)
     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTTTAAAAAAAATCAAATAAATAAAAACaaaaaccaa                                                                                                                                       >  2:4596379/1‑143 (MQ=255)
     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTCTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                       <  1:1709184/150‑1 (MQ=255)
      ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTTTTTTTAAAAAAATTAAAATAAAAAAAAACAATatagatt                                                                                                                                      >  2:2720201/1‑144 (MQ=255)
        cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                                    <  1:521191/150‑1 (MQ=255)
         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAAAAAAGACAATAAAGATGTaa                                                                                                                                   >  2:3031913/1‑150 (MQ=255)
         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAATAAAAAAAGACAAAAACGCTGtta                                                                                                                                   >  2:3031929/1‑148 (MQ=255)
             ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTTAAATTATTAAAATAAAAAAATACAAAAAAGATTTAAATCa                                                                                                                               >  1:1605641/1‑150 (MQ=255)
                tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAACAAAAAAAAGATGTaaaaaacaa                                                                                                                            >  2:3305656/1‑144 (MQ=255)
                     cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                       <  2:1490661/150‑1 (MQ=255)
                     cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                       <  2:1490373/150‑1 (MQ=255)
                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATTAAAAAATTAAAATAAAAAAATAAAAAATAGATTTAAATAACAAAGTTTTAAtt                                                                                                                 >  1:861825/1‑150 (MQ=255)
                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTTTTTATAAAATAAATAAAATAAAAAAAAACAAAAACCATTTAAATAAAAAAATTTAAAta                                                                                                                 >  1:5007927/1‑149 (MQ=255)
                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAATAAATAAAAAAAAAAAAAAAAAAAAAGATTTAAATCAAAAATTTTAAAtt                                                                                                                 >  1:3904546/1‑150 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTAAAACAAATAAACACAAAAACCAGTTAAAGCACTAAGATTAAAtatagcgc                                                                                                           >  2:5171261/1‑143 (MQ=255)
                                gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAATAAATCAATAAAAAAACAATTTAAAAAAAAAAGTTTAAATTtcgcgc                                                                                                           >  2:5171328/1‑150 (MQ=255)
                                        taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAATACAAAACAGATGTAAAGCACAAAGTTTTAAATTCGCGCCaaaaaat                                                                                                   >  2:2505726/1‑145 (MQ=255)
                                                    aaagaaTTTTTTTGGATTGTATTTTTCGTTCGACATTTGCTTATAAATAACATTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATAAATAAAATAGCTATTTAAAGCACTAATTTTAATTTGCGCGCCAATCATGGCGCGCACAAGc                                                                                       <  1:5139681/149‑1 (MQ=255)
                                                     aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTTTTTTTAAAAAAAATAAAAAAAAAAAAAAAAAAAAAGAAGTAAAGAACAAAATTTTAATATTGCGCCACAAATGGCGCGcccaccca                                                                                      >  1:3307868/1‑146 (MQ=255)
                                                           ttatttGTAATTTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                                <  2:1501495/150‑1 (MQ=255)
                                                           ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                                <  2:1501489/150‑1 (MQ=255)
                                                               ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTAAAAAAAAAATAAAACAAAAAAATAAAAAATCCAAGTAAACACCAAAGTTTAAAATTGCCGCCAATAAAGGCGCGCCCAAACAATAAAAACaa                                                                            >  2:1527493/1‑150 (MQ=255)
                                                               ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAAAAATAAAGAAAAAATAGCTGTAAAGCACTAAGTTTAAAATTAGCGCCAAACATGGCGCGCCCAAACAATAAAACaaa                                                                            >  2:4164440/1‑150 (MQ=255)
                                                                      ttATTTTTCGTTCGACTTTTGATTATAAATAAAAGTCTTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAATAGCTATGTAAAGCAATATGTTTAATTTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                                     <  1:240409/149‑1 (MQ=255)
                                                                      gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                                     <  2:4457390/150‑1 (MQ=255)
                                                                      gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAAAAATACAAAATAGATGTAAATAACAAATTTTAAATTTAGCGCCAAAAAATGCGCGCACAAACAATAAAAACaaacaaaaa                                                                     >  2:2002934/1‑143 (MQ=255)
                                                                      gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAATTAAAATAAAAAAAAAAAAAATAGATTTAAAACAAAAAGTTTAAATATAGCGCCAAAAAATGCGCGCACAAGATATAAAAACaaacaaaat                                                                     >  2:2443304/1‑143 (MQ=255)
                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAAAAATAAATAAAATATCGATTTTTATAAATAAGTTTAAATTTGGGGCCAAAAATGTGGGGCGCAAGCTATAATAAAAAACAAAtttata                                                                 >  2:5021474/1‑147 (MQ=255)
                                                                           tttGGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                <  2:4920102/150‑1 (MQ=255)
                                                                              cgttcgACCTTTGCTTATACATAACATTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                             <  2:1416004/150‑1 (MQ=255)
                                                                              cgttcgACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                             <  2:1416003/150‑1 (MQ=255)
                                                                                             tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                              <  1:1596115/150‑1 (MQ=255)
                                                                                                   taaccttCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                        <  1:2841671/144‑1 (MQ=255)
                                                                                                          cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                                 >  1:5149386/1‑150 (MQ=255)
                                                                                                             ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                              >  1:2013191/1‑150 (MQ=255)
                                                                                                             ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                              <  1:435445/150‑1 (MQ=255)
                                                                                                             ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                              <  1:3428788/150‑1 (MQ=255)
                                                                                                             ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                              <  1:433766/150‑1 (MQ=255)
                                                                                                                      ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                     >  1:745911/1‑150 (MQ=255)
                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                   >  1:4220050/1‑150 (MQ=255)
                                                                                                                         taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                  >  2:410333/1‑150 (MQ=255)
                                                                                                                           aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATCATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                >  1:1840195/1‑150 (MQ=255)
                                                                                                                             ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa              <  2:861825/150‑1 (MQ=255)
                                                                                                                             ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCCCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa              >  2:2906453/1‑150 (MQ=255)
                                                                                                                              taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGTCTATATCTCCTGATAGAATTGCGGTATTGACTATCaa             >  1:4179155/1‑150 (MQ=255)
                                                                                                                              taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa             >  1:1351381/1‑150 (MQ=255)
                                                                                                                              taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa             >  1:1351432/1‑150 (MQ=255)
                                                                                                                              taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa             >  1:1351442/1‑150 (MQ=255)
                                                                                                                              taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGGATTGACTATcca             >  2:157558/1‑148 (MQ=255)
                                                                                                                                atacAAAGAAATAAACACAATAGCGATGAAAAGCACTAAGTTTAAATTGCGCGCCAATCAGGGCCCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt           <  1:4816435/147‑1 (MQ=255)
                                                                                                                                  tCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc         <  1:5171328/150‑1 (MQ=255)
                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAATTGc         >  2:3348936/1‑150 (MQ=255)
                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc         <  1:5171261/150‑1 (MQ=255)
                                                                                                                                  tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc         >  2:3347235/1‑150 (MQ=255)
                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc        >  1:2360523/1‑150 (MQ=255)
                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc        >  1:2228338/1‑150 (MQ=255)
                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc        >  2:2116494/1‑150 (MQ=255)
                                                                                                                                    aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc       >  1:3905733/1‑150 (MQ=255)
                                                                                                                                     aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGACTTGCGGTATTGACTATCAAATGCCCt      >  2:2082185/1‑150 (MQ=255)
                                                                                                                                     aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt      >  1:1501236/1‑150 (MQ=255)
                                                                                                                                     aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt      >  2:2081088/1‑150 (MQ=255)
                                                                                                                                       gAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa    <  2:4330960/150‑1 (MQ=255)
                                                                                                                                        aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTCTCAAATGCCCTGAt   >  2:2644821/1‑150 (MQ=255)
                                                                                                                                        aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt   >  1:2507893/1‑150 (MQ=255)
                                                                                                                                         aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt  >  2:3612670/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATT  >  CP009273/2463775‑2464061

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