Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 377,492 G→A intergenic (+181/+43) yaiL → / ← frmB DUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913377,4920GA83.3% 43.4 / 2.7 25intergenic (+181/+43)yaiL/frmBDUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (14/6);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAACCCCTGATCA  >  NC_000913/377254‑377724
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
gCAGGATGAAGACGACCCGTAGGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTTCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCg                                                                                                                                                                                                                              <  1:137955/251‑1 (MQ=255)
                                                                                                                                    gggtacaCCGGATAAGGCGTCCCCGCCCGCTCCGGCAGGCGGGCCCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa                                                                                          <  2:147800/245‑1 (MQ=255)
                                                                                                                                                  atttgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTATTGATAtt                                                                              >  1:534026/24‑251 (MQ=18)
                                                                                                                                                     tgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCAt                                                                           >  1:513249/21‑251 (MQ=18)
                                                                                                                                                           cccccccacccGGCCGTCGTGCGCCATGCCTGATGCGACGTTTTTCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGTGTTCACGCCGCATCCGGCATTCGTGCTATTATCAACGCATATTCAGTTTATTGTCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTTATAACGGATTAACCTCTTGATATTCATCTCCTgg                                                                    <  2:246382/242‑1 (MQ=255)
                                                                                                                                                             taaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                                                                   >  1:144981/13‑251 (MQ=21)
                                                                                                                                                               aacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGACGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGACAATATTCTCGCCAATAAAACTAGAGACAAAATAATAGCTGAGATCATAACCCTCGTAATAACGGCTTAACGACTTGGTATTCCTCTCCttggcg                                                                 >  1:81996/11‑246 (MQ=17)
                                                                                                                                                                 catagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                                                                   >  1:36749/9‑247 (MQ=25)
                                                                                                                                                                       taaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCCCGGGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc                                                         >  1:352085/3‑251 (MQ=25)
                                                                                                                                                                       aaaggggtgcgcagaTGCCTTATGCGACGTTATGCGCGTTTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGGGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc                                                         <  2:473379/237‑1 (MQ=25)
                                                                                                                                                                         aGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGGAGTGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGAAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCCTAACCCTCGTGATAAAGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTcc                                                        >  2:187822/1‑249 (MQ=25)
                                                                                                                                                                          gTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTTTGATCCTAACCCTCGTGATAACGGATTAACGTCTTGATATTAATCTCCTGGCAGATCTTTTCa                                                        >  2:359808/1‑249 (MQ=25)
                                                                                                                                                                            cGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGAtt                                                   <  2:23646/251‑1 (MQ=255)
                                                                                                                                                                               gcgcagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTg                                                 <  1:323473/245‑1 (MQ=255)
                                                                                                                                                                               gcgcagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGACAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCACTAAAACTGGAAACAATATAATAACTGTGATCAATACCCTCGTGATAAAGGTTTAACGTCTTAATATTCATCTCCTGGCATATCCTTTCAAGAtttt                                                 >  2:212093/7‑250 (MQ=21)
                                                                                                                                                                                  cagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGACTTGATATTAATCTCCTGGCAGATCTTTTCAAGATTTgg                                                >  1:32312/4‑249 (MQ=37)
                                                                                                                                                                                                gACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGGGGTAGGCAAAAAGCTCGCCACTAAAACAGGAGACAAAATAATAGCAGTGATCATAACCCTCGTGATAACGGATTAACGGCTTGATATTCATCTCCTGGCAGGTCTTTTAAAAATTTTGAGTCCGCAGCGGtttt                                >  2:393610/1‑249 (MQ=25)
                                                                                                                                                                                                 ccgctattggcgttttttCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTTGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATACCCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTg                               <  2:470939/234‑1 (MQ=38)
                                                                                                                                                                                                            gTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAatca                    <  1:250256/251‑1 (MQ=255)
                                                                                                                                                                                                               ttagccggCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAatcatc                  <  2:257692/244‑1 (MQ=255)
                                                                                                                                                                                                                 aTCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACt               >  1:117054/1‑251 (MQ=255)
                                                                                                                                                                                                                         tacgtgttatgggCTAAGTGTAGCCGGGATAAGCCTTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCAGATCCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAAccc         <  1:309378/239‑1 (MQ=255)
                                                                                                                                                                                                                          acgttttatgggCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAACCCCTg      <  1:273327/242‑1 (MQ=255)
                                                                                                                                                                                                                            ggtttatgggCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAAGTCTTGCTATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTCAAAATCATCACTAAACCCCTGAt    >  2:145538/8‑251 (MQ=39)
                                                                                                                                                                                                                               ttatgggCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAAAGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGATGTGATCATAACCCTCGGGAAAACGGATTAACGACTTGATATTCATCTCCTGGCAGATCTTTACAAGATTTTGAATCTGAAGCAGTTTTGCGTAAAAATCCTCACTCACCCCCTGATCa  >  2:302204/5‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
GCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACTCAACCCCTGATCA  >  NC_000913/377254‑377724

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