Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 13.7% 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
*NC_0009132,470,4270CA13.7% 675.6 / 23.4 241Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (67/131);  new base A (36/4);  total (106/135)
Fisher's exact test for biased strand distribution p-value = 3.06e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470280‑2470576
                                                                                                                                                    |                                                                                                                                                     
gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaa                                                                                                                                                     <  1:512741/150‑1 (MQ=255)
gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaa                                                                                                                                                     <  1:2102780/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaaa                                                                                                                                                    >  1:2492284/1‑150 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:3896897/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  2:2541886/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  2:2541733/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  2:1189656/150‑1 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  1:4641346/150‑1 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaataaat                                                                                                                                               >  1:2960671/1‑150 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAt                                                                                                                                               >  1:2702707/1‑150 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTTTAAAATAATTAAAAGAAAtaa                                                                                                                                             >  1:3637067/1‑148 (MQ=255)
         ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                            <  1:3326717/150‑1 (MQ=255)
           ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                          <  1:1880425/150‑1 (MQ=255)
                  atatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGc                                                                                                                                   >  2:830151/1‑150 (MQ=255)
                   tatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                                                                                                  <  2:2338080/150‑1 (MQ=255)
                    ataacAATTTTTTTTTTCATTTTTTGGGGTATGCTTTTTACTTTTAATCAACCAATTATTTTGTAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                                                                                                  <  2:1248690/146‑1 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                              <  2:3806107/150‑1 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAAAAAAGACAAAAGAGAttt                                                                                                                              >  1:2743866/1‑148 (MQ=255)
                          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           <  2:2920501/150‑1 (MQ=255)
                            tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                                                                                                         <  1:2416677/150‑1 (MQ=255)
                            tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTA‑TTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  1:4582055/150‑1 (MQ=255)
                             ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAATTAAAAGAAATAAAGACAATAGCGATGTAAATCa                                                                                                                        >  1:4062793/1‑150 (MQ=255)
                             ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  1:4826634/150‑1 (MQ=255)
                                 gtgtCATTTTTCGGGGAATGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                    <  2:4620326/150‑1 (MQ=255)
                                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                    <  2:5009187/150‑1 (MQ=255)
                                  tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAg                                                                                                                   >  2:2138896/1‑150 (MQ=255)
                                   gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTGAAAGAAATCAAGTCAATAGCGATGTAGAGCACTAAGt                                                                                                                  >  1:3281252/1‑150 (MQ=255)
                                   gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAATAAAGAAAATAGCGATGTAAATAACAAAGt                                                                                                                  >  1:4569152/1‑150 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATATCGATGTAAAGCACTAAGtt                                                                                                                 >  2:3420350/1‑150 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAATAAATAAAAACAAAAGCGATGTAAAACACAAAGtt                                                                                                                 >  2:4470251/1‑150 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAATAAAAAAAATCAATAAAAAATCCAAGTAAATCAAAAAGtt                                                                                                                 >  2:4895887/1‑150 (MQ=255)
                                     cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                <  2:579543/150‑1 (MQ=255)
                                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTAATAAATAAATAAAAGAAATAAATAAAAAAACGATGTAAAGCACAAAGtttt                                                                                                               >  1:4003000/1‑149 (MQ=255)
                                       tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                                                                                              <  1:929368/150‑1 (MQ=255)
                                        ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  2:82748/150‑1 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAATACAATATCGATGTAAATCACGAAGTTTAAATTtcg                                                                                                        >  1:2976148/1‑150 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                        <  2:773322/150‑1 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                        <  2:5046229/150‑1 (MQ=255)
                                               ggAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                                                                                                      <  1:4693157/150‑1 (MQ=255)
                                               ggAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                                                                                                      <  1:1743965/150‑1 (MQ=255)
                                                gAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                                                                                                     >  2:2984661/1‑150 (MQ=255)
                                                 aaGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                                                                                                    <  1:2816966/150‑1 (MQ=255)
                                                   ggATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAAATAAAAGAAATAAATACAATATCGATGTAAATCACTAAGTTTAAATTTCGCGCCaa                                                                                                  >  1:394439/1‑150 (MQ=255)
                                                   ggATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGAAAAAAGCCATGTAAAGCACTAAGTTTAAATTTCGCGACaa                                                                                                  >  2:1101326/1‑150 (MQ=255)
                                                   ggATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                  >  2:3322091/1‑150 (MQ=255)
                                                      tataGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCatgat                                                                                               <  1:3321442/150‑5 (MQ=255)
                                                      tataGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCatgat                                                                                               <  2:343511/150‑5 (MQ=255)
                                                         agagATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                            <  1:2134299/150‑1 (MQ=255)
                                                         agagATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                            >  2:1224509/1‑150 (MQ=255)
                                                         agagATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGGAAAGCACTAAGGTTAAATTTCGCGCCAATCATgg                                                                                            >  2:2038489/1‑150 (MQ=255)
                                                           agaTTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAAGAAAATATCGATGTAAATCACTAAGTTTAAATTGCGCGCCAATAATGgcg                                                                                          >  2:30369/1‑150 (MQ=255)
                                                                 aagcaagcATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAAAAAATAAAAAAAAAAACGATGTAAAGCACAAAATTTAAATTTCGCGCCAATCATGTCGCGccca                                                                                    >  2:1058803/1‑150 (MQ=255)
                                                                  agcaagcaTTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                   <  2:3616046/150‑1 (MQ=255)
                                                                  agcaagcaTTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATAAAGACAATATAGATGTAAAGCAAAAAGTTTAAATTTAGCGCCCATAATGGCGCGCACaa                                                                                   >  2:1333201/1‑150 (MQ=255)
                                                                   gcaagcaTTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaattt                                                                                  <  2:1356950/150‑5 (MQ=255)
                                                                      agcaTTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                               <  2:4536008/150‑1 (MQ=255)
                                                                           tta‑tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                          <  2:5038275/150‑1 (MQ=255)
                                                                               tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATATCGATGTAAAGCAATAAGTTTAAAATGAGCCCCAAACAAGGCGCGCACAAACAAAAATACCa                                                                       >  2:3130989/1‑150 (MQ=255)
                                                                                       gTATTTTGCGTTCGACCTTTGCTTATAAATAACAGTCTTTTTTTTACTTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAATTTGCGCGCAAATCTTGGCGCGCACAAGCTATAATACCAACCTAAtt                                                               <  1:1571277/150‑1 (MQ=255)
                                                                                          tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAATAAATAAAGAAAAAATCCATGTAAAGCACTAAGTTTAAATTTAGCGCCAAAAAAGGCGCGCACAAGCTATAATACCAAACAAAATtct                                                            >  2:3473499/1‑150 (MQ=255)
                                                                                           ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  2:4393527/150‑1 (MQ=255)
                                                                                           ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  2:315827/150‑1 (MQ=255)
                                                                                           ttttCGTTAGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  2:3185769/150‑1 (MQ=255)
                                                                                                tttcgACCTTTGCTTATACATAACAGACTTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGACAATAGCGATGTAAAGCACTACGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                      <  1:1953472/149‑1 (MQ=255)
                                                                                                 ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATAAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCAAGGCGCGCACAAGCAATAATACCAACCTAATTTatcctccac                                                     >  2:4458986/1‑148 (MQ=255)
                                                                                                 ttcgACATTTTCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGAAAATCGCGATGTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcctc                                                     <  1:1345862/150‑1 (MQ=255)
                                                                                                  tcgACCTGTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                    <  2:3290278/150‑1 (MQ=255)
                                                                                                          tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  1:2778786/150‑1 (MQ=255)
                                                                                                           gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                           <  2:4707144/150‑1 (MQ=255)
                                                                                                              tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAATCGAAGTAAATCACTAAGTTTAAATTTCGCGCCAAACATTGCGCGCACAAGCTATAAAAACAACCAAAATTATACACCTCATAGAGTGACtat                                        >  2:2267666/1‑150 (MQ=255)
                                                                                                                  cATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAACAAATAAAGACAAAAGCGATTTAAAGCACAAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCAATAAAACCAACCTAATTTATCCTCCGCTTAGAGTGACTatataa                                    >  2:1618265/1‑148 (MQ=255)
                                                                                                                     aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                 <  1:4528443/150‑1 (MQ=255)
                                                                                                                       cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                               <  2:2206531/150‑1 (MQ=255)
                                                                                                                        aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAACCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                              >  2:3820302/1‑150 (MQ=255)
                                                                                                                         gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                             <  2:1363309/150‑1 (MQ=255)
                                                                                                                            gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          >  2:1340095/1‑150 (MQ=255)
                                                                                                                            gTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          <  1:1798349/150‑1 (MQ=255)
                                                                                                                               tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:1817285/150‑1 (MQ=255)
                                                                                                                               tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:4374968/150‑1 (MQ=255)
                                                                                                                                 tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:3764973/150‑1 (MQ=255)
                                                                                                                                 tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  2:227925/150‑1 (MQ=255)
                                                                                                                                    aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                  >  1:5090679/1‑150 (MQ=255)
                                                                                                                                      tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  2:1431760/150‑1 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  2:2497846/150‑1 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  2:4277652/1‑150 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  2:4169700/150‑1 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  2:1716734/1‑150 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  2:3925180/150‑1 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  1:767793/1‑150 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  1:4376682/1‑150 (MQ=255)
                                                                                                                                         ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             >  1:206399/1‑150 (MQ=255)
                                                                                                                                           aaaaTAATTCAAAGAAATCAAGACAATAGCGAGGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCGGATAGAAATGCGGTATTGACTa           >  2:3004893/1‑150 (MQ=255)
                                                                                                                                            aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  1:4344603/1‑150 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         <  1:1546674/150‑1 (MQ=255)
                                                                                                                                             aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTATTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         <  2:3605408/150‑1 (MQ=255)
                                                                                                                                                aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  1:1108027/143‑1 (MQ=255)
                                                                                                                                                aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  1:5020727/150‑1 (MQ=255)
                                                                                                                                                  ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  2:3874304/150‑1 (MQ=255)
                                                                                                                                                  ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    >  2:2407135/1‑150 (MQ=255)
                                                                                                                                                    cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  >  1:1026064/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTA‑TTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470280‑2470576

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: ‑