Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,231 A→T 18.0% intergenic (+56/+75) 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,2310AT18.0% 94.3 / 13.7 39intergenic (+56/+75)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (17/15);  new base T (4/3);  total (21/18)
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.56e-01

TGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCA  >  NC_002947/1013084‑1013376
                                                                                                                                                   |                                                                                                                                                 
tGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTAt                                                                                                                                                 >  2:188470/1‑150 (MQ=255)
 gACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATc                                                                                                                                                >  1:583808/1‑150 (MQ=255)
     tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                            >  1:791517/1‑150 (MQ=255)
            cgcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATTAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGTGCTGCTTTGCAGCCCTTCGCGTGCACGCCCGCTCCCACAGAAATCACATGACCtcc                                                                                                                                     <  2:76907/150‑3 (MQ=39)
             gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                    <  2:583808/150‑1 (MQ=255)
                   ctgtatgcgctggaggcgatgccggcgaagatcgaccgtattcgtacgcgaatctggaccagtatcaaagccggaaattgaaatggggccgcctagcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACCAGACCTGGAGGTCAt                                                                                                                              >  2:692947‑M1/99‑150 (MQ=255)
                   ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCATAGGAATCACCAGACCTGGAGGTTAt                                                                                                                              <  1:696651/150‑1 (MQ=255)
                   ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                                                              >  2:346136/1‑150 (MQ=255)
                   ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                                                              <  1:288157/150‑1 (MQ=255)
                    cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGAAGCTTTGTTGTCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                             <  2:308355/150‑1 (MQ=255)
                            ctggaggcgatgccggcgaagatcgaccgtattcgtacgcgaatctggaccagtatcaaagccggaaattgaactgcggccgcctagcgGCCCGTCTCGTGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                     >  1:495945‑M1/90‑150 (MQ=255)
                                gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                 <  2:517733/150‑1 (MQ=21)
                                 cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                            >  1:368665/1‑122 (MQ=255)
                                 cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                            <  2:368665/122‑1 (MQ=255)
                                    tGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                             <  1:108810/134‑1 (MQ=255)
                                    tGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                             >  2:108810/1‑134 (MQ=255)
                                                  tcgaccgtattcgtacgcgaatctggaccagtatcaaagccggaaattgaaatggggccgcctagcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                     >  2:563874‑M1/68‑128 (MQ=255)
                                                  tcgaccgtattcgtacgcgaatctggaccagtatcaaagccggaaattgaaatggggccgcctagcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATCTGATACCt                                                                                                                     <  1:563874‑M1/61‑1 (MQ=255)
                                                   ggccggccATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCATACCTGGAGGTCATGTGATACCTATAGGAGCGTGCGTGCCCTCGaa                                                                                              >  2:61183/1‑150 (MQ=38)
                                                    gaccgtattcgtacgcgaatctggaccagtatcaaagccggaaattgaaatggggccgcctagcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                   >  1:363169‑M1/66‑128 (MQ=255)
                                                      cGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                           <  1:408823/150‑1 (MQ=255)
                                                              aaGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGACCGCTGTGGGAGTGGGCGTGCCCGCGAAGGGctgcaaag                                                                                   <  1:188470/150‑8 (MQ=14)
                                                                         tACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCAc                                                                                                                                              <  2:87459/80‑1 (MQ=255)
                                                                         tACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCAc                                                                                                                                              >  1:87459/1‑80 (MQ=255)
                                                                                     ccATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGATCTGTAGGGCATCTGATACCTGTGGGAGCGGGCGGGCCCGCGAAGTGCCGCTAGGCGGCTCCATTTCAATTTCCGGCt                                                            >  1:774791/1‑150 (MQ=38)
                                                                                      cATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCtt                                                           >  1:165406/1‑150 (MQ=255)
                                                                                      aagccggaaattgaaatggggccgcctagcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGa                                                                                                                          >  2:574061‑M1/32‑87 (MQ=255)
                                                                                      aagccggaaattgaaatggggccgcctagcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGa                                                                                                                          <  1:574061‑M1/56‑1 (MQ=255)
                                                                                            ggCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATAc                                                     <  2:165406/150‑1 (MQ=255)
                                                                                                          gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                       <  1:346136/150‑1 (MQ=255)
                                                                                                                  gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAAt                               >  1:627907/1‑150 (MQ=255)
                                                                                                                    aGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAgg                                                                                                                                  >  2:55785/1‑49 (MQ=255)
                                                                                                                    aGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAgg                                                                                                                                  <  1:55785/49‑1 (MQ=255)
                                                                                                                        CTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc                         >  2:627907‑M1/1‑150 (MQ=255)
                                                                                                                        CTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc                         >  2:791517‑M1/1‑150 (MQ=255)
                                                                                                                                            CCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAggc                                                                                  >  1:576024‑M1/1‑73 (MQ=255)
                                                                                                                                            CCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAggc                                                                                  <  2:576024‑M1/73‑1 (MQ=255)
                                                                                                                                              ACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc                         >  2:363169‑M1/1‑128 (MQ=255)
                                                                                                                                               CAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa  <  1:421241‑M1/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGACCTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCA  >  NC_002947/1013084‑1013376

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