Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,011,483 T→A 14.3% intergenic (‑418/‑109) PP_0871 ← / → prfC glycine betaine/carnitine/choline ABC transporter permease/peptide chain release factor 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,011,4830TA14.3% 64.5 / 5.1 28intergenic (‑418/‑109)PP_0871/prfCglycine betaine/carnitine/choline ABC transporter permease/peptide chain release factor 3
Reads supporting (aligned to +/- strand):  ref base T (14/10);  new base A (2/2);  total (16/12)
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.58e-01

CTGCGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGCGCC  >  NC_002947/1011344‑1011625
                                                                                                                                           |                                                                                                                                              
cTGCGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAACCTGACTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGCACCTGTAGGCGCGGGTTTACCCGCGAAg                                                                                                                                      >  1:138283/1‑150 (MQ=255)
                  tagTGATAAAGAGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGccc                                                                                                                    >  2:692387/3‑149 (MQ=255)
                             cGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCa                                                                                                         <  2:240859/150‑1 (MQ=255)
                                    cTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTg                                                                                                  <  2:275422/150‑1 (MQ=255)
                                      gTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGt                                                                                                <  2:138283/150‑1 (MQ=255)
                                           gtttgtttTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGAtttt                                                                                                   <  2:95415/142‑1 (MQ=255)
                                           gtttgtttTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGAtttt                                                                                                   >  1:95415/1‑142 (MQ=255)
                                           gtttgtttTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAgg                                                                                           >  1:22345/1‑150 (MQ=255)
                                                ttttgtttGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGCACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGCCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCaa                                                                                      >  2:549960/1‑150 (MQ=255)
                                                ttttgtttGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCaa                                                                                      >  1:580241/1‑150 (MQ=255)
                                                ttttgtttGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                        <  1:369843/100‑2 (MQ=255)
                                                ttttgtttGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                        >  2:369843/1‑99 (MQ=255)
                                                      ttgtGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAg                                                                                <  2:580241/150‑1 (MQ=255)
                                                                ctctTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTAccc                                                                                                                                            >  1:449065/1‑80 (MQ=255)
                                                                ctctTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTAccc                                                                                                                                            <  2:449065/80‑1 (MQ=255)
                                                                ctctTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACt                                                                      >  1:327709/1‑150 (MQ=255)
                                                                 tctTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTg                                                                     >  2:253396/1‑150 (MQ=255)
                                                                                   cTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCa                                                   <  2:22345/150‑1 (MQ=255)
                                                                                          gAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGc                                            >  1:157183/1‑150 (MQ=255)
                                                                                                   aCCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCa                                               <  1:688681/138‑1 (MQ=255)
                                                                                                   aCCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCa                                               >  2:688681/1‑138 (MQ=255)
                                                                                                                        gTACCTGTAGGGGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAg              >  2:125871/1‑150 (MQ=255)
                                                                                                                        gTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATa                                                                        <  2:391729/92‑1 (MQ=255)
                                                                                                                        gTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATa                                                                        >  1:391729/1‑92 (MQ=255)
                                                                                                                            cTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAGGTcgc          <  1:253396/150‑1 (MQ=255)
                                                                                                                                cttcgcgGGTAAACCCGCGAAGGGGCCGGCCCAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAg      >  1:421183/4‑150 (MQ=255)
                                                                                                                                    gcgGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGcgcc  >  1:395247/1‑150 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGcgcc  <  2:421183/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CTGCGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGCGCC  >  NC_002947/1011344‑1011625

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: