Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 18.0% 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,4270CA18.0% 562.0 / 39.9 220Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (54/118);  new base A (41/5);  total (97/123)
Fisher's exact test for biased strand distribution p-value = 7.94e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470278‑2470576
                                                                                                                                                     |                                                                                                                                                     
ttGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTc                                                                                                                                                       <  1:2346901/150‑1 (MQ=255)
  gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaa                                                                                                                                                     <  1:3744106/150‑1 (MQ=255)
   tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    >  2:3527875/1‑150 (MQ=255)
   tagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaaa                                                                                                                                                    <  1:5517871/150‑1 (MQ=255)
     gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAAAAATTAAAAGa                                                                                                                                                  >  1:4913947/1‑150 (MQ=255)
       aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  2:5710408/150‑1 (MQ=255)
          tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAtaa                                                                                                                                             >  1:3809651/1‑148 (MQ=255)
          tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAtaa                                                                                                                                             >  1:1918400/1‑148 (MQ=255)
           ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAATAAATCaa                                                                                                                                            >  2:4601211/1‑150 (MQ=255)
                  gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                                                                                                     <  1:522868/150‑1 (MQ=255)
                      atCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                                                                                                 <  2:5436087/150‑1 (MQ=255)
                       tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTGTAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  2:1603547/150‑1 (MQ=255)
                       actcAATTTTATGTGACAATTTTCGTGTAAGGATATAGACATAAAACCAACCATTTAATTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGATTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                                <  1:3850576/149‑1 (MQ=255)
                        ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAAGAGATg                                                                                                                               >  2:4187610/1‑150 (MQ=255)
                         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                                                                                              >  2:3398937/1‑150 (MQ=255)
                         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAATAGCGAttt                                                                                                                              >  1:3581984/1‑148 (MQ=255)
                         tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAAAAAAAAAAAAAACGATGt                                                                                                                              >  1:17031/1‑150 (MQ=255)
                            aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                                                                                           <  2:5586439/150‑1 (MQ=255)
                            aTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGAAAATAGAGATTTaaa                                                                                                                           >  2:4147539/1‑150 (MQ=255)
                             ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAGAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                          <  2:5125767/150‑1 (MQ=255)
                             ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                          <  1:5924895/150‑1 (MQ=255)
                             ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAATACAATAGCGATGTaaat                                                                                                                          >  2:1324566/1‑149 (MQ=255)
                                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGtt                                                                                                                 >  2:5839666/1‑150 (MQ=255)
                                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAATAGAGATGTAAAGCACTAAGtt                                                                                                                 >  1:5495283/1‑150 (MQ=255)
                                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATAGCGATGTAAATCACTAAGtt                                                                                                                 >  2:3374058/1‑150 (MQ=255)
                                      tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAAATCAAATAAATAAAGAAAAAAGCGATGTAAAGCACTAAGtt                                                                                                                 >  1:2744192/1‑150 (MQ=255)
                                        aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATTAAATAATTCAAATAAATAAAGACAAAAGCGGTGTAAAGGACTAAGtttt                                                                                                               >  1:5680043/1‑149 (MQ=255)
                                        aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAATAAATAAAGAAAAAATCGAAGTAAAGCACTAAGTTTa                                                                                                               >  1:2012144/1‑150 (MQ=255)
                                        aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAATAAAAACAAAAGCGATTTAAAGAACGAAGTTTa                                                                                                               >  1:2901563/1‑150 (MQ=255)
                                         tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAATACAAAAGCGATTTAAAGAACTAAGTTTaa                                                                                                              >  2:5185434/1‑150 (MQ=255)
                                          ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                             <  2:2789817/150‑1 (MQ=255)
                                               ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTTAAGCACGAAGTTTAAATTgag                                                                                                        >  2:4818846/1‑148 (MQ=255)
                                               ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAAATAAAAGAAATAAAGAAAATAGCGATGTAAAGCAAAAAGGTTAAATTgag                                                                                                        >  1:765941/1‑148 (MQ=255)
                                                   aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGAAAAAAGCGATGTAAAGCACTAAGTTTAAATTTCGCGcc                                                                                                    >  2:3222325/1‑150 (MQ=255)
                                                   aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAAGAAATAAAGACAAAATCGATGTAAAGCACAAAGTTTAAAATGCGCGcc                                                                                                    >  2:1818549/1‑150 (MQ=255)
                                                     ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAAAAATAAAAAATCGATGTAAATCAATAATTTTAAATTTAGCGccca                                                                                                  >  1:1506675/1‑148 (MQ=255)
                                                      gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                                                                                                 <  1:2211856/150‑1 (MQ=255)
                                                       atatAGAGATTTAAGCAAGCTTTTATTTGGAATGTATTTTTCGTTCGTCCTTTGCTTATACATATCAGTCGTTTTTTTAATTTATAAAATTATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCatga                                                                                                <  2:5063209/150‑4 (MQ=255)
                                                       atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAATAATTAAAAGAAATAAAGAAAATAGCCATGTAAAGCACAAAGTTTAAATTTCGCGCCaaac                                                                                                >  1:4525380/1‑148 (MQ=255)
                                                        tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATAAATACAAAATCGATGTTAAGCACTAATTTTTAATTTAGCGCCaaaaa                                                                                               >  1:992281/1‑147 (MQ=255)
                                                         ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                                                                                              <  2:5461427/150‑1 (MQ=255)
                                                         ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAAAAAAAAACAAAAAATCCAAGTAAAACAAAAAGTTTAAAATTAGCGCCaaaaaa                                                                                              >  2:526783/1‑146 (MQ=255)
                                                           agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAATAATTGAAAAAAAAAAAAACAAAAGCGATGTAAAGCACAAAGTTTGAATTTCGCGCCACTCAAgg                                                                                            >  1:5192393/1‑150 (MQ=255)
                                                            gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATAAATACAAAATCGCTGTAAATAACAAATTTTAATTTTCGCGCCGATCATGgc                                                                                           >  1:1000988/1‑150 (MQ=255)
                                                            gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAATAAATCAAGAAAAAAGAGAGTTAAAGCACAAATTTTAAATTTGGCGCCCAACATGgc                                                                                           >  1:1781263/1‑150 (MQ=255)
                                                                   aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATCAAAGAAAAAAAGACAATATCGAAGTAAAGCACAAAGTTTAAATTGCGCGCCAAACATGGCGCGccca                                                                                    >  2:227961/1‑150 (MQ=255)
                                                                   aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaat                                                                                    <  2:4334887/150‑3 (MQ=255)
                                                                      caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATAGCGATGTAAAGGACTAAGTTTAAATTGCGCGCCAATCATTGCGCGCCCAAGc                                                                                 >  2:1177049/1‑150 (MQ=255)
                                                                       aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt                                                                                <  2:1065626/150‑1 (MQ=255)
                                                                        agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatttaaa                                                                               <  2:2181780/150‑8 (MQ=255)
                                                                        agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaatttaaa                                                                               <  2:1080422/150‑8 (MQ=255)
                                                                         gcaTTTATTTGGAATGTATTTTTCGTTCGACCTGTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtat                                                                              <  2:1969993/150‑1 (MQ=255)
                                                                         gaatttttttGTATTTTTTTTTTCGTTCGACATTTGATTATACAAAACTTTCTTTTTTTTAATTTATAAAATAATTAATATAAATAAAGAAATTAGCTATTTAAAGCACTATGTTTAAATTGCGCGCAAATCTTGGCGCGCACAAGCtat                                                                              <  1:3996338/148‑1 (MQ=255)
                                                                           atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAATAAAAAAAAAAAAAGACAAAAAACATGTTAAGAACAAAGTTTAAATTTAGCGCCACAGAAGGCGCGCACaaccaaaaa                                                                            >  1:281203/1‑143 (MQ=255)
                                                                                    gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATATCGATGTAAAGCACTAAGTTTAAATTTGGCGCCAAAAAATGCGCCCCCAAGCTATAATAACaaaca                                                                   >  2:1581754/1‑147 (MQ=255)
                                                                                         tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                              <  2:4512859/150‑1 (MQ=255)
                                                                                          aTTTTTCGTGCGACCTTTGCTTATACATATCCGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACTATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtc                                                             <  2:1399552/150‑1 (MQ=255)
                                                                                           tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtct                                                            <  1:3454055/150‑1 (MQ=255)
                                                                                            ttttCGTTCGACCTTTGCTTATACATAACCGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATTTAAAGCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                           <  1:317269/150‑1 (MQ=255)
                                                                                             tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                          <  2:5664463/150‑1 (MQ=255)
                                                                                                 gttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATTATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCACCCTAATTTctcctcat                                                      <  2:2503682/150‑3 (MQ=255)
                                                                                                  ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAACAACCTAAATTcaactccgc                                                     >  2:5515928/1‑148 (MQ=255)
                                                                                                  ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATCAATAAAAAAGCGATGTAAATCACTAAGTTTAAATTGCGCGCCAAACATGGCGCGCACAAGCTATAATACCCACCTAATTTaaccacctc                                                     >  2:4924935/1‑150 (MQ=255)
                                                                                                     gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCttc                                                  >  2:4479540/1‑149 (MQ=255)
                                                                                                      aCCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCAAATTTCTCCTCCTCGTag                                                 >  2:5918429/1‑150 (MQ=255)
                                                                                                       ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                                <  1:3781771/150‑1 (MQ=255)
                                                                                                         tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATCAAGACAAAAGCGATGTAAATCACAAATTTTAAATTTCGCGCCAATCATGGCGCGCCCAAGCTATAATACCAACCAAAATTATCCTCCTCTTAGAGt                                              >  2:1890178/1‑150 (MQ=255)
                                                                                                            gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATCGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                           <  1:4876714/150‑1 (MQ=255)
                                                                                                            gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                           <  1:559900/150‑1 (MQ=255)
                                                                                                              ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAAAACCAACCTAATTTCTCCTCCGCTTAGAGTGAcaa                                         >  2:3305465/1‑148 (MQ=255)
                                                                                                              ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                         <  1:1461968/150‑1 (MQ=255)
                                                                                                                ataCATAACAGTCGTGTTTTTAATTTATAAAATAAATCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCCCAAGCTATAAAAACAACCAAATTTCTCCTCCTATTAGAGTGACtata                                       >  2:697787/1‑150 (MQ=255)
                                                                                                                   cATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAAAAAAGACAAAAGCCATGTAAAGCAAAAAGTTTAAATTTCGCGCCAAAAAAGGCGCGCACAAGCAAAAAAAACAACCAAAAATCTCCACCACATAGAGTGACTatatat                                    >  2:5240330/1‑148 (MQ=255)
                                                                                                                     tAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTcc                                  <  2:4254951/150‑1 (MQ=255)
                                                                                                                           tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            <  1:1752780/150‑1 (MQ=255)
                                                                                                                           tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                            >  2:4783775/1‑150 (MQ=255)
                                                                                                                             gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          <  2:2901563/150‑1 (MQ=255)
                                                                                                                             gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          <  1:660641/150‑1 (MQ=255)
                                                                                                                             gTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                          <  2:4880652/150‑1 (MQ=255)
                                                                                                                              tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                         >  2:240757/1‑150 (MQ=255)
                                                                                                                                tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTg                       >  1:4832607/1‑150 (MQ=255)
                                                                                                                                 ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                      <  1:2148810/150‑1 (MQ=255)
                                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  1:1173382/150‑1 (MQ=255)
                                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                     <  2:3973227/150‑1 (MQ=255)
                                                                                                                                  tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAGGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGcn                     <  1:5185713/150‑2 (MQ=255)
                                                                                                                                   ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                    >  2:3632823/1‑150 (MQ=255)
                                                                                                                                   ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                    <  2:2535715/150‑1 (MQ=255)
                                                                                                                                   ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                    >  1:3774539/1‑150 (MQ=255)
                                                                                                                                   ttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCgg                    <  1:3337886/150‑1 (MQ=255)
                                                                                                                                    tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                   >  2:304388/1‑150 (MQ=255)
                                                                                                                                     aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                  <  1:1324566/150‑1 (MQ=255)
                                                                                                                                       tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt                <  1:301283/150‑1 (MQ=255)
                                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  2:1498119/150‑1 (MQ=255)
                                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  2:3867970/150‑1 (MQ=255)
                                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              >  1:497531/1‑150 (MQ=255)
                                                                                                                                         tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  1:857514/150‑1 (MQ=255)
                                                                                                                                          ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc             <  2:2375385/150‑1 (MQ=255)
                                                                                                                                           taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt            <  1:369488/150‑1 (MQ=255)
                                                                                                                                              aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACtctc         >  1:4093152/1‑150 (MQ=255)
                                                                                                                                              aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc         >  1:728205/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470278‑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: