Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,357,941 C→A 7.9% L1073M (CTG→ATG)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,357,9410CA 7.9% 163.9 / 5.4 51L1073M (CTG→ATG) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base A (2/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

TATGACGGCGGCACCGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATG  >  NC_000913/3357800‑3358088
                                                                                                                                             |                                                                                                                                                   
tATGACGGCGGCACCGGCGCAAGTCCGCTTTCGTCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGCGCTTGTTGAAACCCAGCAGGCTCTAGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATTGCGGCCTTaaatc                                                                                                                                              >  4:385688/1‑147 (MQ=255)
tATGACGGCGGCACCGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAAc                                                                                                                                              <  3:114656/149‑1 (MQ=255)
tATGACGGCGGCACCGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAAc                                                                                                                                              >  3:98890/1‑149 (MQ=255)
      ggcggcACCGGCGCAAGTCCGCCTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGgtgt                                                                                                                                        <  3:359409/149‑1 (MQ=255)
        cggcACCGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCg                                                                                                                                      <  4:458380/149‑1 (MQ=255)
            aCCGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGatat                                                                                                                                  >  4:188034/1‑149 (MQ=255)
              cGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATatca                                                                                                                                <  1:13508/149‑1 (MQ=255)
              cGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATatca                                                                                                                                <  1:499349/149‑1 (MQ=255)
              cGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATatca                                                                                                                                <  2:236466/149‑1 (MQ=255)
               ggCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATatcat                                                                                                                               >  3:145566/1‑149 (MQ=255)
                  gcAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCaa                                                                                                                            <  3:188034/149‑1 (MQ=255)
                          gCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGat                                                                                                                                    <  1:143037/133‑1 (MQ=255)
                          gCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGat                                                                                                                                    >  2:143037/1‑133 (MQ=255)
                               cATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATtctc                                                                                                               <  2:190736/149‑1 (MQ=255)
                                              cAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAgcttc                                                                                                <  4:479358/149‑1 (MQ=255)
                                                        gTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCg                                                                                                              <  2:160376/125‑1 (MQ=255)
                                                        gTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCg                                                                                                              >  1:160376/1‑125 (MQ=255)
                                                         tGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCAc                                                                                     >  1:4548/1‑149 (MQ=255)
                                                              gCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGcc                                                                                <  1:477660/149‑1 (MQ=255)
                                                                  gggCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGAt                                                                            <  1:535388/149‑1 (MQ=255)
                                                                      ttgttgAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGAtggtg                                                                        <  3:184857/149‑1 (MQ=255)
                                                                       tgttgAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGAtggtgg                                                                       >  2:344775/1‑149 (MQ=255)
                                                                        gttgAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGgc                                                                      >  2:205132/1‑149 (MQ=255)
                                                                                ccAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCt                                                              <  2:4548/149‑1 (MQ=255)
                                                                                      ggCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAat                                                        <  2:231884/149‑1 (MQ=255)
                                                                                          ttGGTTGCTAACGGCTTTCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCGGAAAACGGGTGTTGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGACACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCt                                                    <  3:385688/148‑1 (MQ=255)
                                                                                            ggTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTAc                                                  <  4:98890/149‑1 (MQ=255)
                                                                                             gTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAatat                                                      <  2:340489/144‑1 (MQ=255)
                                                                                             gTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAatat                                                      >  1:340489/1‑144 (MQ=255)
                                                                                                    aCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGc                                          >  4:312785/1‑149 (MQ=255)
                                                                                                    aCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGc                                          >  1:5686/1‑149 (MQ=255)
                                                                                                      ggCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCa                                        >  4:281101/1‑149 (MQ=255)
                                                                                                      ggCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTACTTGCCa                                        >  3:166533/1‑149 (MQ=255)
                                                                                                              tCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAgcttcggcttcg                                                                                         >  1:552840/1‑92 (MQ=255)
                                                                                                              tCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAgcttcggcttcg                                                                                         <  2:552840/92‑1 (MQ=255)
                                                                                                              tCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGaaca                                >  2:272374/1‑149 (MQ=255)
                                                                                                                acaAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGaacaac                              >  2:592961/1‑149 (MQ=255)
                                                                                                                 caAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACt                             <  4:166533/149‑1 (MQ=255)
                                                                                                                     aTCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTgcgc                         >  1:505624/1‑149 (MQ=255)
                                                                                                                         gTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACg                     >  3:414060/1‑149 (MQ=255)
                                                                                                                          tttGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACgg                    >  2:248033/1‑149 (MQ=255)
                                                                                                                          tttGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACgg                    <  3:217948/149‑1 (MQ=255)
                                                                                                                             gCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGgtg                 <  1:272374/149‑1 (MQ=255)
                                                                                                                                 gTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGc             >  2:347517/1‑149 (MQ=255)
                                                                                                                                    gATGGCGGCATGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCa                                                                                      <  4:403202/73‑1 (MQ=255)
                                                                                                                                    gATGGCGGCATGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCa                                                                                      <  4:49731/73‑1 (MQ=255)
                                                                                                                                    gATGGCGGCATGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCa                                                                                      >  3:49731/1‑73 (MQ=255)
                                                                                                                                    gATGGCGGCATGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCa                                                                                      >  3:403202/1‑73 (MQ=255)
                                                                                                                                       ggcggcCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAAc          <  4:562857/146‑1 (MQ=255)
                                                                                                                                       ggcggcCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAAc          >  3:562857/1‑146 (MQ=255)
                                                                                                                                          ggcCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCGGATGGTGGCGCTCGGCTGTAAATATCTGCGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGa    >  3:367283/1‑149 (MQ=255)
                                                                                                                                            ccTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATg  >  2:565190/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TATGACGGCGGCACCGGCGCAAGTCCGCTTTCATCGGTGAAATACGCAGGCTGTCCGTGGGAGCTGGGGCTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATG  >  NC_000913/3357800‑3358088

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