Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,463,906 C→A 16.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,9060CA16.9% 221.1 / 14.7 87Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (23/47);  new base A (14/3);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 2.86e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  CP009273/2463766‑2464055
                                                                                                                                            |                                                                                                                                                     
ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAATaaaaa                                                                                                                                              >  2:4284113/1‑148 (MQ=255)
  ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTTTTTATAAAaaaaataaaaaaaataaa                                                                                                                                            >  1:1989432/1‑150 (MQ=255)
     tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAAtaaataaaaaa                                                                                                                                         >  2:3343885/1‑147 (MQ=255)
              tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  1:4284113/150‑1 (MQ=255)
              tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAATTAAAATAAAACAACACaaaaaccat                                                                                                                                >  2:1047721/1‑143 (MQ=255)
                   aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           <  1:3123592/150‑1 (MQ=255)
                    ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAAAAAATAAAAGAAATAAAGAAAAAATCGATTTaaaa                                                                                                                          >  1:3949045/1‑149 (MQ=255)
                      ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  1:4488255/150‑1 (MQ=255)
                       tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAATAAAAAAATAAAAAAGCGATTTTAAGCAc                                                                                                                       >  2:391800/1‑150 (MQ=255)
                        aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTTAAATAATTAAAAAAAAAAAATACAAAATCCCTGTAAATCaat                                                                                                                      >  2:3371300/1‑148 (MQ=255)
                             tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                 <  1:4538180/150‑1 (MQ=255)
                             tCATTTTTCGGGGAAGGATATAGAGATGTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                 <  1:2155256/150‑1 (MQ=255)
                                      ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAATCAAGACAAAAGAGATTTAAATCACAAATTTTAAATTtcg                                                                                                        >  2:657730/1‑150 (MQ=255)
                                             gATATAGAGATTTAAGCAAGCATTTATTTGGAATTTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                                                 <  1:2421674/150‑1 (MQ=255)
                                             gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                                                 <  2:265552/150‑1 (MQ=255)
                                                ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAATAAAAAATAGATTTTAATAAAAAAGTTTTAATTTCGCGCCaataat                                                                                              >  2:1137868/1‑150 (MQ=255)
                                                ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAAACAATAACGATTTTAATAACTAAGTTTTAATTTGGCGCCAAACAt                                                                                              >  2:1904075/1‑150 (MQ=255)
                                                      aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg                                                                                        <  1:1790755/150‑1 (MQ=255)
                                                         tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAAACAAGAACAAAGCGATGTTAAGCACGAAGTTTAAATTTAGCGGCAAAAAGGgcgccccc                                                                                     >  2:4246435/1‑146 (MQ=255)
                                                         tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAAGAAATAAACACAAAATCGATTTAAAGCACGAATTTTTAATTGCTCGCCAAAAAAGGCGCGcac                                                                                     >  2:4246892/1‑150 (MQ=255)
                                                          aagcaagcATTTATTTGGAATGTATTTTACGTACGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca                                                                                    <  2:146738/150‑1 (MQ=255)
                                                                 caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAAAAAATAAAAAAAAAAAAAACAAAATAGATGTAAAAAACAAAGTTTAAATTTAGAGCCAAAAATGGCGCGCACaaccaaaa                                                                             >  1:1851137/1‑144 (MQ=255)
                                                                     tatttGTATTGTATTTTTCGTTCGACCTTTTCTTATAAATAAATGTCTTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGAAAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATAc                                                                         <  1:588342/150‑1 (MQ=255)
                                                                         tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATAAAAGAAATAAAAACAAAAAAGATTTAAATAAATAAATTTAAATTTAGCGCCAATAAATGCGCGCCCAAGCTATAAAAACaaa                                                                     >  2:4175070/1‑149 (MQ=255)
                                                                             aTGTATTTTTCTTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGGGCGCACAAGCTATAATACCAACCTaa                                                                 <  2:4227217/150‑1 (MQ=255)
                                                                             aTGTATTTTTCGTTCTACCTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTaa                                                                 <  2:4425735/150‑1 (MQ=255)
                                                                                    tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAATAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAAATTCGCGCCCAAAATGGCGCGCACACACAATAAAACCAACCAAAattaaca                                                          >  2:1947430/1‑146 (MQ=255)
                                                                                      tcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                        <  1:2283987/150‑1 (MQ=255)
                                                                                       cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  2:791465/150‑1 (MQ=255)
                                                                                        gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                      <  2:2657699/150‑1 (MQ=255)
                                                                                         ttcgACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcctc                                                     <  1:2281897/150‑1 (MQ=255)
                                                                                         ttcgACCTTTGATTATAAATAACATTCTTTTTTTTAATTTATAAAATAATACAAATAAATCAAGAAAATAGCTATGTAAAGCACTATTTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcctc                                                     <  1:2930411/150‑1 (MQ=255)
                                                                                          tcgACATTTGATTATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCAATAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                    <  1:1137868/150‑1 (MQ=255)
                                                                                           cgACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                   <  2:2013223/150‑1 (MQ=255)
                                                                                           cgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                   <  2:3382639/150‑1 (MQ=255)
                                                                                           cgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                   <  2:3382332/150‑1 (MQ=255)
                                                                                           cgACATTTGCTTATACATAACATTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                   <  1:3046338/150‑1 (MQ=255)
                                                                                                  tGCTTATACATAACAGTCGTGTTTTTAATTTTTAAAAAAAATAAAAAAAAACAAGAAAAAAGCGAAGTAAAGAACAAAGATTAAATTTCGCGCCAATCATGGCGCGCACAACCAAAAAAAACAAACAAAATTATCCTACACTAAGAGTgg                                            >  2:2041264/1‑149 (MQ=255)
                                                                                                             acctgTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAACCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                 <  1:4050264/146‑1 (MQ=255)
                                                                                                             aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATCGCTATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                 <  1:4049518/150‑1 (MQ=255)
                                                                                                             aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                 <  2:4549346/150‑1 (MQ=255)
                                                                                                             aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                 <  2:4549335/150‑1 (MQ=255)
                                                                                                               ccgTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATGGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                               <  1:3224979/148‑1 (MQ=255)
                                                                                                               cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                               <  2:2927604/150‑1 (MQ=255)
                                                                                                               cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                               <  2:2924777/150‑1 (MQ=255)
                                                                                                                  tcttttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            <  1:1808808/147‑1 (MQ=255)
                                                                                                                  tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATCGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            <  1:328278/150‑1 (MQ=255)
                                                                                                                  tCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            <  2:2364110/150‑1 (MQ=255)
                                                                                                                    ttttttttATTTTAAAAAATAATACAAATAAATCAAGACAATCGCTATTTAAAGAAATAATTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          <  1:3979140/149‑1 (MQ=255)
                                                                                                                    gTTTTTTTGATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTAAATCTCCTGATAGaa                          >  2:1975877/1‑150 (MQ=255)
                                                                                                                      ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                        <  2:1851137/150‑1 (MQ=255)
                                                                                                                      ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                        <  2:1530171/150‑1 (MQ=255)
                                                                                                                      ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                        <  2:1852074/150‑1 (MQ=255)
                                                                                                                        ttttattttaTAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  1:4439540/150‑1 (MQ=255)
                                                                                                                        ttttAATTTATAAAATAATTCAAATAAATCAAGAAAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  1:4439420/150‑1 (MQ=255)
                                                                                                                        ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAACGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  2:2538002/150‑1 (MQ=255)
                                                                                                                        ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  2:2538145/150‑1 (MQ=255)
                                                                                                                            aaTTTATAAAATAATTCAAATAAAACAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                  <  1:1891640/150‑1 (MQ=255)
                                                                                                                              tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  2:2989080/150‑1 (MQ=255)
                                                                                                                              tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  2:2988954/150‑1 (MQ=255)
                                                                                                                                tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              >  1:4293313/1‑150 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            >  2:403474/1‑150 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            >  2:403470/1‑150 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            <  2:601186/150‑1 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            <  2:601263/150‑1 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            >  2:2355800/1‑150 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            <  2:601274/150‑1 (MQ=255)
                                                                                                                                  aaaaaTAATTCAAATAAATCAAGAAAATAGCTATGAAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            <  1:1184634/149‑1 (MQ=255)
                                                                                                                                  aaaaaTAATAAAAATAAATCAATAAAATCGCGATTTAAACAAATAAGTTAAATTTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            <  1:1184773/149‑1 (MQ=255)
                                                                                                                                   aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTCTTGACTa           >  1:4192621/1‑150 (MQ=255)
                                                                                                                                   aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa           >  1:4636387/1‑150 (MQ=255)
                                                                                                                                   aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa           >  1:3395547/1‑150 (MQ=255)
                                                                                                                                   aaaaTAAGTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa           >  2:2549659/1‑150 (MQ=255)
                                                                                                                                    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:2368458/1‑150 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATGCCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         <  2:2653148/150‑1 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         >  1:2393458/1‑150 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         >  1:2334229/1‑150 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         >  2:4354766/1‑150 (MQ=255)
                                                                                                                                      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGTTAGAATTGCGGTATTGACTATCa        >  1:404365/1‑150 (MQ=255)
                                                                                                                                      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  2:1957748/1‑150 (MQ=255)
                                                                                                                                      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  2:1957873/1‑150 (MQ=255)
                                                                                                                                      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  1:1417837/1‑150 (MQ=255)
                                                                                                                                      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  1:816996/1‑150 (MQ=255)
                                                                                                                                      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  1:816949/1‑150 (MQ=255)
                                                                                                                                       taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa       >  2:2134241/1‑150 (MQ=255)
                                                                                                                                        aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      >  2:4102672/1‑150 (MQ=255)
                                                                                                                                        aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  2:3767705/150‑1 (MQ=255)
                                                                                                                                           tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc   >  1:862901/1‑150 (MQ=255)
                                                                                                                                            cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  >  2:3735362/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  CP009273/2463766‑2464055

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: