Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,951,644 G→A 10.8% P633S (CCG→TCG)  mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,951,6440GA10.8% 95.5 / 3.2 37P633S (CCG→TCG) mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base A (2/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

ACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCAT  >  NC_002947/1951509‑1951787
                                                                                                                                       |                                                                                                                                               
aCCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCg                                                                                                                                      <  7:62510/147‑1 (MQ=255)
   tATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATg                                                                                                                                            >  6:18006/1‑138 (MQ=255)
   tATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATg                                                                                                                                            <  5:18006/138‑1 (MQ=255)
       tgCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCca                                                                                                                              <  2:182450/148‑1 (MQ=255)
        gCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCTGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCcaca                                                                                                                            <  6:257972/149‑1 (MQ=255)
         cTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCGGCTAAAACCGCCCAAAACAAAACGGGTTTGGATGCCAGCCATCCGCCCATCCCGCCGGCATTGGCACCGGTACCCcaaa                                                                                                                            >  7:200480/1‑146 (MQ=255)
                agagTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAgg                                                                                                                    >  2:281413/1‑149 (MQ=255)
                              tttGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTc                                                                                                       <  6:182345/148‑1 (MQ=255)
                                                ggcaggcaAGCGGATGCAATCATTTCCCCAGGATCAACCGGCGCACTCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCg                                                                                      <  3:124092/147‑1 (MQ=255)
                                                 gaaggcaACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAACCCGTTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGa                                                                                   <  1:336253/147‑1 (MQ=255)
                                                  caggcaACCCGCTGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGAt                                                                                  <  6:26423/149‑1 (MQ=255)
                                                        aCCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAAccc                                                                            <  1:119957/149‑1 (MQ=255)
                                                          ccGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCcgc                                                                          <  6:250104/149‑1 (MQ=255)
                                                                 aaTCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGc                                                                    >  7:315442/1‑148 (MQ=255)
                                                                    cATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGc                                                                                        >  2:155290/1‑125 (MQ=255)
                                                                    cATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGc                                                                                        <  1:155290/125‑1 (MQ=255)
                                                                         cccccGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCGCACATTCCCAGGCCGCCACCAGCTAGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGc                                                           >  5:89298/1‑149 (MQ=255)
                                                                             cGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGc                                                               >  8:39418/1‑141 (MQ=255)
                                                                             cGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGc                                                               <  7:39418/141‑1 (MQ=255)
                                                                             cGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGtt                                                       <  1:281413/149‑1 (MQ=255)
                                                                                       ggCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAAc                                             <  6:18342/149‑1 (MQ=255)
                                                                                        ttgccaTCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGATGTTCACATTCAAcc                                            >  6:269941/6‑149 (MQ=255)
                                                                                               ccAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTccc                                                                                                                                <  4:171051/58‑1 (MQ=255)
                                                                                               ccAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTccc                                                                                                                                >  3:171051/1‑58 (MQ=255)
                                                                                                 aaaCCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGc                                    >  3:284807/1‑148 (MQ=255)
                                                                                                   aCCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCaa                                 <  1:204908/149‑1 (MQ=255)
                                                                                                           ggATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCa                                      <  8:161134/136‑1 (MQ=255)
                                                                                                           ggATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCa                                      >  7:161134/1‑136 (MQ=255)
                                                                                                              tGGAAAGCCTGCTCCACTGCCCGCGGCAATGGCAGCTGTTCCCCACATACCCAGCCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTTTTCACATACAACCCGCCCAGCCAAGGGCAGGGcac                      <  6:89298/149‑1 (MQ=255)
                                                                                                                            cgcTGCCCGCGACAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcagg                          >  4:249307/3‑131 (MQ=255)
                                                                                                                            cgcTGCCCGCGACAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcagg                          >  4:244806/3‑131 (MQ=255)
                                                                                                                            cgcTGCCCGCGACAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcagg                          <  3:249307/129‑1 (MQ=255)
                                                                                                                            cgcTGCCCGCGACAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcagg                          <  3:244806/129‑1 (MQ=255)
                                                                                                                              cTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCcgc      >  3:233737/1‑149 (MQ=255)
                                                                                                                                 cccGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCa   >  4:218655/1‑149 (MQ=255)
                                                                                                                                  ccGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCAt  >  1:49422/1‑149 (MQ=255)
                                                                                                                                  ccGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCAt  >  8:54661/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                               
ACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCGGTTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGCGCAT  >  NC_002947/1951509‑1951787

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: