Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,882,699 C→A 11.4% intergenic (+30/+18) serA → / ← PP_5156 bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,882,6990CA11.4% 120.9 / 5.9 44intergenic (+30/+18)serA/PP_5156bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/23);  new base A (2/3);  total (18/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 = 8.11e-01

GACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCC  >  NC_002947/5882555‑5882846
                                                                                                                                                |                                                                                                                                                   
gACCAAGGAGAAAGGCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCtggt                                                                                                                                                 <  6:252913/149‑1 (MQ=255)
    aaCGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGcc                                                                                                                                             <  2:270325/149‑1 (MQ=255)
                ggTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAAGAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCttttttt                                                                                                                                  <  4:146782/148‑1 (MQ=255)
                 gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCtttttcata                                                                                                                                <  4:326142/149‑5 (MQ=39)
                    aCGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTAttt                                                                                                                              <  1:336899/148‑1 (MQ=255)
                    aCGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTAttt                                                                                                                              <  6:21201/148‑1 (MQ=255)
                     cGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCa                                                                                                                            <  8:113544/149‑1 (MQ=255)
                     cGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCa                                                                                                                            <  5:126313/149‑1 (MQ=255)
                       tGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttttt                                                                                                                                    <  1:93208/139‑1 (MQ=255)
                       tGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttttt                                                                                                                                    >  2:93208/1‑139 (MQ=255)
                        ggTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCa                                                                                                                            >  1:328147/1‑146 (MQ=255)
                        ggTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCa                                                                                                                            <  2:328147/146‑1 (MQ=255)
                        ggTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACtt                                                                                                                         >  7:244356/1‑149 (MQ=255)
                               gATGTTGAGGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGt                                                                                                                  >  2:320203/1‑149 (MQ=255)
                                        gCCGAGTACTCTGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGCCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATcttctt                                                                                                         <  1:162803/149‑1 (MQ=255)
                                                 tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTtcgactcgac                                                                                                >  5:215687/1‑149 (MQ=255)
                                                 tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCC‑TTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                               >  6:232798/1‑149 (MQ=255)
                                                         ggCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCt                                                                                        >  2:63091/1‑149 (MQ=255)
                                                          gcgcAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATTAAAAAGGGAGGCCCTGGTGTCCTCCCTTTTTTTATTTCACTTCAACCTTGATCTTGTTGGACTCGACGATAGGc                                                                                         <  4:8110/147‑1 (MQ=255)
                                                              aGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCTGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAAt                                                                                   <  8:244356/149‑1 (MQ=255)
                                                                gAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAgg                                                                                          <  6:307087/140‑1 (MQ=255)
                                                                gAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAgg                                                                                          >  5:307087/1‑140 (MQ=255)
                                                                gAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATgg                                                                                 >  6:237129/1‑149 (MQ=255)
                                                                gAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATgg                                                                                 >  7:297583/1‑149 (MQ=255)
                                                                   aaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcac                                                                              <  6:215687/149‑1 (MQ=255)
                                                                     cctgcgccAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                             <  3:152785/141‑1 (MQ=255)
                                                                          aaCAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATAc                                                                       <  1:320203/149‑1 (MQ=255)
                                                                           aCAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACt                                                                      >  3:238668/1‑149 (MQ=255)
                                                                              ggTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGt                                                                   <  4:238668/149‑1 (MQ=255)
                                                                                tCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCa                                                                 <  6:257875/149‑1 (MQ=255)
                                                                                                ttCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTTAATGGCACATGATACTGGTCACCCAGAATCAGTTGCa                                                 <  1:63091/149‑1 (MQ=255)
                                                                                                ttCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTg                                                   >  7:59570/1‑147 (MQ=255)
                                                                                                ttCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCa                                                 <  5:237129/149‑1 (MQ=255)
                                                                                                    cgcgTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGgt                                             >  1:125452/1‑149 (MQ=255)
                                                                                                        tCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGc                                         <  8:59570/149‑1 (MQ=255)
                                                                                                         ccTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCa                                                             >  1:319872/1‑128 (MQ=255)
                                                                                                         ccTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCa                                                             <  2:319872/128‑1 (MQ=255)
                                                                                                                tAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAgg                                               <  6:234280/135‑1 (MQ=255)
                                                                                                                       gCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAgg                          >  3:346358/1‑149 (MQ=255)
                                                                                                                                  aaaaaaGGGAGGCCACCAGGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCggg                                >  7:164350/2‑132 (MQ=39)
                                                                                                                                  aaaaaaGGGAGGCCACCAGGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCggg                                <  8:164350/131‑1 (MQ=39)
                                                                                                                                    aaaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGttt             >  3:175050/1‑149 (MQ=255)
                                                                                                                                     aaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTc            >  8:334585/1‑149 (MQ=255)
                                                                                                                                       aGGGAGGCCCTGGTGGCCTCCC‑TTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGt         <  5:232798/149‑1 (MQ=255)
                                                                                                                                               ccTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGcc  >  2:26943/1‑149 (MQ=255)
                                                                                                                                               ccTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGcc  <  8:297583/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCC  >  NC_002947/5882555‑5882846

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