Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,360,643 T→A 15.4% intergenic (+27/‑533) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,360,6430TA15.4% 72.9 / 6.8 26intergenic (+27/‑533)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (10/12);  new base A (2/2);  total (12/14)
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.86e-01

GACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGA  >  NC_000913/3360496‑3360767
                                                                                                                                                   |                                                                                                                            
gACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACt                                                                                                                              >  1:813909/1‑148 (MQ=255)
       aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCtc                                                                                                                       <  4:781445/148‑2 (MQ=255)
       aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTg                                                                                                                       <  2:813909/148‑1 (MQ=255)
                             gcgATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                             <  2:400137/120‑1 (MQ=255)
                             gcgATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                             >  1:400137/1‑120 (MQ=255)
                                atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                              <  4:89625/132‑1 (MQ=255)
                                atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                              >  3:89625/1‑132 (MQ=255)
                                atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAg                                                                                                 <  1:923303/145‑1 (MQ=255)
                                atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAg                                                                                                 >  2:923303/1‑145 (MQ=255)
                                                  aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGa                                                                                                                   <  1:940360/109‑1 (MQ=255)
                                                  aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGa                                                                                                                   >  2:940360/1‑109 (MQ=255)
                                                         ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtg                                                                                                                           >  4:291268/1‑94 (MQ=255)
                                                         ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtg                                                                                                                           <  3:291268/94‑1 (MQ=255)
                                                         ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                      <  4:450405/99‑1 (MQ=255)
                                                         ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                      >  3:450405/1‑99 (MQ=255)
                                                                   aTTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                             <  3:691560/82‑1 (MQ=255)
                                                                   aTTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                             >  4:691560/1‑82 (MQ=255)
                                                                        cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                  >  1:129991/1‑120 (MQ=255)
                                                                        cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                  <  2:129991/120‑1 (MQ=255)
                                                                                 tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                          >  2:640810/1‑135 (MQ=255)
                                                                                 tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                          <  1:640810/135‑1 (MQ=255)
                                                                                                   tatGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttg                                                                                        >  2:212446/1‑87 (MQ=255)
                                                                                                   tatGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttg                                                                                        <  1:212446/87‑1 (MQ=255)
                                                                                                        aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                 <  2:486211/137‑1 (MQ=255)
                                                                                                        aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                 >  1:486211/1‑137 (MQ=255)
                                                                                                                     ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                             >  3:386492/1‑112 (MQ=255)
                                                                                                                     ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                             <  4:386492/112‑1 (MQ=255)
                                                                                                                                        aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                           >  3:507159/1‑79 (MQ=38)
                                                                                                                                        aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                           <  4:507159/79‑1 (MQ=38)
                                                                                                                                         ttCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTc                                                                                   >  2:110761/1‑54 (MQ=255)
                                                                                                                                         ttCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTc                                                                                   <  1:110761/54‑1 (MQ=255)
                                                                                                                                                 aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtga  <  2:301873/127‑1 (MQ=255)
                                                                                                                                                 aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtga  >  1:301873/1‑127 (MQ=255)
                                                                                                                                                   |                                                                                                                            
GACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGA  >  NC_000913/3360496‑3360767

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