Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,126,316 G→A E54K (GAA→AAA)  yceH → DUF480 domain‑containing protein YceH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,126,3160GA86.2% 56.0 / 1.4 29E54K (GAA→AAA) yceHDUF480 domain‑containing protein YceH
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (13/12);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAA  >  NC_000913/1126088‑1126552
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
aaTGGCGCGATCACGTACTGACTGCAATAACTACTCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCGTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAAGACCCGCACTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGATCCGTGAACTGGTCATGAATCTGAGTGAATTCAAAGTGCAGGATCAGCTGGATaa                                                                                                                                                                                                                        <  2:268490/251‑1 (MQ=255)
              gTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGGCATGAATCTGAGTGAATCCAACGTGCAGGAACAGCTGGAGAATCTGGTCAAACGTc                                                                                                                                                                                                          >  2:188782/1‑251 (MQ=255)
                            aaCTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGACCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACAACGCAGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAGACGAACCGTGAAACGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGACAAACGTCATTATCTACGcac                                                                                                                                                                                             >  1:33123/1‑250 (MQ=255)
                              cTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGCCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGt                                                                                                                                                                                          >  1:354010/1‑251 (MQ=255)
                               tACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTACAGGAACAGCTGAATAATCTGGTCAAAAGTCATTATCTACGCACAGtt                                                                                                                                                                                         >  1:5476/1‑250 (MQ=255)
                                     aGACTGGACGCCCGTCCGCTAAGGAGAGTCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGATCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACCGTGAGCGGt                                                                                                                                                                                   <  2:260475/251‑1 (MQ=255)
                                          ggACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTg                                                                                                                                                                               <  1:227749/250‑1 (MQ=255)
                                            aCGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAACCGACATTATCTACGCACAGTGAGCGGTTTTGGt                                                                                                                                                                             >  2:53088/1‑250 (MQ=255)
                                                ccgcccGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCg                                                                                                                                                                        <  1:220913/247‑1 (MQ=255)
                                                                          atatCAACTTACGGCCCTTGAAGCTCGCTTGATGGGTTGCTTGCTGGAAAAGCATGTCATCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGTCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTTAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGttttt                                                                                                                                              <  2:133073/251‑1 (MQ=255)
                                                                                aCTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAAtt                                                                                                                                        >  1:36450/1‑251 (MQ=255)
                                                                                     cggCCCCTGAAGCTCGCGTTATTTGTTGCTTGCTGGAAAAGCAGTTGACCACTACGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGTATAATCTGGTCAAACGTCATGATCTACGCACTGTGAGCGGTTTTGGTAAGCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGaa                                                                                                                                   <  2:338073/249‑1 (MQ=255)
                                                                                         acTTGCAGCTCGCGTGATTGGATGCTTGCTGGAAAAGCAGGTGACCACACCGGCGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAACCGTCATTATCTACGCACAGTGAGCGGTTTTGATAATCGGGTCACCAAATATGAGCAACATTTTTGAAATTCAGAATTTg                                                                                                                               >  2:15074/2‑251 (MQ=255)
                                                                                            tGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTgg                                                                                                                              <  1:114980/249‑1 (MQ=255)
                                                                                                                    gctGGAAAAGCAGGTGACCTCGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACTCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGaa                                                                                                     <  1:343010/250‑1 (MQ=255)
                                                                                                                      tGGAAAAGCAGGAGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGGGTAGTCACGGCCTGTCATCAGAAACCGACCCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAAAAGCTGGATAATCTGGTAAAAAGACATTATCTACGAAAAGTGAGCGGTTTTGGTTATCGGGTCACCAAATATGAGCAAAGTTTTTGTAATTCAGAATTTGGCGATCTAAAAcagcgagaagcggacgtg                                                                                                  >  1:358552/1‑234 (MQ=255)
                                                                                                                               aGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGACAACGAATCGTGAACCGGTCATGAATCTGAGTGTATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAAGATGAGCAACGTGTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATc                                                                                         <  2:141570/251‑1 (MQ=255)
                                                                                                                                ggTGACCACGCCGGAGCAATACCCGCTCGCAGGCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGGGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGCACAGCTTGATAATCTTGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTGGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTTATTCAGAATGTGGCGATCTGAAACTGTGCGCAGCGGAAGTGGCGTTAATca                                                                                        <  2:147463/251‑1 (MQ=255)
                                                                                                                                     ccACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACg                                                                                   <  1:53088/251‑1 (MQ=255)
                                                                                                                                                             caggcAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAAACGGTCATGAAGCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTTGTCAATCGACATTATCTACGCACAGTGAGCGGTTGTGGTAATCGGGTCAACAAATATGAGCAACGTTTGTGTAAGTCAGAATTTGGCGTTCTGAAATTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTTCGTGGTGCCCAGAc                                                            <  2:57174/246‑1 (MQ=255)
                                                                                                                                                                                                  gAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATg                      >  1:267155/1‑251 (MQ=255)
                                                                                                                                                                                                   aaCCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCCAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCCCCAAATATGGGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGATGCCCAGACGCCAGGTGAGCTGCGCAGCCGCGTCGCGCGAATGTATg                      >  1:81077/1‑250 (MQ=255)
                                                                                                                                                                                                    aCCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGa                     >  1:106990/1‑250 (MQ=255)
                                                                                                                                                                                                     ccGTGAACCGGTCATGAATCTGAGAGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAAt                   >  1:148025/1‑251 (MQ=255)
                                                                                                                                                                                                        tGAACCGGTCATGAATATGAGTGAATCCGAAGTGCAGGAACAGCTGGCTAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGCTCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGGGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTc                 >  2:21165/1‑250 (MQ=255)
                                                                                                                                                                                                            ccGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTATAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCgt            >  1:348378/1‑250 (MQ=255)
                                                                                                                                                                                                                          gAGTGAATCCAAAGTGCATGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCATCCGCGTCGCGCTAATGTATGAATTCAGCGATATGGCGGaa  <  1:322253/247‑1 (MQ=255)
                                                                                                                                                                                                                                 tCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAgg                                                      >  1:236986/1‑188 (MQ=255)
                                                                                                                                                                                                                                 tCCGAAGGTCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTTGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAgg                                                      <  2:236986/188‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
AATGGCGCGATCACGTACTGACGGCATTAACTACCCCAGACTGGACGCCCGGCCGCTAAGGAGAATCTCATGAAATATCAACTTACTGCCCTTGAAGCTCGCGTGATTGGTTGCTTGCTGGAAAAGCAGGTGACCACGCCGGAGCAATACCCGCTCTCAGTCAATGGTGTAGTCACGGCCTGTAATCAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAA  >  NC_000913/1126088‑1126552

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