Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 371,542 C→T Q102* (CAA→TAA)  mhpC → 2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,5420CT96.7% 87.0 / ‑6.1 30Q102* (CAA→TAA) mhpC2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (16/13);  total (17/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 = 1.00e+00

CAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGT  >  NC_000913/371406‑371673
                                                                                                                                        |                                                                                                                                   
cAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTaa                                                                                                                                   >  2:614632/1‑139 (MQ=255)
             cTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCa                                                                                                                      >  1:615501/1‑139 (MQ=255)
                  agagGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATc                                                                                                                 <  2:841985/139‑1 (MQ=255)
                              tCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGc                                                                                                     >  1:773989/1‑139 (MQ=255)
                                     aTCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGa                                                                                              <  1:614632/139‑1 (MQ=255)
                                        ctgctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCACAATCCACCTGCTGGGCAACTCGATgg                                                                                           >  1:381464/1‑139 (MQ=255)
                                               aTTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCa                                                                                    <  2:615501/139‑1 (MQ=255)
                                               aTTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCa                                                                                    <  2:754433/139‑1 (MQ=255)
                                                ttGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCAt                                                                                   >  1:148054/1‑139 (MQ=255)
                                                 tGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTgg                                                                                                                               >  2:76586/1‑94 (MQ=255)
                                                 tGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTgg                                                                                                                               <  1:76586/94‑1 (MQ=255)
                                                         ttGGGGCAAGAGCGATTCGGCCGGTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCtgtg                                                                          >  1:56224/1‑139 (MQ=255)
                                                              gCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGtt                                                                     <  2:381464/139‑1 (MQ=255)
                                                                 agagCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGtt                                                                     <  1:59343/136‑1 (MQ=255)
                                                                 agagCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGtt                                                                     >  2:59343/1‑136 (MQ=255)
                                                                  gagCGATTCGGTCGTTAATAGTGGTTAGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCGAAAATCAACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGAGTGCAcc                                                                 >  2:783416/1‑139 (MQ=255)
                                                                      gATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTa                                                             <  2:482018/139‑1 (MQ=255)
                                                                                         gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCgg                                          <  2:773989/139‑1 (MQ=255)
                                                                                            cgcgATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaa                                       >  1:227661/1‑139 (MQ=255)
                                                                                                               aCGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCAcc                                                                                                             <  2:307940/50‑1 (MQ=255)
                                                                                                               aCGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCAcc                                                                                                             >  1:307940/1‑50 (MQ=255)
                                                                                                               aCGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                    >  1:412925/1‑139 (MQ=255)
                                                                                                                cGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCaa                                                                                                   <  2:633035/59‑1 (MQ=255)
                                                                                                                cGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCaa                                                                                                   >  1:633035/1‑59 (MQ=255)
                                                                                                                      cTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACgg             >  2:433991/1‑139 (MQ=255)
                                                                                                                           aaGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcggc        <  1:433991/139‑1 (MQ=255)
                                                                                                                              cGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCa                                                                                    <  2:596141/60‑1 (MQ=255)
                                                                                                                              cGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCa                                                                                    >  1:596141/1‑60 (MQ=255)
                                                                                                                              cGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATg     >  1:749010/1‑139 (MQ=255)
                                                                                                                                 ggtggATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgt  >  2:461901/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGT  >  NC_000913/371406‑371673

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