Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,372,641 C→A 10.8% S156* (TCA→TAA)  gdhB → NAD‑specific glutamate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,372,6410CA10.8% 122.6 / 3.0 37S156* (TCA→TAA) gdhBNAD‑specific glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base A (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

ACCGAGCTCAACCGCCGTGGCTACAGCATCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAGGCCAAGC  >  NC_002947/2372508‑2372769
                                                                                                                                     |                                                                                                                                
cccGAGCTCAACCGCCGTGGCTACAGCATCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCAc                                                                                                                                <  2:534459/135‑1 (MQ=255)
      ttgaaCCGCCGTGGCGACAGCATCCACACCCTGCATAGCTGCGTGCTCTGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGt                                                                                                                          <  1:265201/133‑1 (MQ=255)
      cTCAACCGCCGTGGCTACAGCATCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGt                                                                                                                          <  2:617455/136‑1 (MQ=255)
        cAACCGCCGTGGCTACAGCATCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTAc                                                                                                                        >  2:6939/1‑136 (MQ=255)
        cAACCGCCGTGGCTACAGCATCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTAc                                                                                                                        >  1:244286/1‑136 (MQ=255)
                           aTCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTgcg                                                                                                     >  2:539777/1‑136 (MQ=255)
                              cacaCCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCa                                                                                                  >  2:610696/1‑136 (MQ=255)
                              cacaCCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCa                                                                                                  >  2:604352/1‑136 (MQ=255)
                               acacCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCaa                                                                                                 >  1:180901/1‑136 (MQ=255)
                               acacCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCaa                                                                                                 >  2:168303/1‑136 (MQ=255)
                                  cccTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATgc                                                                                              <  1:87682/136‑1 (MQ=255)
                                      gCAACCCACCGTGCGCAGCGGGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATgccgcc                                                                                          <  1:609397/136‑1 (MQ=255)
                                          accaccGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGc                                                                                      >  2:396040/1‑136 (MQ=255)
                                          accaccGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTAACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGc                                                                                      >  2:109837/1‑136 (MQ=255)
                                          accaccGTGCTCAGCGTGCGCAGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTAACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGc                                                                                      >  2:432928/1‑136 (MQ=255)
                                                gTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGCGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAg                                                                                >  2:484758/1‑136 (MQ=255)
                                                       gCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGAc                                                                         <  1:512367/136‑1 (MQ=255)
                                                               cGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGGCCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCgaga                                                                 >  2:312680/1‑136 (MQ=255)
                                                                gTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGAt                                                                <  1:647424/136‑1 (MQ=255)
                                                                     gcCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGc                                                           >  2:488418/1‑136 (MQ=255)
                                                                            gCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCt                                                    <  1:6939/136‑1 (MQ=255)
                                                                            gCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCt                                                    <  1:49143/136‑1 (MQ=255)
                                                                                  tGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGa                                              <  1:396040/136‑1 (MQ=255)
                                                                                   gTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAg                                             >  1:218061/1‑136 (MQ=255)
                                                                                                 ccAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGCCCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGggtggtggt                               <  1:312680/136‑1 (MQ=255)
                                                                                                                gCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAAcc                <  2:10226/136‑1 (MQ=255)
                                                                                                                  gAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTAc                                         <  2:527083/109‑1 (MQ=255)
                                                                                                                  gAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTAc                                         >  1:527083/1‑109 (MQ=255)
                                                                                                                  gAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGa              >  1:592267/1‑136 (MQ=255)
                                                                                                                      gCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGaa          >  1:4382/1‑136 (MQ=255)
                                                                                                                       cGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAg         >  1:669202/1‑136 (MQ=255)
                                                                                                                        gTGCGCTACGAGTAACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAgg        <  1:109837/136‑1 (MQ=255)
                                                                                                                        gTGCGCTACGAGTAACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAgg        <  1:432928/136‑1 (MQ=255)
                                                                                                                             cTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAGGCCAAg   <  1:99912/136‑1 (MQ=255)
                                                                                                                             cTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAGGCCAAg   <  1:610696/136‑1 (MQ=255)
                                                                                                                             cTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAGGCCAAg   <  1:604352/136‑1 (MQ=255)
                                                                                                                              tACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAGGCCAAGc  >  2:110850/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                
ACCGAGCTCAACCGCCGTGGCTACAGCATCCACACCCTGCAAACCACCGTGCTCAGCGTGCGCCGTGGCGCCAAGGGCGAGCTGTTGGAGCTGCTGCCCAAGGGAACCCAGGGCGAAGGCGTGCGCTACGAGTCACTGATGTACCTGGAGATCGACCGCTGCGCCAATGCCGCCGAGCTGACAGTGCTGACCCGCGAGATCGAGCAGGTGCTGGCCGAGGTACGGGTGGTGGTGGCCGACTTCGAACCGATGAAGGCCAAGC  >  NC_002947/2372508‑2372769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑