Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,174 G→T 10.7% G185C (GGC→TGC)  corC → Co(2+)/Mg(2+)efflux heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,451,1740GT10.7% 124.1 / 4.3 44G185C (GGC→TGC) corCCo(2+)/Mg(2+)efflux heat shock protein
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base T (2/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGAGTTCTCCGATGACGAGTTCGACAC  >  NC_002947/5451034‑5451319
                                                                                                                                            |                                                                                                                                                 
ccTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCATATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATc                                                                                                                                           <  5:198485/149‑1 (MQ=255)
 cTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCTTCGATATCg                                                                                                                                          >  4:157558/1‑149 (MQ=255)
    cGTGCCCGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGaa                                                                                                                                        <  6:201528/148‑1 (MQ=255)
         ccGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGATGAc                                                                                                                                     >  3:47521/1‑146 (MQ=255)
          cGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCTTCGATATCGCcgacgatc                                                                                                                                 <  3:157558/149‑3 (MQ=255)
             gTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCGGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCac                                                                                                                               >  7:114704/1‑148 (MQ=255)
                cAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGa                                                                                                                                      >  4:176096/1‑138 (MQ=255)
                cAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGa                                                                                                                                      <  3:176096/138‑1 (MQ=255)
                 aaGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAg                                                                                                                                  <  6:1979/141‑1 (MQ=255)
                 aaGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAg                                                                                                                                  >  5:1979/1‑141 (MQ=255)
                      ccTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAgg                                                                                                                     <  1:10296/149‑1 (MQ=255)
                       cTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGa                                                                                                                                      <  4:335741/131‑1 (MQ=255)
                       cTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGa                                                                                                                                      >  3:335741/1‑131 (MQ=255)
                                ttGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCac                                                                                                                               <  3:115920/129‑1 (MQ=255)
                                ttGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCac                                                                                                                               >  4:115920/1‑129 (MQ=255)
                                             tCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATGGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAAcaaatg                                                                                              >  6:189176/1‑144 (MQ=255)
                                                      accacaaccacaTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCt                                                                                                            <  2:348900/126‑1 (MQ=255)
                                                      accacaaccacaTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCt                                                                                                            >  1:348900/1‑126 (MQ=255)
                                                      accacaaccacaTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGt                                                                                     >  5:331391/1‑149 (MQ=255)
                                                      accacaaccacaTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGt                                                                                     >  2:184316/1‑149 (MQ=255)
                                                      accacaaccacaTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGt                                                                                     <  1:184316/149‑1 (MQ=255)
                                                        cacaaccacaTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGa                                                                                   >  8:61450/1‑149 (MQ=255)
                                                                caTGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTCCCTgg                                                                           <  4:47521/149‑1 (MQ=255)
                                                                     ccATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCaaa                                                                      <  8:114704/149‑1 (MQ=255)
                                                                       aTCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAgc                                                                    >  4:329254/1‑149 (MQ=255)
                                                                        tCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCg                                                                                       <  3:24750/129‑1 (MQ=255)
                                                                        tCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCg                                                                                       >  4:24750/1‑129 (MQ=255)
                                                                          gTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCt                                                                 <  7:61450/149‑1 (MQ=255)
                                                                           tCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCt                                                                 >  3:314347/1‑148 (MQ=255)
                                                                             atagacgaATATGGCGGCGTGGCGGGCCTGGTCAACATTGAGGACGTGCTGGAGAAGATCGTCGGCGATATCGAAGACGAGCACGACGACGAGGAAGACAGCAACATCAAACCACACCCCAGCGTGGATATCAAGGTCAAACCGCATAc                                                              >  2:89882/4‑149 (MQ=255)
                                                                              tCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTAcc                                                             >  5:168527/1‑149 (MQ=255)
                                                                                 acgaATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCTGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTAcccc                                                           >  8:231637/1‑148 (MQ=255)
                                                                                              cgTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCa                                             <  6:168527/149‑1 (MQ=255)
                                                                                                           gTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGttcttc                                  <  3:329254/147‑1 (MQ=255)
                                                                                                               ccATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGc                            >  6:51573/1‑149 (MQ=255)
                                                                                                                cATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGc                            >  3:151646/1‑148 (MQ=255)
                                                                                                                cATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGc                            >  3:156830/1‑148 (MQ=255)
                                                                                                                 aTTGAGGACGTGCTGGAGCAGATCGTCTGCGATATCGAAGACGAGCACGACGTCGAGGAAGAAAGCTACATCAAACCACTGCCCAGCGTTGATTTCCTGGTCAAAGAGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGa                          <  1:89882/149‑1 (MQ=255)
                                                                                                                  ttGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGAAAGCTACATCAAACCACTGCCAAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGAg                         >  7:331265/1‑149 (MQ=255)
                                                                                                                      ggACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCaa                                            >  2:31741/1‑126 (MQ=255)
                                                                                                                      ggACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGAGTtct                     <  8:331265/149‑1 (MQ=255)
                                                                                                                      ggACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGACAAAGCGCTTACCCCGATCGAGAACTTCaa                                            <  1:31741/126‑1 (MQ=255)
                                                                                                                           tGCTGGAGCAGATCGTCTGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGAGTTCTCCGAt                <  7:231637/149‑1 (MQ=255)
                                                                                                                                 aGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGAGTTCTCCGATGACGAg          <  6:331391/149‑1 (MQ=255)
                                                                                                                                         gtcgGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCGTCGACAGCGAGTTATCCGATGACGAGTTCGacac  >  7:128333/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CCTTCGTGCCCGAGTCCAAGCGCCTGAACGTGTTGCTGCGTGAATTCCGCGCCAACCACAACCACATGGCCATCGTCATCGACGAATATGGCGGCGTGGCGGGCCTGGTCACCATTGAGGACGTGCTGGAGCAGATCGTCGGCGATATCGAAGACGAGCACGACGTCGAGGAAGACAGCTACATCAAACCACTGCCCAGCGGTGATTTCCTGGTCAAAGCGCTTACCCCGATCGAGAACTTCAACGAGTTCTTCGACAGCGAGTTCTCCGATGACGAGTTCGACAC  >  NC_002947/5451034‑5451319

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: ‑