Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 340,089 G→A intergenic (+346/‑76) yahG → / → yahI DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913340,0890GA93.8% 89.5 / ‑3.3 32intergenic (+346/‑76)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (16/14);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCC  >  NC_000913/339974‑340221
                                                                                                                    |                                                                                                                                    
tACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAg                                                                                                               <  2:971547/139‑1 (MQ=255)
  cGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAgg                                                                                                              >  2:1301474/1‑139 (MQ=255)
     gTGGCTCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAgg                                                                                                              <  2:996762/135‑1 (MQ=255)
     gTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGAGAAGCGTAGCGCATCAgg                                                                                                              >  1:996762/1‑135 (MQ=255)
            gagaATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCg                                                                                                    <  1:1301474/139‑1 (MQ=255)
                                       cgcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                                                        >  2:1078639/1‑139 (MQ=255)
                                            cAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGAc                                                                   >  1:803468/1‑139 (MQ=255)
                                                   ttGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCa                                                                                                                   >  2:684811/1‑84 (MQ=255)
                                                   ttGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCa                                                                                                                   <  1:684811/84‑1 (MQ=255)
                                                   ttGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaa                                                               >  2:276728/1‑136 (MQ=255)
                                                   ttGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaa                                                               <  1:276728/136‑1 (MQ=255)
                                                             tGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAcc                                                                                >  2:97792/1‑109 (MQ=255)
                                                             tGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAcc                                                                                <  1:97792/109‑1 (MQ=255)
                                                             tGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAagag                                                  >  1:1198802/1‑139 (MQ=255)
                                                                       ggcgTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATc                                                                                   <  1:107326/96‑1 (MQ=255)
                                                                       ggcgTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATc                                                                                   >  2:107326/1‑96 (MQ=255)
                                                                       ggcgTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCg                                        <  1:1078639/139‑1 (MQ=255)
                                                                         cgTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGtt                                      >  1:651816/1‑139 (MQ=255)
                                                                             aaCGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCg                                        >  1:34213/1‑133 (MQ=255)
                                                                             aaCGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCg                                        <  2:34213/133‑1 (MQ=255)
                                                                             aaCGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCa                                  <  2:803468/139‑1 (MQ=255)
                                                                                   ttAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAgg                                                                                                              <  2:1059861/57‑1 (MQ=255)
                                                                                   ttAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAgg                                                                                                              >  1:1059861/1‑57 (MQ=255)
                                                                                                cGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCa                                                                                     <  2:1250448/70‑1 (MQ=255)
                                                                                                cGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCa                                                                                     >  1:1250448/1‑70 (MQ=255)
                                                                                                cGAGTGGCGCGAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTa                >  1:612214/1‑139 (MQ=255)
                                                                                                      gcgcGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCa                   >  2:1226068/1‑130 (MQ=255)
                                                                                                      gcgcGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCa                   <  1:1226068/130‑1 (MQ=255)
                                                                                                         cgAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaaca                                                             >  2:166681/1‑85 (MQ=255)
                                                                                                         cgAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaaca                                                             <  1:166681/85‑1 (MQ=255)
                                                                                                         cgAGAATTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATa       >  2:708165/1‑139 (MQ=255)
                                                                                                              aTTGTAAGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTCTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgcc  <  2:1198802/139‑1 (MQ=255)
                                                                                                                    |                                                                                                                                    
TACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAT‑CCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCC  >  NC_000913/339974‑340221

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