Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 639,654 T→C 7.1% intergenic (+146/‑99) ahpC → / → ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913639,6540TC 7.1% 177.6 / 3.6 56intergenic (+146/‑99)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base C (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATT  >  NC_000913/639510‑639802
                                                                                                                                                |                                                                                                                                                    
tttCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCgg                                                                                                                                                  >  2:550731/1‑149 (MQ=255)
   ccTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAc                                                                                                                                               <  1:115656/149‑1 (MQ=255)
    cTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAcc                                                                                                                                              >  4:275697/1‑149 (MQ=255)
        gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCACCCGaa                                                                                                                                          <  4:282695/149‑2 (MQ=255)
        gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCCGCAGCACCCGaa                                                                                                                                          <  3:307009/149‑2 (MQ=255)
         tCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                         >  4:106028/1‑149 (MQ=255)
         tCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                         >  1:45745/1‑149 (MQ=255)
           tttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtttt                                                                                                                                       >  2:49420/1‑149 (MQ=255)
               aCGCCATAGCGGCGTTGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                  >  2:243837/1‑149 (MQ=255)
                cGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                                                                                                                  <  2:45173/149‑1 (MQ=255)
                               ggCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCa                                                                                                                   >  3:446008/1‑149 (MQ=255)
                                gCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCaa                                                                                                                  <  3:86538/149‑1 (MQ=255)
                                     gCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                             <  1:23717/149‑1 (MQ=255)
                                           tCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAGCTTGCATcatg                                                                                                       <  1:43266/149‑4 (MQ=255)
                                           tCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg                                                                                                       <  3:275697/149‑1 (MQ=255)
                                           tCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg                                                                                                       <  3:94963/149‑1 (MQ=255)
                                             aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGt                                                                                                     >  4:88645/1‑149 (MQ=255)
                                              ccccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTa                                                                                                    <  4:140469/149‑1 (MQ=255)
                                                        ttacttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAgg                                                                                          >  2:75520/1‑149 (MQ=255)
                                                           cttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCg                                                                                       >  1:127800/1‑149 (MQ=255)
                                                            ttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                             <  1:250443/126‑1 (MQ=255)
                                                            ttGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                             >  2:250443/1‑126 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                         >  4:503980/1‑95 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                                                                                                                         <  3:503980/95‑1 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                             <  2:31096/123‑1 (MQ=255)
                                                               tgtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGc                                                                                                             >  1:31096/1‑123 (MQ=255)
                                                                gtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGc                                                                                  >  1:181163/1‑149 (MQ=255)
                                                                 tAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCa                                                                                 <  2:127800/149‑1 (MQ=255)
                                                                    gCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatga                                                                              >  3:187489/1‑149 (MQ=255)
                                                                     cTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgat                                                                             <  4:187489/149‑1 (MQ=255)
                                                                        ccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgatgat                                                                          >  4:260981/1‑149 (MQ=255)
                                                                          ggggATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                        <  3:415690/149‑1 (MQ=255)
                                                                                tCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAg                                                                  <  2:263375/149‑1 (MQ=255)
                                                                                 cacaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGc                                                                 >  2:97845/1‑149 (MQ=255)
                                                                                   caGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGccc                                                               >  4:202574/1‑149 (MQ=255)
                                                                                    aGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCa                                                              <  3:88645/149‑1 (MQ=255)
                                                                                    aGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCa                                                              <  2:181163/149‑1 (MQ=255)
                                                                                     ggCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAg                                                             >  2:131740/1‑149 (MQ=255)
                                                                                             gccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaaa                                                     >  2:432370/1‑149 (MQ=255)
                                                                                              ccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAAc                                                    <  2:45745/149‑1 (MQ=255)
                                                                                                    tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttc                                                                                                                                      >  3:525360/1‑61 (MQ=255)
                                                                                                    tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttc                                                                                                                                      <  4:525360/61‑1 (MQ=255)
                                                                                                          gcgAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGAcaca                                        >  2:145978/1‑149 (MQ=255)
                                                                                                                tACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATg                                  >  3:449549/1‑149 (MQ=255)
                                                                                                                    tCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAata                                    <  4:104913/143‑1 (MQ=255)
                                                                                                                    tCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAata                                    >  3:104913/1‑143 (MQ=255)
                                                                                                                                cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                  <  1:243837/149‑1 (MQ=255)
                                                                                                                                cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                  >  4:472936/1‑149 (MQ=255)
                                                                                                                                 tAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                         >  2:390497/1‑141 (MQ=255)
                                                                                                                                 tAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                         <  1:390497/141‑1 (MQ=255)
                                                                                                                                  aaTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCt                >  3:141692/1‑149 (MQ=255)
                                                                                                                                   aTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCtt               >  1:570871/1‑149 (MQ=255)
                                                                                                                                    ttcttcGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTa              >  3:326962/1‑149 (MQ=255)
                                                                                                                                      aatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGAAGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTAcc            >  4:332048/3‑149 (MQ=255)
                                                                                                                                        tcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCtt          >  2:390876/1‑149 (MQ=255)
                                                                                                                                          gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTga        >  1:221606/1‑149 (MQ=255)
                                                                                                                                                tgcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAAtt  >  1:258866/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATT  >  NC_000913/639510‑639802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: