Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,570,512 C→T 100% intergenic (‑23/+133) gadC ← / ← gadB L‑glutamate:4‑aminobutyrate antiporter/glutamate decarboxylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,5120CT80.2% 38.3 / 1.9 20intergenic (‑23/+133)gadC/gadBL‑glutamate:4‑aminobutyrate antiporter/glutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (13/3);  total (15/5)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.

TTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTA  >  NC_000913/1570381‑1570644
                                                                                                                                   |                                                                                                                                    
ttGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAAcc                                                                                                                           >  2:289476/1‑143 (MQ=255)
                    tAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAg                                                                                                       >  1:369154/1‑143 (MQ=255)
                     aaaCAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGc                                                                                                      >  1:276565/1‑143 (MQ=255)
                      aaCAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCg                                                                                                     >  2:158955/1‑143 (MQ=255)
                       aCAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGt                                                                                                    >  1:251752/1‑143 (MQ=255)
                       aCAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGGATTTCTGTCGGAACCGCACCTGTTTTGTTCTGGAGCGt                                                                                                    >  1:163841/1‑143 (MQ=255)
                               tACCCTCGATGCTGTTATTGCAAAGTATCCACGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGaa                                                                                            <  2:154712/143‑1 (MQ=255)
                                                                            ctTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGACATTGAACAGACATTGGAAATTg                                               >  1:286358/1‑143 (MQ=255)
                                                                                  tttACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCAGa                                         >  1:391774/1‑143 (MQ=255)
                                                                                  tttACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTTGAAATTGTCCCGa                                         >  2:244124/1‑143 (MQ=255)
                                                                                               aCTGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGttt                            >  1:319584/1‑143 (MQ=255)
                                                                                                 tggtgcaGCCATATTATCCCCCTACCACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATGCAGAGCAATATTGTCCGTTGCTAATGCCAGTGAACGGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCg                          <  1:144261/137‑1 (MQ=255)
                                                                                                 tGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGTAAATTGTCCCGCAACGGGTCCGTTTCg                          >  1:178553/1‑143 (MQ=255)
                                                                                                 tGATGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCg                          <  1:158955/143‑1 (MQ=255)
                                                                                                   aTGTAGCCATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCGGAAACGGGTTCGGTTCggg                        >  2:28860/1‑143 (MQ=255)
                                                                                                         ccATATTATCCCCCTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGGACAGACTTTCGAAATTGTCCCGAAACGGGTTCGTTTCGGGACaccg                  >  1:313712/1‑143 (MQ=255)
                                                                                                                 tCCACCTAAAACGGTCTTTCTGTCGGAACCGCACCTGTTTTGTTCCGAAGCGTATGCAGAACAATCTTTTCCGTTGCTAATGCCAGTCAACAGACTTTGGAAATTGTCCCGACACGGGTTCGTTTAGGGACaccgttaccgct          >  2:285737/1‑141 (MQ=255)
                                                                                                                    cccTCAGACGGTATTTCTGGCGGAACCGCACCTGTTGTGTTCTGAAGCGTATTCAGAACAATAGTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTaaa       <  1:348141/143‑1 (MQ=255)
                                                                                                                    cccTAAAACGGTATTTCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTaaa       >  2:163691/1‑143 (MQ=255)
                                                                                                                         agtcGGGATTCCTGTCGGGACCGCAACTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTa  <  1:120412/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
TTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTA  >  NC_000913/1570381‑1570644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: