Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,505,814 C→T 22.2% A69A (GCG→GCA yhfS ← FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,505,8140CT22.2% 77.1 / 16.0 36A69A (GCG→GCAyhfSFNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base C (15/13);  new base T (4/4);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATC  >  NC_000913/3505675‑3505954
                                                                                                                                           |                                                                                                                                            
taaGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCt                                                                                                                                    <  1:52002/150‑1 (MQ=255)
taaGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCt                                                                                                                                    <  1:237684/150‑1 (MQ=255)
         gCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGaa                                                                                                                           >  2:188152/1‑150 (MQ=255)
         gCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGaa                                                                                                                           >  2:157384/1‑150 (MQ=255)
         gCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGaa                                                                                                                           >  2:157380/1‑150 (MQ=255)
         gCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGaa                                                                                                                           >  2:157378/1‑150 (MQ=255)
          ccccATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGaaa                                                                                                                          >  2:120371/1‑150 (MQ=255)
          ccccATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGaaa                                                                                                                          >  2:120373/1‑150 (MQ=255)
          ccccATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGaaa                                                                                                                          >  2:120374/1‑150 (MQ=255)
               tCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCAt                                                                                                                     >  1:7796/1‑150 (MQ=255)
               tCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCAt                                                                                                                     >  1:7798/1‑150 (MQ=255)
                    tCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATcagcc                                                                                                                >  1:214876/1‑150 (MQ=255)
                             aCCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGt                                                                                                       >  1:117211/1‑150 (MQ=255)
                               ccGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTc                                                                                                     <  2:7796/150‑1 (MQ=255)
                               ccGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTc                                                                                                     <  2:7798/150‑1 (MQ=255)
                                  tgtCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCcac                                                                                                  <  2:242616/150‑1 (MQ=255)
                                  tgtCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCcac                                                                                                  <  1:137786/150‑1 (MQ=255)
                                  tgtCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCcac                                                                                                  <  1:137770/150‑1 (MQ=255)
                                                  ggCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGGGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAAcgcgc                                                                                  >  2:96338/1‑150 (MQ=255)
                                                            gCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCa                                                                        <  1:157378/150‑1 (MQ=255)
                                                            gCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCa                                                                        <  1:188152/150‑1 (MQ=255)
                                                            gCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCa                                                                        <  1:157384/150‑1 (MQ=255)
                                                            gCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGTGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCa                                                                        <  1:157380/150‑1 (MQ=255)
                                                              aCCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGACTCGCGCCAGTCCCCGCGCCCTGCACCCGCGCCGCAGCCGGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAgn                                                                      <  1:152875/150‑2 (MQ=255)
                                                              aCCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAgn                                                                      <  1:152877/150‑2 (MQ=255)
                                                                              ccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTCAGCCTa                                                      >  1:72868/1‑150 (MQ=255)
                                                                              ccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTa                                                      >  1:57738/1‑150 (MQ=255)
                                                                                                          gcgAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGcc                          >  1:131330/1‑150 (MQ=255)
                                                                                                          gcgAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGcc                          >  1:131313/1‑150 (MQ=255)
                                                                                                                     cAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGc               <  2:210349/150‑1 (MQ=255)
                                                                                                                        tCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGAc            >  1:7731/1‑150 (MQ=255)
                                                                                                                           ccGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaaa         >  2:221168/1‑150 (MQ=255)
                                                                                                                           ccGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaaa         <  2:221291/150‑1 (MQ=255)
                                                                                                                            cgcgCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAAt        <  2:214876/150‑1 (MQ=255)
                                                                                                                                cccTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTa    >  2:170373/1‑150 (MQ=255)
                                                                                                                                  cTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATc  <  1:237945/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
TAAGCGTCAGCCCCATCTGCTCAATAATAACCCGTGTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATC  >  NC_000913/3505675‑3505954

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