Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 25.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,4270CA25.7% 387.9 / 45.5 180Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (35/93);  new base A (47/5);  total (82/98)
Fisher's exact test for biased strand distribution p-value = 2.25e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470281‑2470576
                                                                                                                                                  |                                                                                                                                                     
tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATaaaa                                                                                                                                                    >  1:1019421/1‑150 (MQ=255)
tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:3786678/150‑1 (MQ=255)
tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:225762/150‑1 (MQ=255)
 agtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGAAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAg                                                                                                                                                   >  2:2854372/1‑150 (MQ=255)
     tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCaaataaat                                                                                                                                               >  1:2314251/1‑150 (MQ=255)
       tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATATAATAATTaaaaaaaaaaa                                                                                                                                             >  1:1021573/1‑147 (MQ=255)
       tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCa                                                                                                                                             >  2:70985/1‑150 (MQ=255)
       tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAtaa                                                                                                                                             >  1:1669028/1‑148 (MQ=255)
         tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                                                                                                           <  1:1313348/150‑1 (MQ=255)
         tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                                                                                                           <  1:170026/150‑1 (MQ=255)
         tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                                                                                                           <  1:170044/150‑1 (MQ=255)
          ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                          <  1:2880717/150‑1 (MQ=255)
              ggAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTTAAATAATTAAAAGAAATaaaaaaaaa                                                                                                                                      >  1:3681847/1‑144 (MQ=255)
               gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                                                                                                     <  1:500838/150‑1 (MQ=255)
                 atatCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                                                                                                   <  1:2138329/150‑1 (MQ=255)
                   atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                 <  1:3017648/150‑1 (MQ=255)
                    tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  2:216474/150‑1 (MQ=255)
                      tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTAAAATAAATaaataaaatatagattt                                                                                                                              >  2:3788232/1‑143 (MQ=255)
                        aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAATACAAAATCGATGTaa                                                                                                                            >  1:623343/1‑150 (MQ=255)
                        aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAATAAAGAAAATAACGATGTaa                                                                                                                            >  2:3096855/1‑150 (MQ=255)
                         aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATGTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           <  1:345268/150‑1 (MQ=255)
                          ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                          <  2:1280119/150‑1 (MQ=255)
                            ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  1:2663397/150‑1 (MQ=255)
                            ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  2:2027181/150‑1 (MQ=255)
                            ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAATTAAAAGAAATCAATAAAAAAGCGAAGTAAATCa                                                                                                                        >  2:2064568/1‑150 (MQ=255)
                            ttATGTGTCAATTTTCGGCGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATATTTCGTTCGACCATTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                                                                                                        <  2:751280/150‑1 (MQ=255)
                              aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                                                                                                      <  2:813995/150‑1 (MQ=255)
                               tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAATAAAAAAAGACAAAAGCGATGTAAAGCAcga                                                                                                                     >  2:1267269/1‑148 (MQ=255)
                                 tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACAAAg                                                                                                                   >  2:3889715/1‑150 (MQ=255)
                                   tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAAAAAATAAAGACAAAATCGATGTAAATCACAAAGtt                                                                                                                 >  1:3315652/1‑150 (MQ=255)
                                   tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAAGCACAAAGtt                                                                                                                 >  2:3920410/1‑150 (MQ=255)
                                     aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATATATTTGGAATCTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTAAAAAATAATTAAAAGAAAAAAAGACAAAAACAATGTAAAGCAAAAAATTTa                                                                                                               >  2:1458449/1‑150 (MQ=255)
                                       ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  2:3601154/150‑1 (MQ=255)
                                            ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                                                                                                        <  1:365802/150‑1 (MQ=255)
                                             gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAATAAAGAAAATAGCGATGTTAAGCACGAAGTTTAAAtttagc                                                                                                       >  1:822133/1‑146 (MQ=255)
                                               gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAATTAAAATAAATAAAAACAAAAGCGATTTAAAGCACTAATTTTAAAATgcgcgc                                                                                                     >  1:665424/1‑150 (MQ=255)
                                                aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAAAAAAATAAAAAAAGACAAAAGAGAAGTAAAGCACTAAGATTAAATTTAGCGcc                                                                                                    >  1:1422585/1‑150 (MQ=255)
                                                  ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                                                                                                  >  2:1257711/1‑150 (MQ=255)
                                                  ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATATCGATGTAAAGCACTAAGTTTAAATTTAGCGCCaa                                                                                                  >  2:3088472/1‑150 (MQ=255)
                                                   gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAAAAAATAAAAAAAAAAAACACAAAAACGTTGTAAATCACAAATTTTAAATTTCGGGcccaa                                                                                                 >  1:750087/1‑147 (MQ=255)
                                                   gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAAAAAAAATAGATGTAAAGCACTAAGTTTAAATTTAGCGCCacg                                                                                                 >  1:2031184/1‑148 (MQ=255)
                                                   gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAATAAATCAAGACAAAAGCGATGTAAATCACTAAGTTTAAATTTGGCGCCAAt                                                                                                 >  2:880988/1‑150 (MQ=255)
                                                    atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAACAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATc                                                                                                >  2:948856/1‑150 (MQ=255)
                                                    atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAAAAAAAAAAAAACAAAAAAAAAACCATGTAAAACACAAAATTTAAAATTAGCCCCAAta                                                                                                >  2:1696095/1‑149 (MQ=255)
                                                     tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                                                                                               <  2:743505/150‑1 (MQ=255)
                                                     tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                                                                                               <  2:413752/150‑1 (MQ=255)
                                                     tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                                                                                               <  2:743587/150‑1 (MQ=255)
                                                       taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAATACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                             >  2:1942788/1‑150 (MQ=255)
                                                       taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAGCGATGGAAAGCACTAAGTTTAAATTTCGCGCCAATCaag                                                                                             >  2:1537372/1‑148 (MQ=255)
                                                       taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATAAATAAAAAATAGATGTAAATCACAAATTTTAAATTTAGCGCCAAACaag                                                                                             >  2:3515/1‑148 (MQ=255)
                                                           gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACATTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                         <  2:2283710/150‑1 (MQ=255)
                                                           gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc                                                                                         <  2:2627092/150‑1 (MQ=255)
                                                                aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAAAAAAAAATAGATGTAAAGCACAAAGTTTAAATTTCGCGCCCAACAAGGCGCGccca                                                                                    >  2:1782256/1‑150 (MQ=255)
                                                                   caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAAAAAAAAAAAAAACGAAGTAAATCACAAACTTTTAATTTAGCGCCAAACAAGGCGCGCCCaacc                                                                                 >  1:2960524/1‑148 (MQ=255)
                                                                    aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAAAAAAGAAAAAATCGATGTAAAGCACAAAGTTTAAAATTCTCGCCCAGCATGGCGCGCACaaaca                                                                                >  1:2331273/1‑147 (MQ=255)
                                                                             tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAATACAATATCGGTGTAAAGCACTAAGTTTATATTGGGCGCCAATCATGGCGCGCCCAAGCGATAATAACa                                                                       >  2:57971/1‑150 (MQ=255)
                                                                              ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCaa                                                                      <  2:545122/150‑1 (MQ=255)
                                                                               tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                     <  2:2351708/150‑1 (MQ=255)
                                                                                     ttttttttCTTTCTACTTTGGATTATAATTAACTGTCTTTTTTTTATTTTATAAAATAATTAAAATAAATCATCACATTAGGTTTTTAATGAAATATTTTTAATTTGCGCGCCAATCTTGGCGCGCACAAGCTATAATACCAACCTAAtt                                                               <  1:1552198/147‑1 (MQ=255)
                                                                                     ttttttttCGTTCGAACTTTGATTATAAATAACTGTTTTTTTTTTATTTTATAAAATAATTAAAATAAATAATGACAATTGCTCTTTAACGCAATATTTTTAAATTGGGCGCCAATAATGGCGCGCACAAGCTATAATACCAACCTAAtt                                                               <  1:1263267/147‑1 (MQ=255)
                                                                                     gTATTTTTCGTTCGACATTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGGGCACAAGCTATAATACCAACCTAAtt                                                               <  1:733119/150‑1 (MQ=255)
                                                                                         ttttCTTTCTACATTTGCTTATAAATAACTGTCTTTTTTTTATTTTATAAAATAATTAAAATAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  1:1458449/150‑1 (MQ=255)
                                                                                         ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  2:2645135/150‑1 (MQ=255)
                                                                                          tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATATATAAAATAATTAAAAGAAAAAAAGAAAAAAGCGATTTAAACCACTAAGTTTAAATTGCGCGCCCATCATGGCGCGCCCAAGCGATAATAACAACCAAATTTctcc                                                          >  2:3367936/1‑150 (MQ=255)
                                                                                            tcgttcgACCTTTGCTTATAAATAACTGTTTTTTTTTTAATTTATAAAATAATAAAAAGAAATCAAGACAATAGCGATGTAAAGCACTATGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctc                                                        <  1:790939/150‑1 (MQ=255)
                                                                                             cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                       <  2:1964992/150‑1 (MQ=255)
                                                                                              gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCGATGAAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcct                                                      <  1:974901/150‑1 (MQ=255)
                                                                                                tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                    <  1:2468666/150‑1 (MQ=255)
                                                                                                    ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                <  1:235607/150‑1 (MQ=255)
                                                                                                      tttGCTTATACATAACTGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGt                                              <  2:1422585/150‑1 (MQ=255)
                                                                                                        tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                            <  1:3375820/150‑1 (MQ=255)
                                                                                                              taCATAACAGTCGTTTTTTTAAATTATAAAAAAATTAAAAGAAAACAAAACAAAAACGATGTAAAGCACAAAGTTTAAAATTAGCGCCAAACAATGCGCGCACAATCAATAAAAACAACCAAAATTATACTCCACGTAGAGTGACaataa                                      >  2:337310/1‑149 (MQ=255)
                                                                                                                  tAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGAAAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCCATCATGGCGCGCACCAGCTAAAATACCAACCTAAATTGTACTCCTCTTAGAGTGACTATATCTcc                                  >  2:2416718/1‑150 (MQ=255)
                                                                                                                  tAACACTCGTTTTTTAAATTTATAAAAAAATACAAATAAAACAAGACAAAAGCGATGTAAAGCAATAAGTTTAAATTGCCCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                  <  1:1883837/150‑1 (MQ=255)
                                                                                                                   aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                 <  1:3196797/150‑1 (MQ=255)
                                                                                                                     cAGTCTTTTTTTTATTTTATAAAATAATTCATAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTACTCCTCTTAGAGTGACTATATCTCCTGa                               <  1:2735732/150‑1 (MQ=255)
                                                                                                                     cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                               <  1:3054623/150‑1 (MQ=255)
                                                                                                                       gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                             <  1:3920410/150‑1 (MQ=255)
                                                                                                                             tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:2736324/150‑1 (MQ=255)
                                                                                                                             tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:3414510/150‑1 (MQ=255)
                                                                                                                             tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:1813227/150‑1 (MQ=255)
                                                                                                                             tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       <  2:1197637/150‑1 (MQ=255)
                                                                                                                              ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  1:1180679/150‑1 (MQ=255)
                                                                                                                              ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  1:3234572/150‑1 (MQ=255)
                                                                                                                              ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  1:3234643/150‑1 (MQ=255)
                                                                                                                               tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:3220781/150‑1 (MQ=255)
                                                                                                                               tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     >  1:1535479/1‑150 (MQ=255)
                                                                                                                               tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     >  1:1426488/1‑150 (MQ=255)
                                                                                                                               tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     >  2:1615264/1‑150 (MQ=255)
                                                                                                                               tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  2:468058/150‑1 (MQ=255)
                                                                                                                                    tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCGTGTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                >  1:1139615/1‑150 (MQ=255)
                                                                                                                                      tataAAATAATTCAAAGAAATCAAGACAATCGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  2:2977627/150‑1 (MQ=255)
                                                                                                                                      tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  1:717316/150‑1 (MQ=255)
                                                                                                                                         aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa           <  1:548361/150‑1 (MQ=255)
                                                                                                                                         aaaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTa           >  1:1617894/1‑150 (MQ=255)
                                                                                                                                            ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa        >  1:3790396/1‑150 (MQ=255)
                                                                                                                                                ttCAAATAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  1:3331290/150‑1 (MQ=255)
                                                                                                                                                ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    >  2:1208393/1‑150 (MQ=255)
                                                                                                                                                ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg    <  2:842229/150‑1 (MQ=255)
                                                                                                                                                  cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  <  2:1021573/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470281‑2470576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: