Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,073,563 A→G 12.1% L433P (CTG→CCG)  gatA ← glutamyl‑tRNA(Gln) amidotransferase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,073,5630AG12.1% 73.1 / 3.6 33L433P (CTG→CCG) gatAglutamyl‑tRNA(Gln) amidotransferase subunit A
Reads supporting (aligned to +/- strand):  ref base A (16/13);  new base G (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-01

AGCCGTTGGGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTCGACATCGTTGAAGGCGGCCATGAAGTCGTTC  >  NC_002947/1073414‑1073712
                                                                                                                                                     |                                                                                                                                                     
aGCCGTTGGGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCa                                                                                                                                                       <  1:790700/150‑1 (MQ=255)
        ggggCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGcc                                                                                                                                               >  1:374988/1‑150 (MQ=255)
             gcgGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGtt                                                                                                                                          >  2:706813/1‑150 (MQ=255)
                   aTGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGGGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGt                                                                                                                                    <  1:502194/150‑1 (MQ=255)
                     gCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGa                                                                                                                                  >  1:10187/1‑150 (MQ=255)
                     gCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGa                                                                                                                                  <  1:586082/150‑1 (MQ=255)
                            ggTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCGGGCCTGCCAggc                                                                                                                                           <  1:116995/134‑2 (MQ=255)
                            ggTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCGGGCCTGCCAggc                                                                                                                                           >  2:116995/1‑133 (MQ=255)
                                        gTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTa                                                                                                                           <  2:62899/138‑1 (MQ=255)
                                        gTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTa                                                                                                                           >  1:62899/1‑138 (MQ=255)
                                              gTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAggc                                                                                                         >  1:197592/1‑150 (MQ=255)
                                              gTGCGCGACGCTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAggc                                                                                                         <  2:508373/150‑1 (MQ=255)
                                                gcgcGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAggcgg                                                                                                       <  1:706813/150‑1 (MQ=255)
                                                       gtggAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCgg                                                                                                >  1:247232/4‑150 (MQ=255)
                                                            gCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGGGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCgc                                                                                           >  2:616582/1‑150 (MQ=255)
                                                               ggcggcCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGGCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCgctgc                                                                                        <  2:247232/150‑1 (MQ=255)
                                                                   gcCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTGGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCtt                                                                                    >  2:157738/1‑150 (MQ=255)
                                                                     ctcccTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTgg                                                                                  <  1:606762/147‑1 (MQ=255)
                                                                     cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATg                                                                                                                                >  1:698545/1‑104 (MQ=255)
                                                                     cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATg                                                                                                                                <  2:698545/104‑1 (MQ=255)
                                                                     cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTgg                                                                                  >  1:441235/1‑150 (MQ=255)
                                                                         cTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCGCTTGgcgtc                                                                              >  1:336663/1‑148 (MQ=255)
                                                                                        cagcagTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTcc                                                                      >  1:325063/1‑143 (MQ=255)
                                                                                        cagcagTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTcc                                                                      <  2:325063/143‑1 (MQ=255)
                                                                                              ttGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGTCGGCACCCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGTCGGGTTTgg                                                         <  2:10187/150‑1 (MQ=255)
                                                                                                            gCAGGCCGTCGACGAAGCTGGCAGGCATCGACAGGGCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGTGGCAATCGGGTCGCTACTCTTGAGGCCGAGCTTCCAGGGCGGGGGTAGCGTGGGCGGGCCgg                                           >  2:155721/1‑149 (MQ=255)
                                                                                                               ggCCGTCTACGAAGCCCGCAGGCATCGACAGGCATGGCAGGGCTGCCTGTTTGCCGGTGATGGTGTAGACGTCTTCCATATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGaa                                        <  1:616582/150‑1 (MQ=255)
                                                                                                                    tCGACGAAGCCCGCAGGAGTAGATAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAgg                                   <  1:683899/150‑1 (MQ=255)
                                                                                                                      gacgaAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTc                                 >  2:233260/1‑150 (MQ=255)
                                                                                                                              ccGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGTGAACCGGGTCGCTGCGCTTGTCGTCAAGCTTCCAGGGCGGGTTGGGCGTGGTCGGGCCGAGAATGAGCTCGACATCGt                         >  1:60417/1‑150 (MQ=255)
                                                                                                                                             aaccTGGCGGGCCTGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGAt                                                                                                             <  2:161547/49‑1 (MQ=25)
                                                                                                                                             aaccTGGCGGGCCTGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGAt                                                                                                             >  1:161547/3‑51 (MQ=25)
                                                                                                                                              gCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGt                                  >  2:452223/1‑125 (MQ=255)
                                                                                                                                              gCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGt                                  <  1:452223/125‑1 (MQ=255)
                                                                                                                                                     aGGCCCGCCAGGTTGGCGGTGACGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGTCCGGGTTTGGCGTGGGCGGGCCGAGAATCAGGTCGACATTGTTGAAGGCGGCCATGAAGTCGTTc  <  1:155721/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCCGTTGGGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTCGACATCGTTGAAGGCGGCCATGAAGTCGTTC  >  NC_002947/1073414‑1073712

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