Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,194 T→A 13.9% intergenic (+22/+109) 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,1940TA13.9% 89.7 / 6.1 36intergenic (+22/+109)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (12/19);  new base A (2/3);  total (14/22)
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.99e-01

GAAAACTTGGCCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGA  >  NC_002947/1013050‑1013331
                                                                                                                                                |                                                                                                                                         
gAAAACTTGGCCATCGACGGCGGCGGTCACCAGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAACAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAACGCCAGAAATCGGGGCTGCATTGCAg                                                                                                                                       >  4:99356/1‑149 (MQ=255)
  aaaCTTGGCCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGcc                                                                                                                                     <  1:224065/149‑1 (MQ=255)
    aCTTGGCCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCAAAGCAGcccc                                                                                                                                   <  1:66320/149‑2 (MQ=255)
        ggCCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTcg                                                                                                                                <  1:217821/148‑1 (MQ=255)
        ggCCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCAAAGCAGccccga                                                                                                                                 <  4:11547/147‑4 (MQ=255)
          ccATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                             <  5:83055/149‑1 (MQ=255)
          ccATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCgg                                                                                                                             <  6:332926/149‑1 (MQ=255)
                     ggcggTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                           <  2:354661/140‑1 (MQ=255)
                     ggcggTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCaagatcgg                                                                                                                   >  1:354661/1‑141 (MQ=255)
                               tGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGaa                                                                                                        >  7:166987/1‑149 (MQ=255)
                               tGACCTACCTGGCGCCGACCCGGATCAACCTGGCGCTGATGCAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCATCGCGGGCACGCCCGCTCCCACAggta                                                                                                        >  1:65317/1‑147 (MQ=255)
                                   ctacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCAc                                                                                                    <  6:51493/149‑1 (MQ=255)
                                     acctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                    >  4:135694/1‑147 (MQ=255)
                                      cctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                  <  3:135694/148‑1 (MQ=255)
                                       ctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGa                                                                                                >  4:60381/1‑149 (MQ=255)
                                        tGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGAc                                                                                               >  5:124450/1‑149 (MQ=255)
                                            gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATGAAGTTCCGCGCTACCCGCGAGCACCATTAAGGGCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCtcc                                                                                           <  3:99356/149‑3 (MQ=38)
                                            gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTgg                                                                                           <  8:166987/149‑1 (MQ=255)
                                                     ggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                  <  6:20262/149‑1 (MQ=255)
                                                          aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                             <  7:87648/149‑1 (MQ=255)
                                                                cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                       <  4:352698/149‑1 (MQ=255)
                                                                cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                       <  6:124450/149‑1 (MQ=255)
                                                                        gAAGAGCGCTGGCCGGACATCAGGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGACATGTGATACCTGTGGGAGCGGGc                                                               <  7:52660/149‑1 (MQ=255)
                                                                        gAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGGCATGTGATACCTGTGGGAGCGGGc                                                               <  3:60381/149‑1 (MQ=255)
                                                                           gagCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                     >  8:274982/1‑92 (MQ=255)
                                                                           gagCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                     <  7:274982/92‑1 (MQ=255)
                                                                           gagCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGc                                                               >  1:53901/1‑146 (MQ=255)
                                                                                                                       ttAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCACCCACAGGAATAACCAGACCAGGAGGGGAAGTGATACCGGTGGGAGCGGGCGGGGCCGCGACGGGGCGCGAGGGGGGCCCAATTGAATTTGCggggtt                >  4:282149/1‑145 (MQ=19)
                                                                                                                           ggCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                         <  7:216367/88‑1 (MQ=255)
                                                                                                                           ggCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                         >  8:216367/1‑88 (MQ=255)
                                                                                                                             ccTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGTGGTCATGTGATACCTGTGGGAGCGGGCCGGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCt                  >  1:34029/1‑141 (MQ=255)
                                                                                                                             ccTGAAATCGGGGCTGCTTTCCAGCCCGTCGCGGGCACGCCCGCTCCCCCAGGAATCACCAGACCTGGAGGTCATGGGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCt                  <  2:34029/141‑1 (MQ=255)
                                                                                                                                   cgaagggCTGCAAAGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgtg                                                                        >  1:332066/5‑81 (MQ=21)
                                                                                                                                      ggggCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCACCGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGa  <  3:282149/148‑1 (MQ=255)
                                                                                                                                      ggggCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                    >  2:110422/1‑66 (MQ=14)
                                                                                                                                      ggggCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                    <  1:110422/66‑1 (MQ=14)
                                                                                                                                                |                                                                                                                                         
GAAAACTTGGCCATCGACGGCGGCGGTCACCTGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGA  >  NC_002947/1013050‑1013331

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: