Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,788,121 T→A 16.0% V134D (GTC→GAC)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,788,1210TA16.0% 58.6 / 5.7 25V134D (GTC→GAC) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base A (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

GTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTA  >  NC_002947/2787975‑2788241
                                                                                                                                                  |                                                                                                                        
gTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAATGCAGTTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGAcgtcgt                                                                                                                       >  1:7727/1‑150 (MQ=255)
          ttGGTGCGGGCACTGCTGCAAGTGGGTGGCCACCCGGCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTg                                                                                                             <  2:41293/150‑1 (MQ=255)
                 tGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGa                                                                                                      <  2:562245/150‑1 (MQ=255)
                         ctgcAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACt                                                                                                              >  2:507413/1‑134 (MQ=255)
                         ctgcAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACt                                                                                                              <  1:507413/134‑1 (MQ=255)
                               gtgggtggCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAAc                                                                                        <  1:722989/150‑1 (MQ=255)
                                       ccACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagttc                                                                                                            <  2:85889/122‑5 (MQ=255)
                                       ccACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagttc                                                                                                            >  1:85889/1‑118 (MQ=255)
                                        cACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagt                                                                                                              <  1:192578/119‑3 (MQ=255)
                                        cACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagt                                                                                                              >  2:192578/1‑117 (MQ=255)
                                         aCCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgc                                                                              >  1:120978/1‑150 (MQ=255)
                                                             gAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGt                                                          >  1:715647/1‑150 (MQ=255)
                                                                       aGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCACGGAAGTCGTACAGGCACTACACCTGATGGCGGTCCTCCTCCCCAACTGGCTCCATATGGCCATCGACGGTGCGCTGTTcc                                                <  2:7727/150‑1 (MQ=255)
                                                                        gCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCa                                               <  2:715647/150‑1 (MQ=255)
                                                                           gATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGa                                            <  2:120978/150‑1 (MQ=255)
                                                                                aGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAg                                       >  2:424304/1‑150 (MQ=255)
                                                                                       ggAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATc                                <  1:412648/150‑1 (MQ=255)
                                                                                               aaTTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGa                        >  1:82348/1‑150 (MQ=255)
                                                                                                 ttCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATc                      >  2:676486/1‑150 (MQ=255)
                                                                                                  tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTAc                                                                                             >  1:670648/1‑78 (MQ=255)
                                                                                                  tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTAc                                                                                             <  2:670648/78‑1 (MQ=255)
                                                                                                      cTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGc                 >  1:37998/1‑150 (MQ=255)
                                                                                                             aCCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCg          >  1:208183/1‑150 (MQ=255)
                                                                                                                     ctcGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTa  <  2:208183/150‑1 (MQ=255)
                                                                                                                                     aaCTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTg                                     >  2:116613/1‑99 (MQ=255)
                                                                                                                                     aaCTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTg                                     <  1:116613/99‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                        
GTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTA  >  NC_002947/2787975‑2788241

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