Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,463,906 C→A 23.8% 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,9060CA23.8% 152.1 / 16.9 75Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (16/39);  new base A (14/6);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 2.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  CP009273/2463758‑2464055
                                                                                                                                                    |                                                                                                                                                      
tGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAAATCa                                                                                                                                                       >  2:4241659/1‑150 (MQ=255)
    gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTATTTTTTAAAAAAAATaaaata                                                                                                                                                   >  1:2441104/1‑148 (MQ=255)
        ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGATTATACATAACAGTCGTTTTTTTTATTTAtaaaaaaaataaaaaaaaaa                                                                                                                                               >  1:4246751/1‑148 (MQ=255)
             tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAataaaataaaaaaaaaa                                                                                                                                          >  2:1696467/1‑143 (MQ=255)
               tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAAAAAAAAAGaaaa                                                                                                                                        >  2:956003/1‑150 (MQ=255)
               tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAATAAATAAAAAAAAAAAAAACaa                                                                                                                                        >  2:2916345/1‑150 (MQ=255)
                ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                                                                                                       <  1:2593804/150‑1 (MQ=255)
                 gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                                                                                                      <  2:4229282/150‑1 (MQ=255)
                 gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAAAAATTAAAAAAAAAAAAAACaaaa                                                                                                                                      >  1:454671/1‑148 (MQ=255)
                     atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                  <  1:2868048/150‑1 (MQ=255)
                     atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                  <  1:2869296/150‑1 (MQ=255)
                                aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                                                                                                       <  1:1967409/150‑1 (MQ=255)
                                  gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                     <  1:676391/150‑1 (MQ=255)
                                   tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                                                                                                    <  2:1433458/150‑1 (MQ=255)
                                         ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                              <  1:812258/150‑1 (MQ=255)
                                          tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                                                                                                             <  1:3686873/150‑1 (MQ=255)
                                            tCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAAATAAAATAAAAAAAAAAAAAAAAGCTGTAAAAAACtaaatttaaattt                                                                                                           >  2:4176047/1‑149 (MQ=255)
                                                ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                                                                                                       <  2:998707/150‑1 (MQ=255)
                                                  aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAAAAAAAAAAAAAGCCAATTAAATCACAAAGTTTAAATTGAGCGcc                                                                                                    >  2:2766077/1‑150 (MQ=255)
                                                   aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                   <  1:496737/150‑1 (MQ=255)
                                                     gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTAAAAATAAATAAAATAAAAAAAAAAAAAATCGAAGTAAAAAACTAAAGTTAAATTTGGCGCCaca                                                                                                 >  2:3197466/1‑148 (MQ=255)
                                                     gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                                                 <  2:1469408/150‑1 (MQ=255)
                                                     gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAAACAAAAGCGATGTAAATAAATAAGTTTAAATTGAGCGCCAAt                                                                                                 >  2:764046/1‑150 (MQ=255)
                                                            agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATAAAAAAAAAATAGAATTAAATCAAAAAGTTTTAATTTAGCGGCAAACAAGgcg                                                                                          >  2:4645283/1‑150 (MQ=255)
                                                                   agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTAATTTTAAAAAAATTAAAATAAATAAATACAAAATCGAGGTAAAACACTAAGTTTAAATTTCGCGCCCAAACAGGCGCGCACaa                                                                                   >  2:2641296/1‑150 (MQ=255)
                                                                   agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAAAAAATAAAAAAAAAAAAAAAAAAAACCATGTAAAAAAAAAAAGTTAAATTTAGCGCCCAACAAAGCGCCCACaa                                                                                   >  1:4782734/1‑150 (MQ=255)
                                                                    gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATCAAATAAATCAAGACAAAAGAGATGTAAAGCACAAAGTTTAAATTGCGCGCCAATCAATGCGCGCACaac                                                                                  >  2:2386608/1‑149 (MQ=255)
                                                                                ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAATAAAGACAATAACGAATTAAAGAACAAAGTTTTAAATTAGCGCCAATAAATGCGCCCACAAACAAAAAAACCaa                                                                      >  2:3485516/1‑150 (MQ=255)
                                                                                ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAAAAATAAATAAAAAAGCGAAGTAAAAAAAGAAGTTTAAATTTCTCGCCCAACAAGTCGCGCCCAACCAagaataaaaa                                                                      >  2:3484289/1‑146 (MQ=255)
                                                                                 tgtattttATTTTTCGTTCGACCTTTGATTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                     <  1:1932529/143‑1 (MQ=255)
                                                                                 tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                     <  1:1578662/150‑1 (MQ=255)
                                                                                  ggATTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTATTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGTTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                    <  2:2582632/150‑1 (MQ=255)
                                                                                  ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTGGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                    <  2:152928/150‑1 (MQ=255)
                                                                                     aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTTTAAAAAAAATAAAAAAAAAAAAAAAAAAAACGAAGTAAATGACCAAGTTTAAAATTTGCGCCAAGAACGGGGCGCCCAACCAAAAAAAACAACCaaa                                                                 >  1:4356341/1‑150 (MQ=255)
                                                                                       ttATTTTTCGTTCGACCTTTGATTATAAATAACTGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCGTTGTAAAGCAATAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                               <  1:2236936/149‑1 (MQ=255)
                                                                                       gTATTTTTCGTTCGACATTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                               <  1:3541199/150‑1 (MQ=255)
                                                                                        tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                              <  2:550221/150‑1 (MQ=255)
                                                                                         aTTTTTCGTTCTACATTTGATTATACATAACATTCTTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                             <  1:1066924/150‑1 (MQ=255)
                                                                                            tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                          <  2:3198516/150‑1 (MQ=255)
                                                                                            tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                          <  2:3197930/150‑1 (MQ=255)
                                                                                             ttcgttctACATTTGATTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTAAAATAAATCAATACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcct                                                         <  1:4349162/150‑1 (MQ=255)
                                                                                             ttcgttcgACCTTTGATTATACATAACTTTCTTTTTTTTAATTTATAAAATATTAAAAATAAATAAAGACAATAGCGTTTTAAAGCAATAATTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcct                                                         <  1:1289896/150‑1 (MQ=255)
                                                                                               cgttctACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  1:4698177/150‑1 (MQ=255)
                                                                                               cgttcgACCTTTTCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  1:4698169/150‑1 (MQ=255)
                                                                                                gttcgACTTTTTATTATACATAACTGTATTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAATAGCTTTTTAAAGAACTATTTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                      <  1:4508345/150‑1 (MQ=255)
                                                                                                gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATGGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                      <  1:4319994/150‑1 (MQ=255)
                                                                                                  tctacttttGATTATACATAAAATTATTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAATAGCGATTTAAAGCAATAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                    <  1:2536181/144‑1 (MQ=255)
                                                                                                   cgACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                   <  2:454671/150‑1 (MQ=255)
                                                                                                   cgACCTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCtt                                                   <  2:2888215/150‑1 (MQ=255)
                                                                                                     acatttGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTag                                                 <  1:3073903/147‑1 (MQ=255)
                                                                                                          tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  2:4319970/150‑1 (MQ=255)
                                                                                                          tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  2:4318847/150‑1 (MQ=255)
                                                                                                             ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                         <  2:2728573/150‑1 (MQ=255)
                                                                                                             ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                         <  2:297751/150‑1 (MQ=255)
                                                                                                                taCATAACAGTCGTTTTTTTAATTTAAAAAATAATTCAAATAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                      <  1:3038940/150‑1 (MQ=255)
                                                                                                                        aGTCGTTTTTTTAATTTATAAAATAATTCAAATAAAACAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCAATAATACCAAACTAATTTCTCCTCCTCTTAGATTTACTATATATCCAGAt                              >  2:2418306/1‑150 (MQ=255)
                                                                                                                        aGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAACAAAAAAATATCGATGTAAAGCACAAAGTTTAAATTTCGCGCCAAAAATTGCGCGCAAAACATATAAAAACAAACAAATTTCTACAACAATTAGAGTTAAAATATAACCTGAt                              >  2:4293806/1‑150 (MQ=255)
                                                                                                                          tcttttttttAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCTATTTAAAGCAATAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            <  1:2132144/147‑1 (MQ=255)
                                                                                                                          acttttttttAATTTATAAAATAATTCAAATAAATAAATAAAATAGCTATGAAAAGCAATAAGTTTAAATTGCGCGCCAATCAGGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            <  1:9857/147‑1 (MQ=255)
                                                                                                                           cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                           >  2:4650759/1‑150 (MQ=255)
                                                                                                                           cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                           >  2:4650725/1‑150 (MQ=255)
                                                                                                                              ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                        <  1:4241659/150‑1 (MQ=255)
                                                                                                                              ttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAtt                        >  2:4633646/1‑150 (MQ=255)
                                                                                                                               tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       >  1:1629364/1‑150 (MQ=255)
                                                                                                                               tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:3475749/150‑1 (MQ=255)
                                                                                                                                 tttAATTTATAAAATAATTAAAATAAATAAATACAATACCGATGTAAATCACTAATTTTAATTCCCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:614506/150‑1 (MQ=255)
                                                                                                                                   tAATTTAAAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                   <  1:3372524/150‑1 (MQ=255)
                                                                                                                                      tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                >  1:963784/1‑150 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  2:2383416/1‑150 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:2017082/1‑150 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:2037725/1‑150 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:2037727/1‑150 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:2037761/1‑150 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:54635/1‑150 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         >  1:885033/1‑150 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         >  1:884895/1‑150 (MQ=255)
                                                                                                                                                 atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt     <  2:3193930/150‑1 (MQ=255)
                                                                                                                                                  ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  1:3099535/150‑1 (MQ=255)
                                                                                                                                                    cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  >  1:534323/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  CP009273/2463758‑2464055

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: