Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,788,115 C→T 14.8% P132L (CCG→CTG)  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,1150CT14.8% 60.2 / 4.3 27P132L (CCG→CTG) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base T (2/2);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

CACGAGTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTA  >  NC_002947/2787970‑2788241
                                                                                                                                                 |                                                                                                                              
cACGAGTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGAc                                                                                                                            <  2:166160/150‑1 (MQ=255)
     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)
                                                                                                                                                 |                                                                                                                              
CACGAGTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTA  >  NC_002947/2787970‑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: