Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,242 G→C 11.9% intergenic (+70/+61) prfC → / ← potF‑I peptide chain release factor 3/putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,013,2420GC11.9% 79.5 / 5.5 34intergenic (+70/+61)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (12/17);  new base C (2/2);  total (14/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 = 5.86e-01

CGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGC  >  NC_002947/1013096‑1013388
                                                                                                                                                  |                                                                                                                                                  
cGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                                  <  4:368822/149‑1 (MQ=255)
cGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                                  <  1:202504/149‑1 (MQ=255)
cGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                                  <  6:359610/149‑1 (MQ=255)
      gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                            <  1:350971/149‑1 (MQ=255)
       tttCAACCTGGCTCTTATGGAATATCGCTGGCCGGACATAAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCTGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCTCTCCCACAGGAATCACCAGACCTGGATGTCAtgt                                                                                                                                           <  3:317037/147‑1 (MQ=38)
       ggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                           <  3:239082/149‑1 (MQ=255)
       ggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                           <  5:21820/149‑1 (MQ=255)
                 gcgcTGATGGAAGAGCGCTGGGCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCTGCCCTTCGCGGGAACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                                  <  7:254985/148‑1 (MQ=21)
                 gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                                  <  5:44182/148‑1 (MQ=255)
                       aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCg                                                                                                                           >  3:122444/1‑149 (MQ=255)
                                cgcTGGCAGGACATCACGTTCCGGTCTACCCGCGAGCACAATTAAGGCCTGAACTCGGGGCTGCTTTGCAGCCCTTCGCGGTCACGCCCGCTCCCACAGGAATCACCCCACCTGGAGTTCATGTGATACCTGTGGGAGGGGGCGTGccc                                                                                                                  <  4:72636/149‑1 (MQ=37)
                                 gcTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCcg                                                                                                                 <  5:119829/149‑1 (MQ=255)
                                            aTCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGc                                                                                                      <  4:122444/149‑1 (MQ=255)
                                              cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                            >  3:324477/1‑109 (MQ=255)
                                              cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                            <  4:324477/109‑1 (MQ=255)
                                                            cccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCAtt                                                                                        <  5:317722/147‑1 (MQ=255)
                                                            cccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCCCACGAATAACCACACCGGGAGGTCATGTGGTACCTGTGGGGGCGGGCGTGCCCGCGAAGGGCCCCGAGGCGGGCCCAttt                                                                                       >  1:145637/1‑148 (MQ=37)
                                                                       cATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGCCACGCCCGCACCCACACGAATCACCACACCAAGAGCACAAGTGATACCAGAGGGAGCGGGAGTGCCCGCGAAAGGCCCCTAGGCGCCCCCATTTAAATTTACGGca                                                                           >  8:259703/1‑148 (MQ=25)
                                                                                 tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCgg                                                                                                                          >  1:201164/1‑92 (MQ=255)
                                                                                 tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCgg                                                                                                                          <  2:201164/92‑1 (MQ=255)
                                                                                             tgctTTGAAGCCCGTCGCGGGCACGCCCGCTACCACATGACTAACCAGACCTGTAGTTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgcg                                                     <  5:205170/149‑1 (MQ=37)
                                                                                                         cTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGt                                         <  7:259703/149‑1 (MQ=255)
                                                                                                             gcgttcACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAGCTGTGTGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGAt                                     <  5:183855/144‑1 (MQ=255)
                                                                                                                              cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGt                                         >  2:89022/1‑128 (MQ=255)
                                                                                                                              cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGt                                         <  1:89022/128‑1 (MQ=255)
                                                                                                                                caGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGcc                              >  5:126603/1‑137 (MQ=21)
                                                                                                                                caGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGcc                              <  6:126603/137‑1 (MQ=21)
                                                                                                                                      tcacatgACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCa            >  3:227163/7‑149 (MQ=37)
                                                                                                                                        aCCAGACCTGGAGGTCATGTGATACCTGTGGGATCGGGCGTGAACGCGAAGTGCCGCTAGGCTGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGACTACGGTCGATCTTCGCCGGCATCGCCTCCATCGCATa          >  7:71441/1‑149 (MQ=255)
                                                                                                                                         ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCTAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa          >  1:323801/1‑148 (MQ=255)
                                                                                                                                         ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa          >  8:32697/1‑148 (MQ=255)
                                                                                                                                         ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa          >  8:26731/1‑148 (MQ=255)
                                                                                                                                         ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTACGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa          >  3:21699/1‑148 (MQ=255)
                                                                                                                                          cAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa          >  2:124769/1‑147 (MQ=255)
                                                                                                                                                  gAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGc  <  7:32697/147‑1 (MQ=255)
                                                                                                                                                  gAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGc  <  7:26731/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGC  >  NC_002947/1013096‑1013388

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