Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,422,793 C→T A203V (GCC→GTC)  ulaD → 3‑keto‑L‑gulonate 6‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,7930CT100.0% 158.0 / NA 49A203V (GCC→GTC) ulaD3‑keto‑L‑gulonate 6‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/24);  total (25/24)

TCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCC  >  NC_000913/4422688‑4422918
                                                                                                         |                                                                                                                             
tcaccgGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACt                                                                                              >  2:480094/1‑139 (MQ=255)
      gAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTgggg                                                                                        <  1:752675/139‑1 (MQ=255)
         gCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTa                                                                                     <  1:480094/139‑1 (MQ=255)
           cTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAg                                                                                   >  2:533927/1‑139 (MQ=255)
           cTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAg                                                                                   <  1:533927/139‑1 (MQ=255)
                gcTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCg                                                                              <  2:745061/139‑1 (MQ=255)
                           tGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTc                                                                   <  1:1215991/139‑1 (MQ=255)
                            gCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCg                                                                                                     <  2:100566/104‑1 (MQ=255)
                            gCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCg                                                                                                     >  1:100566/1‑104 (MQ=255)
                                   ttCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                          <  2:815953/92‑1 (MQ=255)
                                   ttCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                          >  1:815953/1‑92 (MQ=255)
                                             ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCAc                                                                                                                          <  1:107768/66‑1 (MQ=255)
                                             ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCAc                                                                                                                          >  2:107768/1‑66 (MQ=255)
                                             ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTTGGGCTAAGGAgggg                                                                             >  2:1213367/1‑111 (MQ=255)
                                             ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCgg                                                                             <  1:1213367/111‑1 (MQ=255)
                                             ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGatat                                                                         >  2:50827/1‑115 (MQ=255)
                                             ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGatat                                                                         <  1:50827/115‑1 (MQ=255)
                                                       cGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGaa                                       >  2:447935/1‑139 (MQ=255)
                                                                cgcgGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCAc                                                                                                                          <  1:1388127/47‑1 (MQ=255)
                                                                cgcgGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCAc                                                                                                                          >  2:1388127/1‑47 (MQ=255)
                                                                    ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCt                                                                                                   >  1:503573/1‑66 (MQ=255)
                                                                    ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCt                                                                                                   <  2:503573/66‑1 (MQ=255)
                                                                    ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                             >  1:1109109/1‑72 (MQ=255)
                                                                    ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                             <  2:1109109/72‑1 (MQ=255)
                                                                    ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGc                                                <  1:837301/117‑1 (MQ=255)
                                                                    ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGc                                                >  2:837301/1‑117 (MQ=255)
                                                                      ccGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccgcc                        >  2:153474/1‑139 (MQ=255)
                                                                         tagtaTCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                     <  1:975494/139‑1 (MQ=255)
                                                                               ccGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgc               >  2:147042/1‑139 (MQ=255)
                                                                               ccGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgc               <  1:447935/139‑1 (MQ=255)
                                                                                cGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGc                                                                                                    <  2:1079742/53‑1 (MQ=255)
                                                                                cGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGc                                                                                                    >  1:1079742/1‑53 (MQ=255)
                                                                                   gATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAgcg                                 >  2:560305/1‑117 (MQ=255)
                                                                                   gATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAgcg                                 <  1:560305/117‑1 (MQ=255)
                                                                                   gATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctggc           >  1:590821/1‑139 (MQ=255)
                                                                                    aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATAtgttgt                                                                    <  1:1576762/81‑1 (MQ=255)
                                                                                    aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATAtgttgt                                                                    >  2:1576762/1‑81 (MQ=255)
                                                                                    aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAAC‑‑TGGGGCTAAGGAGCGGATATGTTGTc                                                                   >  2:895650/1‑80 (MQ=255)
                                                                                    aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAAC‑‑TGGGGCTAAGGAGCGGATATGTTGTc                                                                   <  1:895650/80‑1 (MQ=255)
                                                                                      gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCa                                                                                                        >  2:1385008/1‑43 (MQ=255)
                                                                                      gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCa                                                                                                        <  1:1385008/43‑1 (MQ=255)
                                                                                      gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGTATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGctggctgg        >  1:275925/1‑139 (MQ=255)
                                                                                      gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                     >  2:1591766/1‑126 (MQ=255)
                                                                                      gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                     <  1:1591766/126‑1 (MQ=255)
                                                                                        cGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTa                                                                                     <  1:766571/60‑1 (MQ=255)
                                                                                        cGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTa                                                                                     >  2:766571/1‑60 (MQ=255)
                                                                                            tctcCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGcc  >  2:178115/1‑139 (MQ=255)
                                                                                                cGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                <  1:1247911/41‑1 (MQ=255)
                                                                                                cGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                >  2:1247911/1‑41 (MQ=255)
                                                                                                         |                                                                                                                             
TCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCC  >  NC_000913/4422688‑4422918

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