Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,018,977 C→A intergenic (+40/+36) matP → / ← ompA macrodomain Ter protein/outer membrane protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9770CA90.0% 51.5 / ‑2.2 20intergenic (+40/+36)matP/ompAmacrodomain Ter protein/outer membrane protein A
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (10/8);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGC  >  NC_000913/1018733‑1019222
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
tGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCgc                                                                                                                                                                                                                                                  >  1:230313/1‑250 (MQ=255)
                   cgcAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                                                                                                                              <  2:356560/251‑1 (MQ=255)
                                                   cTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAACAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAc                                                                                                                                                                                               >  1:361045/1‑250 (MQ=255)
                                                   cTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACgagatcggaagacacacgtctgaactcc                                                                                                                                                                                              >  1:187783/1‑227 (MQ=255)
                                                   cTGGCAACGGCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACg                                                                                                                                                                                                                         <  2:187783/224‑1 (MQ=255)
                                                                   gggcTTGCGCAGCGTCGCGGTAAAACGCGTGCTGAAACGATTGTACATCTGATGTAAGATGCGGAAAACAAAGAGAAATACGCGAAGAAAATGTCTTCACTGAAGCAAGATCTGCAGGCATTGCTGGGTAAGGAATAACGGACGAAAGTCAGTTCAATTTCCTAAAGGCAAAAAAAAACCCGGAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCtt                                                                                                                                                                              <  2:172488/248‑1 (MQ=255)
                                                                                          aaCGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacg                                                                                                                                                        >  2:335216/1‑250 (MQ=255)
                                                                                            cGCTTTCTGAAACGAATTTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACGAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTGTTTCTACCAGACGAGAACGTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATc                                                                                                                                                     <  2:140700/251‑1 (MQ=255)
                                                                                                ttCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgg                                                                                                                                                  <  2:228475/250‑1 (MQ=255)
                                                                                                  cTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                                                                                                                               >  2:93504/1‑251 (MQ=255)
                                                                                                      aaCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGCGAAATACGCGAATAAACTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAACCAgg                                                                                                                                           >  2:298956/1‑251 (MQ=255)
                                                                                                                  gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt                                                                                                                               >  2:169975/1‑251 (MQ=255)
                                                                                                                             aTGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAGGGAATCACTGACGAAAGTCAGTTCAATTTACTCAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGCCGCGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGCGCAGCACGCTg                                                                                                                    >  1:414088/1‑251 (MQ=255)
                                                                                                                                cGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                                                                                                                 <  1:173218/251‑1 (MQ=255)
                                                                                                                                       cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTg                                                                                                          >  1:126638/1‑251 (MQ=255)
                                                                                                                                                                                  cTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCAc                                                               >  1:369246/1‑251 (MQ=255)
                                                                                                                                                                                                         aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGAAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTAGCCCATACCACGTGCGGAGATCTTGTCTGCCggg                                        >  1:388687/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCttt                                >  1:49711/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaaGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCATCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGACAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTgg                              <  1:88582/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCGGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATca                <  1:217430/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaaaaCCCGCAGCAGCGGGGTTTTTCTACCAGAAGAGAACTTAAGCATGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCATTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGTAGATCAGGTAATCAACCACAGACTGAGc  >  1:116997/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
TGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGC  >  NC_000913/1018733‑1019222

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