Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,570,600 T→C intergenic (‑111/+45) 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,6000TC97.1% 113.5 / NA 36intergenic (‑111/+45)gadC/gadBL‑glutamate:4‑aminobutyrate antiporter/glutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (15/19);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTC  >  NC_000913/1570465‑1570731
                                                                                                                                       |                                                                                                                                   
tACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCg                                                                                                                                 <  2:492215/140‑1 (MQ=255)
   ctgtctgtACTGATGTAGCCATTTTTTCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGaaa                                                                                                                              <  1:245899/140‑1 (MQ=255)
         gtACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGtt                                                                                                                        <  1:69683/140‑1 (MQ=255)
                  tAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAACGCCAGTGATCAGGCTTCGGAAATTGTCCCGAAACGGGTTCGTTTCggg                                                                                                               >  2:236210/1‑140 (MQ=255)
                    gCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGac                                                                                                             >  2:616039/1‑140 (MQ=255)
                      cATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGacac                                                                                                           >  1:7204/1‑140 (MQ=255)
                           tatCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACaccgtta                                                                                                      <  1:295890/140‑1 (MQ=255)
                                  cTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTaa                                                                                               >  1:642839/1‑140 (MQ=255)
                                  cTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTaa                                                                                               <  2:448762/140‑1 (MQ=255)
                                   tAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTaa                                                                                               >  1:154361/1‑139 (MQ=255)
                                   tAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTaa                                                                                               <  2:154361/139‑1 (MQ=255)
                                           tATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCa                                                                                      >  1:57503/1‑140 (MQ=255)
                                             ttCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAgg                                                                                    <  2:475164/140‑1 (MQ=255)
                                              tCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGt                                                                                   <  1:509014/140‑1 (MQ=255)
                                              tCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGt                                                                                   <  1:441248/140‑1 (MQ=255)
                                                  gTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGt                                                                               >  2:381720/1‑140 (MQ=255)
                                                   tCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGtt                                                                              <  2:388137/140‑1 (MQ=255)
                                                           gCACCTGTTTTGTTCTGAAGCGTATTCAGCACAATATCTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctg                                                                      <  1:226935/140‑1 (MQ=255)
                                                            cACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctgt                                                                     >  1:8230/1‑140 (MQ=255)
                                                                tGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctgttctg                                                                 <  1:381720/140‑1 (MQ=255)
                                                                          tGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATa                                                       <  1:535391/140‑1 (MQ=255)
                                                                            aaGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATAcc                                                     >  2:26543/1‑140 (MQ=255)
                                                                              gCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCt                                                   <  2:7204/140‑1 (MQ=255)
                                                                                         aCAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCggg                                        >  2:544105/1‑140 (MQ=255)
                                                                                           aaTATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTg                                      >  1:599539/1‑140 (MQ=255)
                                                                                                 ttCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCt                                <  1:538048/140‑1 (MQ=255)
                                                                                                         cTAATGCCAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGCGCAATACTCTGAAGGTTCGGGCGATCGCTGAGATAtt                        >  1:209228/1‑140 (MQ=255)
                                                                                                               ccAGTGAATAGACTTTGGAAATTGCTCCGAACCGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGACAATGCCCTGCAGTTTCAGGTGATCGCTGAGATATTTCAggg                  <  2:480671/140‑1 (MQ=25)
                                                                                                               ccAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAggg                  <  1:538173/140‑1 (MQ=255)
                                                                                                               ccAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAggg                  <  1:616039/140‑1 (MQ=255)
                                                                                                               ccAGTGAACAGACTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAggg                  <  2:642839/140‑1 (MQ=255)
                                                                                                                 aGTGAACAGACTTTGGAAATTGACCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAg                <  2:271311/140‑1 (MQ=255)
                                                                                                                          aCTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAg       >  2:438501/1‑140 (MQ=255)
                                                                                                                          aCTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCAGATGATCGCTGAAATATTTCAGGGAGGCTTTGTAg       >  2:296698/1‑140 (MQ=255)
                                                                                                                           cTTTGGAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGt      <  1:666803/140‑1 (MQ=255)
                                                                                                                               ggAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTc  >  2:537150/1‑140 (MQ=255)
                                                                                                                               ggAAATTGCCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGGGATCGCTGAGATATTTCAGGGCGGCTTTGTAGTCTTc  >  1:299995/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                   
TACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTC  >  NC_000913/1570465‑1570731

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