Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,710,250 G→A intergenic (‑86/‑170) rseD ← / → nadB rpoE leader peptide/L‑aspartate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,710,2500GA88.2% 40.9 / ‑0.5 17intergenic (‑86/‑170)rseD/nadBrpoE leader peptide/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (5/10);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

CGCTCATCCGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATC  >  NC_000913/2710006‑2710461
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                   
cgctcATCCGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCttt                                                                                                                                                                                                               <  1:350734/251‑1 (MQ=255)
                 gTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCACAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAAt                                                                                                                                                                                              >  2:60514/1‑251 (MQ=255)
                        caaaccacatcaccACGCGCTAGCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAATGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCAGTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATTCCCATTTAACTGCCTGAAGAGTAACCCAACATACTTTTTATTTCACCACCTAAtcctcccc                                                                                                                                                                                       <  2:130756/239‑3 (MQ=255)
                                            tATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt                                                                                                                                                                   >  1:118690/1‑251 (MQ=255)
                                            tATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt                                                                                                                                                                   >  2:155256/1‑251 (MQ=255)
                                                                aGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc                                                                                                                                               >  2:97219/1‑251 (MQ=255)
                                                                            agcaagcaAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAAc                                                                                                                                   <  1:183705/251‑1 (MQ=255)
                                                                                                       ttttGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCAGGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCGCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGgtg                                                                                                         <  2:43300/250‑1 (MQ=255)
                                                                                                                                               gCTGTAGACGGTTCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTTTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGtt                                                                <  2:118690/251‑1 (MQ=255)
                                                                                                                                               gCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGtt                                                                >  2:55008/1‑251 (MQ=255)
                                                                                                                                                 tGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGntt                                                               >  2:435577/1‑250 (MQ=255)
                                                                                                                                                 tGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAAAATCGCATTATCTGAGTCGAAATGACCATTTAACTGCCTGAAGAGTAACCCAACGGCTTTTTTATTTCACCAGCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGAGTTTATCCGCACAACATGATGCTATACTGACCAAACCATGttt                                                               >  2:482/1‑250 (MQ=255)
                                                                                                                                                     gACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTaa                                                          <  1:425580/251‑1 (MQ=255)
                                                                                                                                                                       aTTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGa                                        <  1:155256/251‑1 (MQ=255)
                                                                                                                                                                          ataAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAAt                                      <  1:444530/250‑1 (MQ=255)
                                                                                                                                                                            aaGTCTACAGAATCTGAACATCGCATTATCTGTGGAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATAct                                   <  1:435577/251‑1 (MQ=255)
                                                                                                                                                                                                              tAGAAATGCCCATTTAACTGCCGTAAGAGTACCCCAACAGCCTTTTTATTTCCCCACCTAACCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTTTTGATTATc  <  2:308120/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                   
CGCTCATCCGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATC  >  NC_000913/2710006‑2710461

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