Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,237 C→A 8.7% intergenic (+62/+69) 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,2370CA 8.7% 137.5 / 2.4 46intergenic (+62/+69)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base A (2/2);  total (24/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 = 9.36e-01

CTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAG  >  NC_002947/1013092‑1013386
                                                                                                                                                 |                                                                                                                                                       
ctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGAc                                                                                                                                                     <  1:165235/150‑1 (MQ=255)
        gACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGt                                                                                                                                             >  1:195451/1‑150 (MQ=255)
               gTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAGGTGa                                                                                                                                      <  2:1139940/150‑1 (MQ=255)
                  aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTT‑CGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAcc                                                                                                                                  >  1:282404/1‑150 (MQ=255)
                    ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                                 <  2:569204/150‑1 (MQ=255)
                      tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAATTGCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGACCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  1:162447/150‑1 (MQ=255)
                           cTGATGGAAGAGCGCTGGCCGGACATCAAGTT‑CGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGc                                                                                                                         <  2:282404/150‑1 (MQ=255)
                            tGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               >  1:81722/1‑144 (MQ=21)
                            tGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               <  2:81722/144‑1 (MQ=21)
                             gATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                >  2:183564/1‑126 (MQ=255)
                             gATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                <  1:183564/126‑1 (MQ=255)
                                        gttGGCCGGACATCAAGTTTCGCGCTACCCGCGAGCACAGTTAAGGACTGAAATCGGGGCTGCGTTGCAGCCCTCCTCGGGCACGCCCACCCCCTCAGGAATCCCCCAGCCTGCCGGTCATGTGATACCTGTGGGAGCCGGCGTGC‑Ccgc                                                                                                            <  2:365955/148‑1 (MQ=255)
                                          tGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGt                                                                                                                   >  1:222769/1‑142 (MQ=255)
                                          tGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGt                                                                                                                   <  2:222769/142‑1 (MQ=255)
                                                        gTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                               >  2:1127586/1‑100 (MQ=255)
                                                        gTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                               <  1:1127586/100‑1 (MQ=255)
                                                         ttCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAggc‑gg                                                                                          >  2:83790/1‑150 (MQ=255)
                                                                    ccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCCCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCAGGCGTGC‑CCGCGAAGAGCCGCTAGGC‑GGCCCCATTTCaa                                                                               >  1:478411/1‑150 (MQ=255)
                                                                            aCCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATa                                                                                                                                    >  1:515389/1‑91 (MQ=255)
                                                                            aCCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATa                                                                                                                                    <  2:515389/91‑1 (MQ=255)
                                                                                ttAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAg                                                                                                        <  2:191832/114‑1 (MQ=255)
                                                                                ttAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAg                                                                                                        >  1:191832/1‑114 (MQ=255)
                                                                                  aaGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGAc                                                                                                                                                     <  1:92716/68‑1 (MQ=255)
                                                                                  aaGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGAc                                                                                                                                                     >  2:92716/1‑68 (MQ=255)
                                                                                     gCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑GCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACt                                                              >  1:684666/1‑150 (MQ=255)
                                                                                      ccTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTg                                                             >  2:589993/1‑150 (MQ=255)
                                                                                      ccTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCATCACGCCGGCATCACCAGTCCTGGCGGTGATGTGATACCTGTGGGGGCTGGGGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTg                                                             <  2:739425/150‑1 (MQ=255)
                                                                                        tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                >  2:181926/1‑67 (MQ=255)
                                                                                        tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                <  1:181926/67‑1 (MQ=255)
                                                                                             ggggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAGCAGACCTGGAGGTGATGGGATACCTGTGGGTGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGa                                                      <  2:478411/148‑1 (MQ=255)
                                                                                                                 ttCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACCGTCg                                  >  2:575096/1‑150 (MQ=255)
                                                                                                                          aCGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGcc                         >  2:248032/1‑150 (MQ=255)
                                                                                                                            gCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCCCGTAAGAATACGGTAGATCTTCTCCgg                       >  1:525837/1‑150 (MQ=255)
                                                                                                                                gCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATc                   <  1:589993/150‑1 (MQ=255)
                                                                                                                                  tCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTg                                                             <  1:879295/106‑1 (MQ=11)
                                                                                                                                  tCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTg                                                             >  2:879295/1‑106 (MQ=11)
                                                                                                                                      acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTc                                                          >  1:870188/1‑105 (MQ=255)
                                                                                                                                      acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTc                                                          <  2:870188/105‑1 (MQ=255)
                                                                                                                                      acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTAc                                             >  2:330872/1‑118 (MQ=255)
                                                                                                                                      acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTAc                                             <  1:330872/118‑1 (MQ=255)
                                                                                                                                      acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGAGCAGGTCGCGGTCCACCACG‑TGGGAGGCTTGGC‑CGGCCCATTTCAATTTCCGGCTTTGATACGGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc             >  1:328108/1‑150 (MQ=255)
                                                                                                                                       caGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa            >  2:1044976/1‑150 (MQ=255)
                                                                                                                                          gAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCtt                                                                     <  2:1180301/90‑1 (MQ=255)
                                                                                                                                          gAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCtt                                                                     >  1:1180301/1‑90 (MQ=255)
                                                                                                                                          gAATCACCAGACCTGGAGGTCATGTGATAACTGTGGGGGGGGGCGTGC‑CTTCGAAGTTGCACCAGGCGGGCGCCATTCTACGTTTCGGGTAT‑TGAGTGGTCCTGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcg         >  1:274919/1‑150 (MQ=255)
                                                                                                                                               aCCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATAc    <  1:83790/150‑1 (MQ=255)
                                                                                                                                                 cAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAg  <  2:274919/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGC‑CCGCGAAGGGCCGCTAGGC‑GGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAG  >  NC_002947/1013092‑1013386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑