Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 337,142 A→G N122S (AAC→AGC)  yahF → putative acyl‑CoA synthetase YahF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913337,1420AG93.3% 83.7 / ‑2.0 30N122S (AAC→AGC) yahFputative acyl‑CoA synthetase YahF
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (15/13);  total (16/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 = 8.88e-01

GTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGG  >  NC_000913/337021‑337266
                                                                                                                         |                                                                                                                            
gtggAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTTATTACGCTCAGCGCTCTGTTTg                                                                                                                   >  1:214810/1‑133 (MQ=255)
gtggAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTg                                                                                                                   >  1:100113/1‑133 (MQ=255)
    agaTTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGc                                                                                                               >  1:66858/1‑133 (MQ=255)
                tGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGcgcg                                                                                                   >  1:115708/1‑133 (MQ=255)
                  ttCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTATTTGCCGCTCGCGAAGCGCGTc                                                                                                 >  1:156653/1‑133 (MQ=255)
                  ttCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGGGATTTCGGGCAGGTATACGTTTGCCGCTCGCAAAGCGCGTc                                                                                                 >  2:238952/1‑133 (MQ=255)
                  ttCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGCGATTTCGGTCAGCGGTTTGTTTGCCCCTCGCGAAGCGCGTc                                                                                                 >  2:229873/1‑133 (MQ=255)
                            aagcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCa                                                                                       >  2:500466/1‑133 (MQ=255)
                            aagcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTCTGCCGCTTGCGTAGTGCGTCAAGCGCTGCa                                                                                       >  1:234216/1‑133 (MQ=255)
                            aagcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCa                                                                                       >  1:340576/1‑133 (MQ=255)
                            aagcgcaaagcgACTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCa                                                                                       >  2:346555/1‑133 (MQ=255)
                               cgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCaaaa                                                                                    <  2:115708/133‑1 (MQ=255)
                               cgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCaaaa                                                                                    <  1:287199/133‑1 (MQ=255)
                                 caaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACg                                                                                  >  2:284380/1‑133 (MQ=255)
                                    agcgGCTCGCACGAGGCGCGGTACGCCACTAGTGCCAGCGCCAAAAAGCATATCCCGGAAAGTAACCGGGCGGAGATTTCGGTCAGCGGTCTGCTTGCCGCACGCGAAACGCGTCCGGCGCTTCAAAACCAtc                                                                               >  2:35902/1‑133 (MQ=255)
                                       gacTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCa                                                                            <  1:346555/131‑1 (MQ=255)
                                               cGAGGCGCGTAACGCGACTATTGGCAGCCCCAGAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATg                                                                    <  2:156653/133‑1 (MQ=255)
                                                           cGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGAt                                                        <  2:100113/133‑1 (MQ=255)
                                                              cACTATTGGCCGCGCAAAAAAGCACATCTCGGGAAGTAACCAGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAACCGCGGCAGGCGCTGCAGAACGATCTCAACGTGATGCTGTTTTCCGATAAc                                                     <  2:214810/133‑1 (MQ=255)
                                                                  aTTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                 <  2:340576/133‑1 (MQ=255)
                                                                      gCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGt                                             <  2:66858/133‑1 (MQ=255)
                                                                          cgctaaaaaGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtc                                                  >  1:402404/5‑124 (MQ=255)
                                                                          cgcTAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtc                                                  <  2:402404/124‑1 (MQ=255)
                                                                          cgcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGaa                                         >  2:328270/1‑133 (MQ=255)
                                                                          cgcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGa                                          <  2:408098/132‑1 (MQ=255)
                                                                          cgcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGa                                          >  1:408098/1‑132 (MQ=255)
                                                                           gcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                 <  1:470551/124‑1 (MQ=255)
                                                                           gcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                 >  2:470551/1‑124 (MQ=255)
                                                                           gcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                 >  2:387296/1‑124 (MQ=255)
                                                                           gcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                 <  1:387296/124‑1 (MQ=255)
                                                                                         cccGGAAAGTAACCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGAGAACGTCTGAGTGGACGATGAAATGGCGCTc                          >  1:158024/1‑133 (MQ=255)
                                                                                                    aCCTGGCGGTGATTTCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCCGTTTATGTTGAACTGGGGCTCAAGCAACTgtc               >  2:243483/1‑131 (MQ=255)
                                                                                                     ccTGGCGGTGCTTTCGGTCAGCGGCCTGTTTGTCGCTCGCGGAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGcc              <  2:197297/133‑1 (MQ=255)
                                                                                                                 ttCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAggg  <  1:328270/133‑1 (MQ=255)
                                                                                                                 ttCGGTCAGCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAggg  <  1:480792/133‑1 (MQ=255)
                                                                                                                         |                                                                                                                            
GTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGG  >  NC_000913/337021‑337266

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑