Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,436 C→A 13.2% 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
*NC_0009132,470,4360CA13.2% 730.7 / 28.1 248Q436K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (86/122);  new base A (36/4);  total (122/126)
Fisher's exact test for biased strand distribution p-value = 5.72e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

CTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGT  >  NC_000913/2470287‑2470585
                                                                                                                                                     |                                                                                                                                                     
 tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAACa                                                                                                                                                      >  1:993836/1‑150 (MQ=255)
 tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCa                                                                                                                                                      >  1:3637053/1‑150 (MQ=255)
  ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                                     <  1:3326638/150‑1 (MQ=255)
         gtatatCTCATTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCTAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAtc                                                                                                                                              <  2:2306896/148‑2 (MQ=255)
         gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATa                                                                                                                                              >  2:2778786/1‑150 (MQ=255)
           atatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGc                                                                                                                                            >  2:830151/1‑150 (MQ=255)
                    ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                                   <  1:4553078/150‑1 (MQ=255)
                    ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                                   <  1:212286/150‑1 (MQ=255)
                     tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                                 <  1:4582055/150‑1 (MQ=255)
                      ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                                 <  1:4826634/150‑1 (MQ=255)
                      ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                                 <  1:701913/150‑1 (MQ=255)
                          gtgtCATTTTTCGGGGAATGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                             <  2:4620326/150‑1 (MQ=255)
                          gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                             <  2:5009187/150‑1 (MQ=255)
                           tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAATTAAAATAAATAAAAAAAATAAACATGTTAAGCACTAAg                                                                                                                            >  1:1838505/1‑150 (MQ=255)
                           tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAg                                                                                                                            >  2:2138896/1‑150 (MQ=255)
                            gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTGAAAGAAATCAAGTCAATAGCGATGTAGAGCACTAAGt                                                                                                                           >  1:3281252/1‑150 (MQ=255)
                            gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                                                                                                           <  2:1309912/150‑1 (MQ=255)
                            gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAATAGAGATGTAAATAACAAAGt                                                                                                                           >  1:1651472/1‑150 (MQ=255)
                             tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAATCGATGTAAAGCACTAAGtt                                                                                                                          >  2:3272369/1‑150 (MQ=255)
                             tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAATAAATAAAAACAAAAGCGATGTAAAACACAAAGtt                                                                                                                          >  2:4470251/1‑150 (MQ=255)
                             tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAATAAAAAAAATCAATAAAAAATCCAAGTAAATCAAAAAGtt                                                                                                                          >  2:4895887/1‑150 (MQ=255)
                              cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                         <  2:579543/150‑1 (MQ=255)
                               aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTTAAAGAAAACAATACAATAGCGTTGGAAATGACTAAGTTTa                                                                                                                        >  2:2661492/1‑150 (MQ=255)
                               aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAAAAATAAAATATCGATTTAAATCAATAAATTTa                                                                                                                        >  2:4362188/1‑150 (MQ=255)
                               aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTAATAAATAAATAAAAGAAATAAATAAAAAAACGATGTAAAGCACAAAGtttt                                                                                                                        >  1:4003000/1‑149 (MQ=255)
                                      ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                                 <  2:5046229/150‑1 (MQ=255)
                                      ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                                 <  2:773322/150‑1 (MQ=255)
                                         gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                              >  2:2984756/1‑150 (MQ=255)
                                            ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAAATAAAAGAAATAAATACAATATCGATGTAAATCACTAAGTTTAAATTTCGCGCCaa                                                                                                           >  1:394439/1‑150 (MQ=255)
                                            ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                           >  2:3322091/1‑150 (MQ=255)
                                               tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCatgat                                                                                                        <  2:343511/150‑5 (MQ=255)
                                                ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                                                                       <  1:2157354/150‑1 (MQ=255)
                                                  agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGGAAAGCACTAAGGTTAAATTTCGCGCCAATCATgg                                                                                                     >  2:2038489/1‑150 (MQ=255)
                                                           agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                            <  2:3616046/150‑1 (MQ=255)
                                                               agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                        <  2:4535977/150‑1 (MQ=255)
                                                                  atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataa                                                                                     <  2:3813840/150‑1 (MQ=255)
                                                                   tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAAAAAAAAAAAACAAAAAGAAGGAAAGCAAAAATTTTAAAATTAGCGCCCATCATGGCGCGCCCAAGACataaa                                                                                    >  1:4780129/1‑149 (MQ=255)
                                                                    ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                                   <  2:5038275/150‑1 (MQ=255)
                                                                    ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                                   <  2:1554039/150‑1 (MQ=255)
                                                                       tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATATCGATGTAAAGCAATAAGTTTAAAATGAGCCCCAAACAAGGCGCGCACAAACAAAAATACCa                                                                                >  2:3130989/1‑150 (MQ=255)
                                                                       tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTTTAAAAAAAATAAAAAAAAAAAAAAAAAAAGCCATTTAAAACACAAAATTTAAAATTCGCGCCAAAAAAGGCGCGCCCCACCGATAAAACCa                                                                                >  2:1345862/1‑150 (MQ=255)
                                                                        ttggtttGTATTTTTCGTTCTACCTTTGATTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATCGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCaa                                                                               <  1:4398379/144‑1 (MQ=255)
                                                                           tattttATTTTTCGTTCGACCTTTGATTATACATAACAGTCTTTTTTTTTATTTATAAAATAATTAAAATAAATAAAGAAAATCGCGATTTAAAGCACTAAGTTTAAATTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCAACCt                                                                            <  1:2661492/145‑1 (MQ=255)
                                                                             aTGAATATTTCGTACCACCTTTCCATATACATAACACACCTATTTTTAATTAAAAAAATAATAAAAACAAAACAACACAATAGCCATGAAAACCACTAATTTAAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTaa                                                                          <  1:5078873/150‑1 (MQ=255)
                                                                               ttATTTTTCTTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                                        <  1:4138081/149‑1 (MQ=255)
                                                                                  tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCAATAATACCAACCTAATTtct                                                                     >  2:3820706/1‑150 (MQ=255)
                                                                                  tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAATAAATAAAGAAAAAATCCATGTAAAGCACTAAGTTTAAATTTAGCGCCAAAAAAGGCGCGCACAAGCTATAATACCAAACAAAATtct                                                                     >  2:3473499/1‑150 (MQ=255)
                                                                                   ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                    <  2:4393527/150‑1 (MQ=255)
                                                                                   ttttCGTTAGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                    <  2:3185769/150‑1 (MQ=255)
                                                                                         ttcgACATTTTCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGAAAATCGCGATGTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcctc                                                              <  1:1345862/150‑1 (MQ=255)
                                                                                          tcgACCTTTGCTTATACATAACTGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                             <  1:1558162/150‑1 (MQ=255)
                                                                                             aCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAAACAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGAGCGCCAATCAAGGCGCGCACAAGCTATAAAAACAACCAAATTTATACACCTCTTag                                                          >  2:4457859/1‑150 (MQ=255)
                                                                                              ccTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                         <  1:310887/150‑1 (MQ=255)
                                                                                              acTTTGCTTATACATAACATTCTTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                         <  1:1906602/149‑1 (MQ=255)
                                                                                                 ttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg                                                      <  1:2817699/150‑1 (MQ=255)
                                                                                                  tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                                     <  1:2778786/150‑1 (MQ=255)
                                                                                                  tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATGGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTCCTCCTCCTCTTAGAGTGa                                                     <  1:2475175/150‑1 (MQ=255)
                                                                                                   gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                    <  2:4707144/150‑1 (MQ=255)
                                                                                                   gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTTAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAAACCAACCAAATTTCTCCTGCTCTGAGAGTGAc                                                    >  2:2273885/1‑150 (MQ=255)
                                                                                                    cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                   <  1:1758527/150‑1 (MQ=255)
                                                                                                      tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGATTGACtat                                                 >  2:3392767/1‑150 (MQ=255)
                                                                                                      tataCATAACAGTCGTTTTTTTAATTTATAAAATAAGTAAAACAAATAAATCCACTCTCGGTGTGAAGCCCTAAGTTTAAATTGCGCGCCCAACGTGGCGCGCACCACCGATAGAAACCACCGAATATCTACACCGCGGAGAGTGAcgaa                                                 >  2:2267603/1‑147 (MQ=255)
                                                                                                           aTAACCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtctc                                            <  1:2063656/150‑1 (MQ=255)
                                                                                                               cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                        <  2:2206531/150‑1 (MQ=255)
                                                                                                                 gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                      <  2:1363270/150‑1 (MQ=255)
                                                                                                                 gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                      <  2:3942409/150‑1 (MQ=255)
                                                                                                                  tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                                     >  2:1762720/1‑150 (MQ=255)
                                                                                                                    gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                                   >  2:1340074/1‑150 (MQ=255)
                                                                                                                     tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                  >  2:804138/1‑150 (MQ=255)
                                                                                                                        ttttAATTTATAAAATAATTCAAATAAATCAAGACAATTGCGATTTAATGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                               <  1:4638055/150‑1 (MQ=255)
                                                                                                                              tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                         <  2:1431760/150‑1 (MQ=255)
                                                                                                                              tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                         >  1:296379/1‑150 (MQ=255)
                                                                                                                                tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                       <  2:2497846/150‑1 (MQ=255)
                                                                                                                                 ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                      >  2:1716734/1‑150 (MQ=255)
                                                                                                                                 ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                      <  2:3925180/150‑1 (MQ=255)
                                                                                                                                 ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                      >  2:4277652/1‑150 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                     >  1:274939/1‑150 (MQ=255)
                                                                                                                                  taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                     >  2:1060518/1‑150 (MQ=255)
                                                                                                                                   aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa                    >  2:1222464/1‑150 (MQ=255)
                                                                                                                                   aaaaTAATTCAAAGAAATCAAGACAATAGCGAGGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCGGATAGAAATGCGGTATTGACTa                    >  2:3004893/1‑150 (MQ=255)
                                                                                                                                    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                   >  1:4344603/1‑150 (MQ=255)
                                                                                                                                    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                   >  2:3066297/1‑150 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                  <  1:3804434/150‑1 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTATTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                  <  2:3605408/150‑1 (MQ=255)
                                                                                                                                        aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                      >  2:1108027/1‑143 (MQ=255)
                                                                                                                                          ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg             >  2:2407135/1‑150 (MQ=255)
                                                                                                                                            cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc           >  1:1026064/1‑150 (MQ=255)
                                                                                                                                             aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc          <  2:2123597/150‑1 (MQ=255)
                                                                                                                                             aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc          >  1:1310113/1‑150 (MQ=255)
                                                                                                                                              aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt         >  2:1247641/1‑150 (MQ=255)
                                                                                                                                              aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt         >  1:444640/1‑150 (MQ=255)
                                                                                                                                              aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt         <  2:3733501/150‑1 (MQ=255)
                                                                                                                                              aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt         >  1:3710233/1‑150 (MQ=255)
                                                                                                                                              aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt         >  2:524096/1‑150 (MQ=255)
                                                                                                                                               aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTGGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg        >  1:3468145/1‑150 (MQ=255)
                                                                                                                                               aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg        >  1:7890/1‑150 (MQ=255)
                                                                                                                                               aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg        >  1:5075400/1‑150 (MQ=255)
                                                                                                                                               aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg        >  1:1096180/1‑150 (MQ=255)
                                                                                                                                                    tCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCg   >  2:2359732/1‑150 (MQ=255)
                                                                                                                                                     cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt  <  1:4930207/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGT  >  NC_000913/2470287‑2470585

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: