Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 354,035 G→T G444W (GGG→TGG)  prpE → propionyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,0350GT94.1% 51.1 / ‑2.6 17G444W (GGG→TGG) prpEpropionyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (6/10);  total (6/11)
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.67e-01

CTGGACGACAGGCCGACGCGTCTGGGAAGCCCCGGTGTGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAAGGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCGGTGTACGCCACCTTTGACTGGGGCATCCGTGACGCTGAC  >  NC_000913/353909‑354160
                                                                                                                              |                                                                                                                             
cTGGACGACAGGCCGACGCGTCTGGGAAGCCCCGGTGTGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTgccgcc                                                                                                                  <  1:562456/140‑1 (MQ=255)
   gacgacAGGCCGACGCGTCTGGGAAGCCCCGGTGTGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCggg                                                                                                               <  1:383210/140‑1 (MQ=255)
                             ccccGGTGTGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGcgacg                                                                                     >  1:454611/1‑140 (MQ=255)
                                cGGTGTGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGacgacg                                                                                  >  2:469082/1‑140 (MQ=255)
                                  gtgtGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGc                                                                                <  1:638749/140‑1 (MQ=255)
                                             tATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAAGGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGaa                                                                     <  2:268723/140‑1 (MQ=255)
                                                      aaCGTGCAGCTGCTTAATGAAGTCACCGGCGAACCGGGTGGCAACAACGCGAAATGGATGCTGGCGGTGGAATGGCCGCTGCCGCCGGGGTGCATTCAGACCATCTGGGGCGACGACGGCCGCTTTGGGAAGACTTACTg                                                            >  2:382224/1‑140 (MQ=255)
                                                        cGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGt                                                          <  2:267228/140‑1 (MQ=255)
                                                            cAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCt                                                      <  2:294506/140‑1 (MQ=255)
                                                               ttGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGtt                                                   >  2:387226/1‑140 (MQ=255)
                                                                     aaTGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTcccg                                             >  2:69015/1‑140 (MQ=255)
                                                                     aaTGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTcccg                                             >  2:421655/1‑140 (MQ=255)
                                                                       tGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTcccgcc                                           <  2:454611/140‑1 (MQ=255)
                                                                          aGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCgg                                        <  1:387226/140‑1 (MQ=255)
                                                                                                   aaCGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCGGTGTACGCCACCTTTGACTGGGGCAt               <  1:69015/140‑1 (MQ=255)
                                                                                                   aaCGAGAAAGGGATGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCGGTGTACGCCACCTTTGACTGGGGCAt               <  1:421655/140‑1 (MQ=255)
                                                                                                               aTGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCGGTGTACGCCACCTTTGACTGGGGCATCCGTGACGCTGa   <  1:469082/140‑1 (MQ=255)
                                                                                                                tGCTGGTGGTGGAATGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCGGTGCACGCCACCTTTGACTGGGGCATCCGTGACGCTGAc  >  2:12825/1‑140 (MQ=255)
                                                                                                                              |                                                                                                                             
CTGGACGACAGGCCGACGCGTCTGGGAAGCCCCGGTGTGCCGATGTATGGCTATAACGTGCAGTTGCTTAATGAAGTCACCGGCGAACCGTGTGGCGTCAACGAGAAAGGGATGCTGGTGGTGGAAGGGCCGCTGCCGCCGGGGTGTATTCAGACCATCTGGGGCGACGACGGCCGCTTTGTGAAGACTTACTGGTCGCTGTTTTCCCGCCCGGTGTACGCCACCTTTGACTGGGGCATCCGTGACGCTGAC  >  NC_000913/353909‑354160

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