Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,122,117 G→A 8.5% T563T (ACC→ACT)  glcB ← malate synthase G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,122,1170GA 8.5% 135.2 / 3.1 47T563T (ACC→ACT) glcBmalate synthase G
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (2/2);  total (26/21)
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.12e-01

CTCCACCCAGCGCACCACGTACCCCAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTGGCTGTACATGTCTGCCATCA  >  NC_000913/3121970‑3122265
                                                                                                                                                   |                                                                                                                                                    
cTCCACCCAGCGCACCACGTACCCCAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCg                                                                                                                                                      <  1:100141/148‑1 (MQ=255)
    ccccAGCGCACCACGTACCCCAGAATCCCCTGCACGTTGGTATCCAGCTCTTGTTGGAGCTCTTGCGCCGACCAGTTAGCGTTTTCAGCCACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCAGTCTg                                                                                                                                                 >  7:111817/2‑149 (MQ=255)
    aCCCAGCGCACCACGTACCCCAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCAGTCTg                                                                                                                                                 >  7:111815/1‑149 (MQ=255)
         gcgcACCACGTACCCCAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGAACTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCaa                                                                                                                                            <  5:29321/149‑1 (MQ=255)
                  gTACCCCAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCtt                                                                                                                                   <  2:109065/149‑1 (MQ=255)
                        cAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGc                                                                                                                             >  7:226454/1‑149 (MQ=255)
                        cAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGc                                                                                                                             >  5:184641/1‑149 (MQ=255)
                              ccccTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTAGCGCCGACAAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTCGGCAATGTTGGCTTGTACGCTCTGt                                                                                                                        >  3:104257/1‑148 (MQ=255)
                                ccTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACg                                                                                                                     >  3:57143/1‑149 (MQ=255)
                                   gCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTa                                                                                                                       >  4:78300/1‑144 (MQ=255)
                                   gCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTa                                                                                                                       <  3:78300/144‑1 (MQ=255)
                                      cGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGttggtt                                                                                                               <  3:231728/149‑1 (MQ=255)
                                       gttgttATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGttggttt                                                                                                              >  3:183272/1‑149 (MQ=255)
                                             aTCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTtggtgg                                                                                                         >  5:41340/1‑148 (MQ=255)
                                                   ctctTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTAGGTCTg                                                                                                                                                 <  3:172533/102‑1 (MQ=255)
                                                   ctctTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTAGGTCTg                                                                                                                                                 >  4:172533/1‑102 (MQ=255)
                                                           cgttctcTTGCGCCGACCAGTTAGAGTTTTCAGCAACCTGATTAGTCAGCAGATCGTCCAGCAGCGGTTCATATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATg                                                                                           <  6:10792/145‑1 (MQ=255)
                                                                tctTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCg                                                                                      >  7:77314/1‑148 (MQ=255)
                                                                tctTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGt                                                                                     <  4:57143/149‑1 (MQ=255)
                                                                  ttGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCACCAGCGGTAAAAATTCAGCATTGAACTCGGTCTGGGCAATGTAGGCTTGTACGCTCTGTACGTTGCTTTGGTGGTAGTGCAGCGCATg                                                                                           >  5:10792/1‑141 (MQ=255)
                                                                    gcgcCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTag                                                                                   >  7:39582/1‑147 (MQ=255)
                                                                    gcgcCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTag                                                                                   >  7:35614/1‑147 (MQ=255)
                                                                    gcgcCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTa                                                                                    <  6:184641/146‑1 (MQ=255)
                                                                      gcCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGttg                                                                                                                  >  6:48538/1‑114 (MQ=255)
                                                                      gcCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGttg                                                                                                                  <  5:48538/114‑1 (MQ=255)
                                                                        cGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGCACTCAGTCTGGGCAATGTTGGCTTGTCCGCGCTGTACGTTGGTTTGGTGGTCGTGCAGCGCATGGAGCGTAGCAGCgg                                                                             <  8:111817/149‑1 (MQ=255)
                                                                        cGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCAGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCgg                                                                             <  8:111815/149‑1 (MQ=255)
                                                                            cAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCATATCGTCCAGCAGAGGTTCAAATTCAGAATTGAACTCGGTCTTGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGATCAGCTCATGGAGCGTAGCAGCGGTTgg                                                                         >  8:129551/1‑149 (MQ=255)
                                                                               ttAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGa                                                                      >  6:90825/1‑149 (MQ=255)
                                                                               ttAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGa                                                                      >  7:71150/1‑149 (MQ=255)
                                                                                 aGCGTTTTCAGCAACCGAAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGAc                                                                     <  1:206578/148‑1 (MQ=255)
                                                                                    ggtgtCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGaa                                                                 <  2:18076/145‑1 (MQ=255)
                                                                                        tCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAAccc                                                              <  5:90825/148‑1 (MQ=255)
                                                                                          aGCAACCGGAATAGTCAGCAGCTCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAgg                                                           <  4:104257/149‑1 (MQ=255)
                                                                                              aCCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGGGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCtgt                                                       <  3:162208/149‑1 (MQ=255)
                                                                                                  gAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCATCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGAAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGttt                                                   <  5:196030/149‑1 (MQ=255)
                                                                                                  gAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGttt                                                   >  5:91295/1‑149 (MQ=255)
                                                                                                       gTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGccc                                               <  6:91295/148‑1 (MQ=255)
                                                                                                          agcagATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGc                                           >  4:126876/1‑149 (MQ=255)
                                                                                                            cagATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCAc                                         >  3:98291/1‑149 (MQ=255)
                                                                                                             agATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACg                                        <  8:39582/149‑1 (MQ=255)
                                                                                                               aTCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACg                                        <  8:35614/147‑1 (MQ=255)
                                                                                                                       gcagcGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCg                              >  1:85961/1‑149 (MQ=255)
                                                                                                                       gcagcGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCACCGGCACGCAGTTGGTCg                              >  4:155799/1‑149 (MQ=255)
                                                                                                                               tCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTg                      >  1:109811/1‑149 (MQ=255)
                                                                                                                                        gCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTGGCTGTACAt             <  8:71150/149‑1 (MQ=255)
                                                                                                                                          aTTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTGGCTGTACATGt           >  4:174680/1‑149 (MQ=255)
                                                                                                                                          aTTGAACTCGGTCTGGGCAATGCTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTGGCTGTACATGt           >  3:147875/1‑149 (MQ=255)
                                                                                                                                                   ggTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTTGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGTTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTGGCTGTACATGTCTGCCATCa  <  7:129551/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCCACCCAGCGCACCACGTACCCCAGAATCCCCTGCACGTTGTTATCCAGCTCTTGTTGGATCTCTTGCGCCGACCAGTTAGCGTTTTCAGCAACCGGAATAGTCAGCAGATCGTCCAGCAGCGGTTCAAATTCAGCATTGAACTCGGTCTGGGCAATGTTGGCTTGTACGCTCTGTACGTTGGTTTGGTGGTAGTGCAGCGCATGGAGCGTAGCAGCGGTTGGTGACGGAACCCAGGCTGTGTTTGCCCCGGCACGCAGTTGGTCGCCCTTCTGGCTGTACATGTCTGCCATCA  >  NC_000913/3121970‑3122265

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