Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 371,588 (G)6→7 coding (350/867 nt) 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,5821.G90.9% 65.9 / 0.3 22coding (344/867 nt)mhpC2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (9/11);  total (10/12)
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.91e-01

TCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAAT  >  NC_000913/371444‑371715
                                                                                                                                           |                                                                                                                                     
tCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAt                                                                                                                                        >  2:105097/1‑139 (MQ=255)
     ctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                                                                                  <  2:181618/139‑1 (MQ=255)
     ctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                                                                                  <  1:105097/139‑1 (MQ=255)
      tgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATgggggg                                                                                                                                 >  2:401833/1‑139 (MQ=255)
          ttGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATgggggggcca                                                                                                                             <  2:193510/139‑5 (MQ=255)
                       ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATgg                                                                                                                                     <  2:386681/118‑1 (MQ=255)
                       ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATgg                                                                                                                                     >  1:386681/1‑118 (MQ=255)
                         cAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGtt                                                                                                               <  2:51613/139‑1 (MQ=255)
                         cAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGtt                                                                                                               <  1:33582/139‑1 (MQ=255)
                         cAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGtt                                                                                                               <  1:401833/139‑1 (MQ=255)
                                 aTTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTa                                                                                                       >  2:245433/1‑139 (MQ=255)
                                                 tGGTTCGCGATCGGATCTTACTGCACGAATCCTGAAAA‑‑GTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCATCCTTAAATGGCCGTAGCGCGtgg                                                                                     >  2:887569/1‑137 (MQ=255)
                                                 tGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATgggggggcca                                                                                                                             >  2:952942/1‑96 (MQ=255)
                                                 tGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATgggggggcca                                                                                                                             <  1:952942/100‑5 (MQ=255)
                                                   gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCg                                                                                     >  2:1390534/1‑139 (MQ=255)
                                                    ttCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                                                                                  >  2:346347/1‑92 (MQ=255)
                                                    ttCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                                                                                  <  1:346347/92‑1 (MQ=255)
                                                     tCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGc                                                                                   >  2:667057/1‑139 (MQ=255)
                                                                        cACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcg                                                                <  2:464792/139‑1 (MQ=255)
                                                                         aCGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                                                                                  >  2:924271/1‑71 (MQ=255)
                                                                         aCGAATCCTGAAAAGCGTGCTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                                                                                  <  1:924271/71‑1 (MQ=255)
                                                                          cGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                                                              >  1:472553/1‑139 (MQ=255)
                                                                            aaTCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                    >  1:895478/1‑66 (MQ=255)
                                                                            aaTCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                    <  2:895478/66‑1 (MQ=255)
                                                                                                          tatCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGAt                              >  2:572778/1‑139 (MQ=255)
                                                                                                           atCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATg                             <  1:245433/139‑1 (MQ=255)
                                                                                                           atCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATg                             >  2:1357156/1‑139 (MQ=255)
                                                                                                            tCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGc                            <  2:82159/139‑1 (MQ=255)
                                                                                                            tCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGc                            <  2:472553/139‑1 (MQ=255)
                                                                                                             cGCCAAAATCCACCTGCTGGGCAACTCGATgggggggcca                                                                                                                             <  1:352012/40‑5 (MQ=255)
                                                                                                             cGCCAAAATCCACCTGCTGGGCAACTCGATgggggggcca                                                                                                                             >  2:352012/1‑36 (MQ=255)
                                                                                                             cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCa                          >  1:466777/1‑139 (MQ=255)
                                                                                                               ccAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCaa                         >  1:1476490/1‑139 (MQ=255)
                                                                                                                  aaaTCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      <  1:572778/139‑1 (MQ=255)
                                                                                                                            gctgGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTaa            <  1:1357156/139‑1 (MQ=255)
                                                                                                                             ctgGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGc          <  2:466777/139‑1 (MQ=255)
                                                                                                                                    aCTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGa    <  1:1390534/139‑1 (MQ=255)
                                                                                                                                     cTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGaa   >  2:1423499/1‑139 (MQ=255)
                                                                                                                                      tCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAAt  <  1:1234226/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
TCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAAT  >  NC_000913/371444‑371715

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