Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 21.2% 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,4270CA21.2% 888.0 / 79.8 369Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (91/187);  new base A (81/8);  total (174/195)
Fisher's exact test for biased strand distribution p-value = 1.41e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470280‑2470576
                                                                                                                                                   |                                                                                                                                                     
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:2822398/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    >  1:5214161/1‑150 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  2:4114742/150‑1 (MQ=255)
 tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  2:4019648/150‑1 (MQ=255)
   gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGa                                                                                                                                                  <  2:3072536/150‑1 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  1:2606439/150‑1 (MQ=255)
     aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  2:3356531/150‑1 (MQ=255)
      tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAaaataaat                                                                                                                                               >  2:5697384/1‑150 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                             <  2:4628510/150‑1 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                             <  1:3345212/150‑1 (MQ=255)
        tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                             <  2:4960857/150‑1 (MQ=255)
            tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGaa                                                                                                                                         >  2:436794/1‑149 (MQ=255)
               ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAAATAAAAGAAATAAATACAAt                                                                                                                                      >  1:3022985/1‑150 (MQ=255)
                gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATAAAGAGATTTAAGCAAGCATTTATTTGGAATATATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAAGAAATAAAGACaaaa                                                                                                                                     >  2:872975/1‑148 (MQ=255)
                 aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAACAAGAAAAAAg                                                                                                                                    >  2:5482976/1‑150 (MQ=255)
                    atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                 <  2:2973710/150‑1 (MQ=255)
                     tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  2:173905/150‑1 (MQ=255)
                      ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAATATAGATg                                                                                                                               >  1:4857957/1‑150 (MQ=255)
                       tctaTTTTTTGTGTGATTTTTGGGGGAAGGATAGAGAGATTTAAGCAAGCATTTTTTTGGTATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                              <  1:4439460/147‑1 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAAAAAATAAAAAAGCGaatt                                                                                                                              >  1:2651461/1‑147 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATAAAGACAATATCGATGt                                                                                                                              >  1:4452411/1‑150 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATCAATAAAAAATCGATGt                                                                                                                              >  2:1678233/1‑150 (MQ=255)
                       tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAAAAAGAAAATATAGATGt                                                                                                                              >  1:4059605/1‑150 (MQ=255)
                        cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                                                                                             <  2:1563062/150‑1 (MQ=255)
                         aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAAGAAAAAAGCGATGTaa                                                                                                                            >  1:2357788/1‑150 (MQ=255)
                          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           <  2:413947/150‑1 (MQ=255)
                          aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           <  2:4315303/150‑1 (MQ=255)
                                 gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                                                                                                    <  2:2833500/150‑1 (MQ=255)
                                  tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                                                                                                   <  1:4641279/150‑1 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGACAATAGCGATGTAAAGAACTAAGtt                                                                                                                 >  1:3776831/1‑150 (MQ=255)
                                    tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAATCAATAAAATAGCGATGTAAAGCACAAAGtt                                                                                                                 >  1:4847931/1‑150 (MQ=255)
                                     cATTTTTCGGGGAATGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                                                                                                <  1:5574634/150‑1 (MQ=255)
                                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAAAAAATAAAATAAAAAAAAAAAAAAGCGATGTAAATAACAAAATTTa                                                                                                               >  1:1979474/1‑150 (MQ=255)
                                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAAAAATTAAAAGAAATAAAGACAAAAGCGATGGAAAAGACGAAGTTTa                                                                                                               >  1:1546160/1‑150 (MQ=255)
                                      aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATAAATACAAAATCGATGTAAATCACAAAGtttt                                                                                                               >  2:3630293/1‑149 (MQ=255)
                                       tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                                                                                              <  1:265523/150‑1 (MQ=255)
                                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGTAAAGCAAGAAGTTTAAAtt                                                                                                           >  1:4326367/1‑150 (MQ=255)
                                          ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAAAAGCGGTGTTAATCAAAAAGTTTTAAtt                                                                                                           >  1:1000067/1‑150 (MQ=255)
                                             ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                        <  1:1384729/150‑1 (MQ=255)
                                                     atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAATAATTTAAAGAAATCAAGTCAATATCGATGTTAAGGACTAAGTTTTAATTTGGCGCCAATc                                                                                                >  2:5130827/1‑150 (MQ=255)
                                                     atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATAAAAATAAATCAAAAAAAAAGAGATTTTAATAACAAAGTTTAAATATAGCGCCaaaa                                                                                                >  1:4688702/1‑148 (MQ=255)
                                                         agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAATAAATAAAAAAAAAAAAAAAAAAAGCGATGTAAATAAAAAAATTTAAATTTAGCGCCaaaaaagg                                                                                            >  1:1794807/1‑144 (MQ=255)
                                                         agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAAGAAATAAAGACAAAAGCGATTTAAATCACAAAGTTTAAATTTCGCGCCAAGCATgg                                                                                            >  2:2737256/1‑150 (MQ=255)
                                                         agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAATAAATAAAGACAAAATCGATTTAAATCACAAAGTTTAAATTTAGCGCCAAACAttg                                                                                            >  2:2270255/1‑148 (MQ=255)
                                                         agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATTAAATAATTAAAATAAATCAATAAAAAAACGATTTAAAGGAAGAAGTTTAAATTTAGCGCCaaaaaatg                                                                                            >  1:57573/1‑144 (MQ=255)
                                                           agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAAAAACGATGTAAATCACTAAGTTTAAATTTAGCGCCCATCAAGgcg                                                                                          >  2:4604430/1‑150 (MQ=255)
                                                           agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTAATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATATCGATGTAAAGAACTAACTTTAAATTTCGCGCCAATCATTgcg                                                                                          >  2:5273260/1‑150 (MQ=255)
                                                             aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATTAATAAAATAAATAAAAAACAAATCGCTTTAAAGCACGAAGTTTAAAATTCGCGCCCAACCTTgcgcg                                                                                        >  2:4876545/1‑150 (MQ=255)
                                                               ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgcc                                                                                      >  2:3127951/1‑149 (MQ=255)
                                                                tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaa                                                                                     <  2:5442118/150‑2 (MQ=255)
                                                                   gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAATCACTAAGTTTAAATTTCGCGCCAATAATGGCGCGCACAAg                                                                                  >  2:907757/1‑150 (MQ=255)
                                                                       gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAAGAAAAAAAGAAAAAAGAGATGTTAAGAAAAAATTTTAAATTTCACGCCAAAAATTGCGCGAACAAGAtaa                                                                              >  2:264801/1‑149 (MQ=255)
                                                                           ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                          <  2:1162874/150‑1 (MQ=255)
                                                                              tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAAAAAAGAAAAAAAAGATGTAAAACAAAAAGTTTTAATTTAGCGCCCAACAAGGGGCGCCCCAGCCATAATAACa                                                                       >  1:5802835/1‑150 (MQ=255)
                                                                               ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAAACCaa                                                                      >  2:1867796/1‑150 (MQ=255)
                                                                                     tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGAGATGTAAATCAATAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATACCaaacaaaa                                                                >  2:570099/1‑144 (MQ=255)
                                                                                      gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATATCGATGTAAATCACAAAGTTTAAATTTCGCGCCAATCATTGCGCGCACAAGCTATAAAACCAACCTAAtt                                                               >  2:3858632/1‑150 (MQ=255)
                                                                                         tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAAAAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAAAAATTGCGCGCACAAGCTATAATAACCACCAAAattat                                                            >  2:2109525/1‑148 (MQ=255)
                                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGTTGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  1:1835352/150‑1 (MQ=255)
                                                                                          ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  2:2486316/150‑1 (MQ=255)
                                                                                           tttCGTTCGACCTTTGCTTATAAATAACAGTCTTTTTTTTAATTTATAAAATATTTCAAATAAATCAATACAATAGCTATGTAAAGCACTAAGTTTAATTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                          <  1:1899888/150‑1 (MQ=255)
                                                                                              cgttcgACCTTTGCTTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  1:2266823/150‑1 (MQ=255)
                                                                                               gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTcacctcct                                                      >  2:5738260/1‑150 (MQ=255)
                                                                                               gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATAAAGAAAAAATCGATTTAAAGCAAAAAGTTTAAAATTAGCGCCAATAATGGCGCGCACAAGCTATAATAACAAACTAATATatactcca                                                      >  2:4178689/1‑149 (MQ=255)
                                                                                                   gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAACACAATAGCGATGTAAAGCACAAAGTTTAAATTGCGCCCCAATAATGGCGCGCACAAGCTATAATAACAAACAAATTTATCCACCGCTTa                                                  >  2:1340647/1‑150 (MQ=255)
                                                                                                     ccTTTGATTATACATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                <  1:4229139/150‑1 (MQ=255)
                                                                                                       tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACAAAGTTTAAATTGCGCGCCAAGCAAGGCGCGCACAAGCTATAATACCAACCAAATTTCTCCTCCTCTTAGAGt                                              >  2:1972942/1‑150 (MQ=255)
                                                                                                         tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTgn                                            <  1:3595672/150‑2 (MQ=255)
                                                                                                         tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  1:2462122/150‑1 (MQ=255)
                                                                                                         tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  2:2939412/150‑1 (MQ=255)
                                                                                                         tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  1:4568562/150‑1 (MQ=255)
                                                                                                         tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  1:1738905/150‑1 (MQ=255)
                                                                                                                   tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                  <  1:3681627/150‑1 (MQ=255)
                                                                                                                   tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                  <  1:2391229/150‑1 (MQ=255)
                                                                                                                       aGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                              <  1:5678394/150‑1 (MQ=255)
                                                                                                                          cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGa                           <  1:1608419/150‑1 (MQ=255)
                                                                                                                           gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          >  1:5468333/1‑150 (MQ=255)
                                                                                                                           gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTACAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          >  2:3173053/1‑150 (MQ=255)
                                                                                                                               ttttAATGTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  2:648222/150‑1 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:671406/150‑1 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:1281186/150‑1 (MQ=255)
                                                                                                                                tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     >  1:291291/1‑150 (MQ=255)
                                                                                                                                  tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                   >  1:523653/1‑150 (MQ=255)
                                                                                                                                  tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                   <  2:761110/150‑1 (MQ=255)
                                                                                                                                   aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                  <  1:4011190/150‑1 (MQ=255)
                                                                                                                                    aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAt                 >  1:500463/1‑150 (MQ=255)
                                                                                                                                    aTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACCATATCTCCTGATAGAATTGCGGTAt                 >  2:4801294/1‑150 (MQ=255)
                                                                                                                                     tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  2:373105/150‑1 (MQ=255)
                                                                                                                                        ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  1:3246109/150‑1 (MQ=255)
                                                                                                                                           aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  2:690593/1‑150 (MQ=255)
                                                                                                                                           aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt          >  1:4270678/1‑150 (MQ=255)
                                                                                                                                            aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         <  1:885781/150‑1 (MQ=255)
                                                                                                                                             ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        <  2:3018858/150‑1 (MQ=255)
                                                                                                                                               aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  1:3191717/150‑1 (MQ=255)
                                                                                                                                               aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa      <  2:4934549/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  1:3135806/150‑1 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    >  2:106021/1‑150 (MQ=255)
                                                                                                                                                 ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    >  1:113869/1‑150 (MQ=255)
                                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  <  1:1660132/150‑1 (MQ=255)
                                                                                                                                                   cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  <  1:5087183/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  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: