Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,852,280 C→T H442Y (CAC→TAC)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,852,2800CT97.3% 116.9 / ‑6.2 37H442Y (CAC→TAC) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/16);  minor base G (0/1);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGCACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCACCAGCAAGCACGGCAACCCCAGCCGTGACTGGTTGAA  >  NC_002947/1852137‑1852416
                                                                                                                                               |                                                                                                                                        
aCGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCTGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCg                                                                                                                                    >  1:695095/1‑150 (MQ=255)
 cGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCga                                                                                                                                   <  1:688088/150‑1 (MQ=255)
                 aTCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGc                                                                                                                   <  1:555064/150‑1 (MQ=255)
                 aTCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGc                                                                                                                   >  1:987484/1‑150 (MQ=255)
                       ggcctggccGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCgcggc                                                                                                             <  1:142584/150‑1 (MQ=255)
                             gccGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAg                                                                                                       >  1:127253/1‑150 (MQ=255)
                                      ttGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCAACCGCGTGTACACAGATTTCTGGCCCAACTGCCGCGGCGCGACGATCAACGGc                                                                                              >  2:379753/1‑150 (MQ=255)
                                               gACAGCGAGCCTGACCGGGTTTATGTGTTCACCCCCCACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTTGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTg                                                                                     <  2:1008774/150‑1 (MQ=255)
                                                aCATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGCCCGTATCGTGc                                                                                    >  1:237154/1‑150 (MQ=255)
                                                   tCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGc                                                                                 >  1:673457/1‑150 (MQ=255)
                                                    cGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCt                                                                                <  2:127253/150‑1 (MQ=255)
                                                       gCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCCCGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCCAAGATCAACGGACGTCTCGAGCCCCCGaa                                                                             >  2:189583/1‑150 (MQ=255)
                                                        ccTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAAc                                                                            <  2:987484/150‑1 (MQ=255)
                                                        ccTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAAc                                                                            >  1:1028101/1‑150 (MQ=255)
                                                            aCCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACa                                                                        <  2:927628/150‑1 (MQ=255)
                                                              cGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGc                                                                      >  2:1034019/1‑150 (MQ=255)
                                                              cGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGc                                                                      >  1:224964/1‑150 (MQ=255)
                                                                ggTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACAcc                                                                                                                                     <  1:559059/85‑1 (MQ=255)
                                                                ggTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACAcc                                                                                                                                     >  2:559059/1‑85 (MQ=255)
                                                                   ttATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGct                                                                 <  2:673457/150‑2 (MQ=255)
                                                                         gtTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTgg                                                           >  1:950920/1‑150 (MQ=255)
                                                                         gtTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTgg                                                           >  2:366893/1‑150 (MQ=255)
                                                                            cACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGa                                                        >  1:607365/1‑150 (MQ=255)
                                                                             accaccGACGGCTACGCGATCGACCTGCCCAAAGGCGCCACGCAATCGGCTCTCGCCTACCGCGTGTCCACCCAGATCGGCCACAACGGCCGCGGCGTGAAGATCCACCGCCGTCTCCTGCCGCCGCACTACAGCCTGCAGACTGGCGAg                                                       <  2:66849/150‑1 (MQ=255)
                                                                             aCCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGgcgc                                                                                                           <  2:16596/98‑1 (MQ=255)
                                                                             aCCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGgcgc                                                                                                           >  1:16596/1‑98 (MQ=255)
                                                                             aCCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAg                                                       <  1:1034019/150‑1 (MQ=255)
                                                                                   gACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTg                                                 <  2:826633/150‑1 (MQ=255)
                                                                                              cATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCAc                                      >  2:842314/1‑150 (MQ=255)
                                                                                               aTCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGCCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCAcc                                     <  2:237154/150‑1 (MQ=255)
                                                                                                  gACCTGCCCAAAGGGGCCACGCTATTGGGGTTGACCGACGTCGTGTACACTGAGAGCGGCCACAACTGCCGCGGCGCGAAGATCAACGTCCGTATCGTGCCGCTGAAATACACCCCACCGACTTGCGAGGATGTATAGATTATCACCagc                                  <  2:358800/150‑1 (MQ=255)
                                                                                                    ccTGCCCAAAGGCGCCCCGCAACTGGGTTTAGCCTACGGGGGGGAGACCGTGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCACCagcaa                                <  2:695095/150‑1 (MQ=255)
                                                                                                        cccAAAGGGGCCACACCATTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGACGGGAGACGACCAACGCGCGTATCGAGCCACCGACCGACGCGCTGCAGACAGGCCAAAAATTGGCGACCATCACCagcaagcaa                            >  2:122723/1‑149 (MQ=255)
                                                                                                                  ccacgccaTTGGACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCACCAGCAAGCACGGCAACCCCa                  >  1:642181/1‑150 (MQ=255)
                                                                                                                        caTTGGACTTCGCCTACCGCAAGTACACCGAGATCGGCCACAACTGCCGCTGGAAAAAGATCAACAGCCATATCGTGACGCTCAACTACAGACTGCAGACTACCGAGCACGTGGAGATCATCACCAGCAAGCACGGCAACCCCAGCCGTg            >  2:421046/1‑150 (MQ=255)
                                                                                                                             gACTTCGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCACCAGCAAGCACGGCAACCCCAGCCGTGACTgg       >  1:1102626/1‑150 (MQ=255)
                                                                                                                                  cGCCTACCGCGTGTACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCACCAGCAAGCACGGCAACCCCAGCCGTGACTGGTTGaa  <  1:366893/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
ACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTGACCGGGTTTATGTGTTCACCCCCGACGGCCACGCCATCGACCTGCCCAAAGGCGCCACGCCATTGGACTTCGCCTACCGCGTGCACACCGAGATCGGCCACAACTGCCGCGGCGCGAAGATCAACGGCCGTATCGTGCCGCTGAACTACAGCCTGCAGACTGGCGAGCAGGTGGAGATCATCACCAGCAAGCACGGCAACCCCAGCCGTGACTGGTTGAA  >  NC_002947/1852137‑1852416

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: