Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,056,273 G→A P11P (CCG→CCA)  fau → 5‑formyltetrahydrofolate cyclo‑ligase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,056,2730GA100.0% 149.5 / NA 46P11P (CCG→CCA) fau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/26);  total (20/26)

ACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAAC  >  NC_000913/3056138‑3056406
                                                                                                                                       |                                                                                                                                     
aCAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaa                                                                                                                                    >  2:1038096/1‑139 (MQ=255)
   gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                                 <  2:800692/139‑1 (MQ=255)
   gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaa                                                                                                                                 <  1:1196744/139‑1 (MQ=255)
     cTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAacaaca                                                                                                                               >  2:143060/1‑139 (MQ=255)
             gcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaac                                                                                                                                >  2:736268/1‑130 (MQ=255)
             gcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaacaac                                                                                                                                <  1:736268/130‑1 (MQ=255)
                            ccTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCg                                                                                                        >  1:105099/1‑139 (MQ=255)
                                 ttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                                   <  2:611475/139‑1 (MQ=255)
                                 ttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                                   <  1:559696/139‑1 (MQ=255)
                                  tATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAccc                                                                                                  <  1:555400/139‑1 (MQ=255)
                                      tGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgat                                                                                              >  1:1163342/1‑139 (MQ=255)
                                          aCCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGatgatg                                                                                          >  2:901779/1‑139 (MQ=255)
                                           ccAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGAtgatga                                                                                         >  2:929681/1‑139 (MQ=255)
                                              gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACtt                                                                                      >  1:622903/1‑139 (MQ=255)
                                                   gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                                 <  1:901779/139‑1 (MQ=255)
                                                      gCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCg                                                                              <  2:1059350/139‑1 (MQ=255)
                                                         aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGa                                                                                                                         >  2:403525/1‑93 (MQ=255)
                                                         aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGa                                                                                                                         <  1:403525/93‑1 (MQ=255)
                                                         aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtg                                                                           >  1:243845/1‑139 (MQ=255)
                                                             ccAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGa                                                                       <  1:1038096/139‑1 (MQ=255)
                                                                          tGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGt                                                          <  2:725382/139‑1 (MQ=255)
                                                                          tGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGt                                                          <  1:1167262/139‑1 (MQ=255)
                                                                          tGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGt                                                          <  1:929681/139‑1 (MQ=255)
                                                                           gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTc                                                         >  1:253810/1‑139 (MQ=255)
                                                                           gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGG‑AATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                                             >  2:666755/1‑118 (MQ=255)
                                                                           gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGG‑AATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCgg                                                                             <  1:666755/118‑1 (MQ=255)
                                                                            aCATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCg                                                        <  2:622903/139‑1 (MQ=255)
                                                                                tATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGa                                                                                                                         >  1:328886/1‑70 (MQ=255)
                                                                                tATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGa                                                                                                                         <  2:328886/70‑1 (MQ=255)
                                                                                 aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCg                                                        <  2:53366/134‑1 (MQ=255)
                                                                                 aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCg                                                        >  1:53366/1‑134 (MQ=255)
                                                                                 aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTa                                                   <  2:83496/139‑1 (MQ=255)
                                                                                  tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaac                                                                                                                                   >  2:188916/1‑58 (MQ=255)
                                                                                  tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGaac                                                                                                                                   <  1:188916/58‑1 (MQ=255)
                                                                                  tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAt                                                  >  1:1088007/1‑139 (MQ=255)
                                                                                                 cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGt                                                    >  1:82196/1‑122 (MQ=255)
                                                                                                 cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGt                                                    <  2:82196/122‑1 (MQ=255)
                                                                                                   cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                        <  2:672764/100‑1 (MQ=255)
                                                                                                   cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                        >  1:672764/1‑100 (MQ=255)
                                                                                                      aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGaa                              >  1:522576/1‑139 (MQ=255)
                                                                                                           ttCGGCAACGTCGTCGTGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGa                         <  2:105099/139‑1 (MQ=255)
                                                                                                                   cgtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGTTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc                 <  2:243845/139‑1 (MQ=255)
                                                                                                                    gtcgtcgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGcc                <  1:270754/139‑1 (MQ=255)
                                                                                                                         cgtGCGTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCa           <  2:522576/139‑1 (MQ=255)
                                                                                                                              gTTAACGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaa      <  2:1163342/139‑1 (MQ=255)
                                                                                                                                  aCGCCAGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaacaac  >  1:257457/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                     
ACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAAC  >  NC_000913/3056138‑3056406

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