Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,940 C→A 8.8% intergenic (+25/+90) 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,9400CA 8.8% 228.3 / 11.8 68intergenic (+25/+90)appA/etkphosphoanhydride phosphorylase/tyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base A (3/3);  total (30/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 = 7.29e-02

GTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCT  >  NC_000913/1041671‑1042212
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
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:450782/275‑1 (MQ=255)
                            tCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaagaa                                                                                                                                                                                                                                                 <  2:193911/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                                                                                                                                                                              >  2:511661/1‑227 (MQ=255)
                                                                                                                                               aCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtct                                                                                                                                                                              <  4:511660/227‑1 (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                                                                                                                                                                                                                                                                             <  7:196120/114‑1 (MQ=255)
                                                                                                                                                                 gagCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGaa                                                                                                                                                                                                                                                                             >  8:196120/1‑114 (MQ=255)
                                                                                                                                                                    cGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACAGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTACCTATAACCGTCATAGTCTTAGCCTTAACTTTAACCGTTTCTGGCCCTTTTAATCCCAcctt                                                                                                         >  4:134691/1‑273 (MQ=255)
                                                                                                                                                                     gAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                                             <  7:574761/110‑1 (MQ=255)
                                                                                                                                                                     gAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGaa                                                                                                                                                                                                                                                                             >  8:574761/1‑110 (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                                                                                                                                                                                                                                     >  7:1131290/1‑130 (MQ=255)
                                                                                                                                                                                         tCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCg                                                                                                                                                                                                                                     <  8:1131290/130‑1 (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                                                                                                                                                                                                                                              >  5:665142/1‑87 (MQ=255)
                                                                                                                                                                                                                           cACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTg                                                                                                                                                                                                                                              <  6:665142/87‑1 (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                                                                                                                                                                                                   >  6:256001/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                  gCATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGcc                                                                                                                                                                                                   <  5:256001/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ttCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                                                                                                                                                             <  7:333936/46‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ttCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                                                                                                                                                             >  8:333936/1‑46 (MQ=255)
                                                                                                                                                                                                                                                                           tACCTGAATGCGATGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCt  <  2:134692/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
GTGGTGAACTGGTGTTTGAACGCTGGCGTCGGCTAAGCGATAACAGCCAGTGGATTCAGGTTTCGCTGGTCTTCCAGACTTTACAGCAGATGCGTGATAAAACGCCGCTGTCATTAAATACGCCGCCCGGAGAGGTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCT  >  NC_000913/1041671‑1042212

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: