Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,909,387 A→T L238* (TTA→TAA)  kdgR ← DNA‑binding transcriptional repressor KdgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,909,3870AT92.3% 83.1 / ‑0.2 26L238* (TTA→TAA) kdgRDNA‑binding transcriptional repressor KdgR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (12/12);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

AAATAACCAGATTCTCACTCCTGAGTAAAATACTCTTATAGCAAGCGCGTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTAAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTT  >  NC_000913/1909172‑1909624
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
aaataaCCAGATTCTCACTCCTGAGTAAAATACTCTTATAGCAAGCGCGTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGga                                                                                                                                                                                                             <  1:55781/250‑1 (MQ=255)
                   ccTGAGTAAAATACTCGTATAGCAAGCGCGTTTAATATAACTCCACCTTCAACTTTAGTTTTCTCATGCTAAATAAAGAAGCTGTCTGTTTGATATTCCGACAGTTTCTTACTCTCAGCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTTTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGa                                                                                                                                                                                         <  1:241180/251‑1 (MQ=255)
                      gAGTAAAATACTCTTATAGCAAGCGCGTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGAc                                                                                                                                                                                        <  1:335442/249‑1 (MQ=255)
                                          aaGCGCGTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTAAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGTTGCTCAAACCGGCAATgcac                                                                                                                                                                                       >  2:89722/1‑227 (MQ=255)
                                          aaGCGCGTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTAAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGAcc                                                                                                                                                                                       <  1:89722/230‑1 (MQ=255)
                                           aGCGCGTTTAATATAACGACACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCg                                                                                                                                                                 <  1:175238/251‑1 (MQ=255)
                                             cgcgTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGt                                                                                                                                                               >  1:249338/1‑251 (MQ=255)
                                                    aaTATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCa                                                                                                                                                         <  1:230570/250‑1 (MQ=255)
                                                               aCCTTCAACTTTATTTTTCTCATGCTAAAGAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGc                                                                                                                                             <  1:5748/251‑1 (MQ=255)
                                                                ccTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGcc                                                                                                                                            >  1:102729/1‑251 (MQ=255)
                                                                                         aaataaaGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCt                                                                                                                     <  1:175526/249‑1 (MQ=255)
                                                                                            taaaGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGAGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTAGGGAAGGAGATGCTCAAACCGGCAATGAACACGCCAAAGCGATCCAATCCCGGTACAGCAATGCATCGCAGCCCTTCTTCCTTCTCATCATTATCTTCGcc                                                                                                                >  2:165314/1‑251 (MQ=255)
                                                                                                  aactgtctgaTTGGTATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAAAGTCGGGAAGAAGATGCTCAAACCGGCAATGACCACGACAAAGCGATCACATACCGGTACCACAATGCAACGCAGCCCTTCTACCTGCTCTTCATGATCTTCGACATAccc                                                                                                          >  2:77968/3‑251 (MQ=255)
                                                                                                      gtctgtTTGATATTCCGACAGGTTCTTACTCTCATCAGCAAGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGAGCGATGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAAc                                                                                                                                                                                                                             >  2:230402/1‑132 (MQ=255)
                                                                                                      gtctggTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAAc                                                                                                                                                                                                                             <  1:230402/132‑1 (MQ=255)
                                                                                                            ttGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATGCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTcgcg                                                                                                 <  1:57319/250‑1 (MQ=255)
                                                                                                                    ccGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGt                                                                                         <  1:235415/250‑1 (MQ=255)
                                                                                                                              cTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCCCCTGGTCCAGAACGGGt                                                                              >  2:284575/1‑251 (MQ=255)
                                                                                                                                     ccatcaGAACGGATAGGCGTGATAACCCATTTGGGCAGAAATTTTGCGCGATGCGGTATGCAACATTGCGACATATGCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCTCCATACCCCTGCTCGCGCACCTGGTACAGAACGGGTAATAACg                                                                       <  2:51010/250‑1 (MQ=255)
                                                                                                                                                 aTAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTgg                                                           >  1:202729/1‑251 (MQ=255)
                                                                                                                                                              cccATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATAGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTc                                              >  1:319028/1‑251 (MQ=255)
                                                                                                                                                                cATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGCGATGGTCCGCTCgg                                            >  1:222624/1‑251 (MQ=255)
                                                                                                                                                                                  gcgcgcTGCGGTATGCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTcc                           <  1:284575/250‑1 (MQ=255)
                                                                                                                                                                                                gCAACATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAAtt            >  2:260593/1‑251 (MQ=255)
                                                                                                                                                                                                  aaCATTGCGACATATTCCTGTTAACGCTCTTCAGAGAAACGCAAAGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGAGTTTATACTCCACGCCCTCAAGAATTTg          >  1:300976/1‑251 (MQ=255)
                                                                                                                                                                                                          gACATATTCCTGTTAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCCTACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGATCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGcttcactt  >  1:206335/1‑251 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
AAATAACCAGATTCTCACTCCTGAGTAAAATACTCTTATAGCAAGCGCGTTTAATATAACGCCACCTTCAACTTTATTTTTCTCATGCTAAATAAAGAACCTGTCTGTTTGATATTCCGACAGGTTCTTACTCTCATCAGAACGGATAGTCGTGATAACCCATTTGGGCAGAAATTTTGCGCGCTGCGGTATGCAACATTGCGACATATTCCTGTAAACGCTCTTCAGAGAAACGCAACGTCGGGAAGGAGATGCTCAAACCGGCAATGACCACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTT  >  NC_000913/1909172‑1909624

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