Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,595,493 C→T G120D (GGT→GAT)  lgoR ← putative DNA‑binding transcriptional regulator LgoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,595,4930CT95.2% 62.5 / ‑4.8 21G120D (GGT→GAT) lgoRputative DNA‑binding transcriptional regulator LgoR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (11/9);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCG  >  NC_000913/4595375‑4595627
                                                                                                                      |                                                                                                                                      
ctgctCGACGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCa                                                                                                                   >  2:59630/1‑140 (MQ=255)
      ggCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAg                                                                                                             >  2:147224/1‑140 (MQ=255)
        cGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGaa                                                                                                           <  2:280093/140‑1 (MQ=255)
               ttGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGttt                                                                                                    >  1:348005/1‑140 (MQ=255)
                tGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTc                                                                                                   >  1:131345/1‑140 (MQ=255)
                        ttGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGc                                                                                           <  1:223598/140‑1 (MQ=255)
                                 gCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGgc                                                                                  <  2:224105/140‑1 (MQ=255)
                                  cTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGgcg                                                                                 <  2:348005/140‑1 (MQ=255)
                                  cTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTggtg                                                                                 <  1:409615/140‑3 (MQ=255)
                                           ccGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTat                                                                        >  1:135084/1‑140 (MQ=255)
                                            cgcgCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTatc                                                                       >  2:242962/1‑140 (MQ=255)
                                                           gAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGcc                                                        <  1:242962/140‑1 (MQ=255)
                                                                      taACGGCCGCATTTTAAAAGATATTCTCGTACAACAACCGGGCTGACATCCGCGGCCGGAGCCAGGTGCAGTTCAGAAAACGTTTCCCCTGGGCGAAACTGGCGCTGGCTTATCATCGAGAAATCAGCCTGCTCAAacac                                             >  2:522130/1‑140 (MQ=255)
                                                                            ccGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACtttgtt                                       >  2:265554/1‑140 (MQ=255)
                                                                            ccGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACttt                                          <  2:101448/137‑1 (MQ=255)
                                                                            ccGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACttt                                          >  1:101448/1‑137 (MQ=255)
                                                                                           tatTCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGTGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGa                        <  1:59630/140‑1 (MQ=255)
                                                                                             ttCTCTTACTACAACCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAg                      <  2:135084/140‑1 (MQ=255)
                                                                                                          aCCGGGCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGc         >  1:388202/1‑140 (MQ=255)
                                                                                                              ggCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGaaa     >  2:135762/1‑140 (MQ=255)
                                                                                                              ggCTGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGaaa     >  2:288473/1‑140 (MQ=255)
                                                                                                                 tGACATCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCg  <  1:288473/140‑1 (MQ=255)
                                                                                                                      |                                                                                                                                      
CTGCTCGGCGTAGGATTGATCGAATTGCTTCATGCTCCACTGGCCGCGCTTTTCGCTGTGAATAAGATTATAACGGCCGAATTTTAAAAGATATTCTCTTACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAACCG  >  NC_000913/4595375‑4595627

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