Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,073,568 C→T 14.2% A431A (GCG→GCA 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,5680CT14.2% 71.4 / 4.6 34A431A (GCG→GCAgatAglutamyl‑tRNA(Gln) amidotransferase subunit A
Reads supporting (aligned to +/- strand):  ref base C (16/13);  new base T (2/3);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

GGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTCGACATCGTTGAAGGCGGCCATGAAGTCGTTC  >  NC_002947/1073422‑1073712
                                                                                                                                                  |                                                                                                                                                
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                                                                                                                                <  2:698545/104‑1 (MQ=255)
                                                             cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATg                                                                                                                                >  1:698545/1‑104 (MQ=255)
                                                             cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTgg                                                                                  >  1:441235/1‑150 (MQ=255)
                                                                 cTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCGCTTGgcgtc                                                                              >  1:336663/1‑148 (MQ=255)
                                                                                cagcagTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTcc                                                                      <  2:325063/143‑1 (MQ=255)
                                                                                cagcagTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTcc                                                                      >  1:325063/1‑143 (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)
                                                                                                                                                  |                                                                                                                                                
GGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTCGACATCGTTGAAGGCGGCCATGAAGTCGTTC  >  NC_002947/1073422‑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: