Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,937 T→G 9.3% intergenic (+22/+93) appA → / ← etk phosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,9370TG 9.3% 209.9 / 14.3 70intergenic (+22/+93)appA/etkphosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base G (3/4);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-02

GCCGCCAGGTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAA  >  NC_000913/1041663‑1042202
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                         
gccgccAGTTGGTGTATAGGTGTTTGAAGCCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAAGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGt                                                                                                                                                                                                                                                                           <  2:277838/275‑1 (MQ=255)
        gtggtgAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                                   <  2:52631/275‑1 (MQ=255)
                           aCGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAc                                                                                                                                                                                                                                                 <  2:401866/274‑1 (MQ=255)
                               tGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACTGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACg                                                                                                                                                                                                                                            <  2:115893/275‑1 (MQ=255)
                                ggCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAcc                                                                                                                                                                                                                                                                       >  4:348195/1‑247 (MQ=255)
                                    tCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGGGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                                   <  2:348196/247‑1 (MQ=255)
                                    tCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaagaa                                                                                                                                                                                                                                       <  2:193911/275‑1 (MQ=255)
                                    tCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaagaa                                                                                                                                                                                                                                       <  2:450782/275‑1 (MQ=255)
                                          aaGCGATAACAGCCAGTGGATTCAGGTTTTGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCt                                                                                                                                                                                                                                 <  2:424371/275‑1 (MQ=255)
                                           aGCGATAACGGCCGGTGGATTCAGGGTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTa                                                                                                                                                                                                                                <  4:136879/275‑1 (MQ=255)
                                                    atccAGTGGATTCAGGTTTCGCTGTTCTTCCAGACTTTACAGGAGATGCGTGACAAAACGCCGCTGTCATAAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGAGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAg                                                                                                                                                                                                                        <  2:301184/272‑1 (MQ=255)
                                                                                       tttACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCGGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACttt                                                                                                                                                                                    >  2:403782/1‑275 (MQ=255)
                                                                                         tACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGCGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACttt                                                                                                                                                                                    >  1:90946/1‑273 (MQ=255)
                                                                                           cgtgagATGGGTGCAAAAACGCGGCATGANTTAAAGCCGCCGCCCGGAGAGGTGANACGTGCCCTGGCAGGATGTGAAGAGCGAAATGCGCGGTGCGTGTGTTCGTTGGCAGGTTTTACGACAATCGTGAATGAACNCCGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGt                                                                                                                                                                                <  3:90946/271‑1 (MQ=255)
                                                                                                 aTGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAc                                                                                                                                                                           >  1:74643/1‑274 (MQ=255)
                                                                                                 aTGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCAACTTTCAGTTTTAct                                                                                                                                                                          >  1:10519/1‑275 (MQ=255)
                                                                                                    cggagtAAAACGACGATGGCATTAAATACGCCGCCCGGAGAGGGGAAACTGCCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGAAGGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTAct                                                                                                                                                                          <  4:500663/267‑1 (MQ=255)
                                                                                                                   ggaaagttAAATACGCCGCCCGGGAAGGGGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATTTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACtat                                                                                                                                                        <  4:363336/269‑1 (MQ=255)
                                                                                                                                cgccgcCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGCCCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGCTTTCCTCTTTCTCGGAGTAACTATCACCGTAATAGtt                                                                                                                                           >  2:240844/1‑275 (MQ=255)
                                                                                                                                  ccgccCGGAGAGGTGAAACTGACCCTGGCAGGTTGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTtat                                                                                                                                         <  4:457972/275‑1 (MQ=255)
                                                                                                                                    gccCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTaa                                                                                                                                                            >  2:363337/1‑254 (MQ=255)
                                                                                                                                                 aaaCTGACCCTGGCAGGATGTGGAGGGCGAAATGCGGAGGGCATGTGTTCGTTGGCAGGTTATACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGAAAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAg                                                                                                                          <  2:250882/275‑1 (MQ=255)
                                                                                                                                                 aaaCGGACCCTGGGGAGAAGTGAAGGGGGAAATGGGCAGGGAATGGGTTCGTGGGAGGGTGTTAGGCAAATCGTGAATTAAGCACTCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCATGTACCTGAATGCTCTGAGGCTGATGACAAACTAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAg                                                                                                                          <  4:240843/275‑1 (MQ=255)
                                                                                                                                                      gACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTaa                                                                                                                                                            >  4:297372/1‑236 (MQ=255)
                                                                                                                                                       aCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtct                                                                                                                                                                    <  4:511660/227‑1 (MQ=255)
                                                                                                                                                       aCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtct                                                                                                                                                                    >  2:511661/1‑227 (MQ=255)
                                                                                                                                                          cTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACtat                                                                                                                                                        <  2:297373/236‑1 (MQ=255)
                                                                                                                                                             gCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTg                                                                                                                                                                                                                                                            >  6:849631/1‑133 (MQ=255)
                                                                                                                                                             gCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTg                                                                                                                                                                                                                                                            <  5:849631/133‑1 (MQ=255)
                                                                                                                                                               aGGAAGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAAAGCATAAAAAAGAGCATTCATTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGcgcg                                                                                                             <  2:231200/274‑1 (MQ=255)
                                                                                                                                                                 gATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAAGAACTATAACCGTAAAAGTTAGAGGCGTAACGGTAAGCGGtcatgg                                                                                                                 >  4:231199/1‑263 (MQ=255)
                                                                                                                                                                    gcgAAGAGCGAAATGCGGAGGGCATGGGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCGCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTCCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTaa                                                                                                        <  4:281538/272‑1 (MQ=255)
                                                                                                                                                                        agagCGAAATGCGCAGGGCATGTGTTCGGTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGAATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACGGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATcaca                                                                                                   <  4:411679/275‑1 (MQ=255)
                                                                                                                                                                         gagCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGaa                                                                                                                                                                                                                                                                   >  8:196120/1‑114 (MQ=255)
                                                                                                                                                                         gagCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGaa                                                                                                                                                                                                                                                                   <  7:196120/114‑1 (MQ=255)
                                                                                                                                                                            cGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACAGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTACCTATAACCGTCATAGTCTTAGCCTTAACTTTAACCGTTTCTGGCCCTTTTAATCCCAcctt                                                                                               >  4:134691/1‑273 (MQ=255)
                                                                                                                                                                             gAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                                   >  8:574761/1‑110 (MQ=255)
                                                                                                                                                                             gAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                                   <  7:574761/110‑1 (MQ=255)
                                                                                                                                                                                     cAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGc                                                                                                                                                                                                                                                                <  5:255761/105‑1 (MQ=255)
                                                                                                                                                                                     cAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGc                                                                                                                                                                                                                                                                >  6:255761/1‑105 (MQ=255)
                                                                                                                                                                                      aGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCg                                                                                                                                                                  >  4:3634/1‑198 (MQ=255)
                                                                                                                                                                                          cATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGt                                                                                                                                                              <  2:3634/198‑1 (MQ=255)
                                                                                                                                                                                                 tCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCg                                                                                                                                                                                                                           <  8:1131290/130‑1 (MQ=255)
                                                                                                                                                                                                 tCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCg                                                                                                                                                                                                                           >  7:1131290/1‑130 (MQ=255)
                                                                                                                                                                                                  cGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCNTAGACCGGAAAAGGCGTTCACGCCGAATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAANCGTAATAGTTATAGCGGTAAATGTAATGGGTGCGGGCGCGTTTAANCNCCCCAATGAGGAAAGACACTTTTATAtt                                                                         >  3:47278/1‑275 (MQ=255)
                                                                                                                                                                                                  cGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCCCGGAGTAACTATAACCGTAATCGTTATAGCCGTAACTGTACGCGGTTCTGCCGCGTTTAATCACACCATTGAGTATAACGCCTTTAATAtt                                                                         >  1:42478/1‑275 (MQ=255)
                                                                                                                                                                                                          ggTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcct                                                                  >  4:413976/1‑274 (MQ=255)
                                                                                                                                                                                                              ttACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGcctgcctg                                                             >  2:165097/1‑275 (MQ=255)
                                                                                                                                                                                                                 cGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                                                                                                                                        <  8:903261/69‑1 (MQ=255)
                                                                                                                                                                                                                 cGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTAc                                                                                                                                                                                                                                                                        >  7:903261/1‑69 (MQ=255)
                                                                                                                                                                                                                      aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTc                                                                                                                                                                                                        >  2:506325/1‑128 (MQ=255)
                                                                                                                                                                                                                      aTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                      >  2:364372/1‑274 (MQ=255)
                                                                                                                                                                                                                        cGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACg                                                    <  2:475054/274‑1 (MQ=255)
                                                                                                                                                                                                                          tgaatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCAcgc                                                                                                                                                                                                    <  4:506324/128‑1 (MQ=255)
                                                                                                                                                                                                                            aatgaaGCACGCATACCGGCGTGCAGTTTGTAATTCATAAAAAAGGGCATTAAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCa                                               <  1:47278/275‑1 (MQ=255)
                                                                                                                                                                                                                            aatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAAc                                                                                                                                                                                                                                      >  8:495264/1‑92 (MQ=255)
                                                                                                                                                                                                                            aatgaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAAc                                                                                                                                                                                                                                      <  7:495264/92‑1 (MQ=255)
                                                                                                                                                                                                                                   cACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTg                                                                                                                                                                                                                                    <  6:665142/87‑1 (MQ=255)
                                                                                                                                                                                                                                   cACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTg                                                                                                                                                                                                                                    >  5:665142/1‑87 (MQ=255)
                                                                                                                                                                                                                                              gCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCtt                             <  2:413977/275‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACaa                                                                                                                                                                                                                                               <  7:170367/56‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACaa                                                                                                                                                                                                                                               >  8:170367/1‑56 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTa                                                                              >  4:224283/1‑206 (MQ=255)
                                                                                                                                                                                                                                                                      aaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAatat                                                                          <  2:224284/206‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gagCATTCAGTTACCGGAATGGCCTGAGGCTGATGACAAACGAAGAACTGGCTAATCCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCa   <  4:364371/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                         agCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCa   >  2:77049/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          gCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCaa  >  2:392923/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          gCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGcc                                                                                                                                                                                         <  5:256001/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGcc                                                                                                                                                                                         >  6:256001/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                             ttCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                                                                                                                                                   <  7:333936/46‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ttCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                                                                                                                                                   >  8:333936/1‑46 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                         
GCCGCCAGGTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAA  >  NC_000913/1041663‑1042202

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: