Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,050,691 G→A intergenic (‑24/+424) gcvT ← / ← ubiI aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,6910GA100.0% 196.7 / NA 60intergenic (‑24/+424)gcvT/ubiIaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/34);  total (26/34)

GATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATC  >  NC_000913/3050556‑3050825
                                                                                                                                       |                                                                                                                                      
gatgTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGTGAGCGCCGCCAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTTTGCCATCTTGTCCTCATTGAATAAGCGGGACTg                                                                                                                                     >  1:784184/1‑139 (MQ=255)
gatgTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTg                                                                                                                                     <  2:1331129/139‑1 (MQ=255)
  tgTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGAc                                                                                                                                   <  2:63668/139‑1 (MQ=255)
                               cGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACg                                                                                                      <  2:1468879/139‑1 (MQ=255)
                                            ccGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCg                                                                                         <  2:1329597/139‑1 (MQ=255)
                                               tGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAg                                                                                      <  2:1608382/139‑1 (MQ=255)
                                               tGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAg                                                                                      >  1:1608382/1‑139 (MQ=255)
                                                ggAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGc                                                                                     <  1:757337/139‑1 (MQ=255)
                                                                aGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGaaa                                                                                                   <  2:1492566/109‑1 (MQ=255)
                                                                aGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGaaa                                                                                                   >  1:1492566/1‑109 (MQ=255)
                                                                aGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCaa                                                                     >  1:1138594/1‑139 (MQ=255)
                                                                     cgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCaa                                                                <  1:74685/139‑1 (MQ=255)
                                                                       cAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACt                                                                                                                                      >  1:1474431/1‑67 (MQ=255)
                                                                       cAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACt                                                                                                                                      <  2:1474431/67‑1 (MQ=255)
                                                                               gtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                      <  2:1138594/139‑1 (MQ=255)
                                                                                   tgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGc                                                  <  1:247397/139‑1 (MQ=255)
                                                                                   tgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGc                                                  <  1:1276778/139‑1 (MQ=255)
                                                                                     ttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCt                                                <  1:1581983/139‑1 (MQ=255)
                                                                                     ttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCt                                                >  2:653267/1‑139 (MQ=255)
                                                                                      tCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCAc                                                                                  <  2:652008/104‑1 (MQ=255)
                                                                                      tCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCAc                                                                                  >  1:652008/1‑104 (MQ=255)
                                                                                       cGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGa                                              >  2:1758876/1‑139 (MQ=255)
                                                                                         tACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGc                                                                  <  1:780378/117‑1 (MQ=255)
                                                                                         tACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGc                                                                  >  2:780378/1‑117 (MQ=255)
                                                                                          aCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACt                                           >  1:1130474/1‑139 (MQ=255)
                                                                                          aCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACt                                           >  2:867280/1‑139 (MQ=255)
                                                                                              aGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacc                                       >  2:110017/1‑139 (MQ=255)
                                                                                                       ttgtGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACa                              <  1:867280/139‑1 (MQ=255)
                                                                                                       ttgtGCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACa                              <  2:1130474/139‑1 (MQ=255)
                                                                                                         gtgCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGa                            <  1:110017/139‑1 (MQ=255)
                                                                                                         gtgCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGa                            <  2:1790129/139‑1 (MQ=255)
                                                                                                          tgCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                           >  2:1805534/1‑139 (MQ=255)
                                                                                                          tgCCATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACt                          >  2:841093/1‑139 (MQ=255)
                                                                                                            ccATCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTg                         >  2:295242/1‑139 (MQ=255)
                                                                                                               tCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCa                                                                                                                         <  1:33388/40‑1 (MQ=255)
                                                                                                               tCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCa                                                                                                                         >  2:33388/1‑40 (MQ=255)
                                                                                                               tCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacca                                      <  2:806904/123‑1 (MQ=255)
                                                                                                               tCTTGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGAAATCGTTCTCTTTTGCCTGAACTTacca                                      >  1:806904/1‑123 (MQ=255)
                                                                                                                  tGTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCa                  <  1:841093/139‑1 (MQ=255)
                                                                                                                   gTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAGCTTACCACCGAAACAGACTGTTAACCa                  <  2:784184/139‑1 (MQ=255)
                                                                                                                   gTCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCa                  <  1:975562/139‑1 (MQ=255)
                                                                                                                    tCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACg                                                                                                      >  1:732025/1‑54 (MQ=255)
                                                                                                                    tCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACg                                                                                                      <  2:732025/54‑1 (MQ=255)
                                                                                                                    tCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                   <  1:1138299/89‑1 (MQ=255)
                                                                                                                    tCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                   >  2:1138299/1‑89 (MQ=255)
                                                                                                                    tCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTaccac                                     >  2:809270/1‑119 (MQ=255)
                                                                                                                    tCCTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTaccac                                     <  1:809270/119‑1 (MQ=255)
                                                                                                                      cTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacca                                      <  2:1651917/116‑1 (MQ=255)
                                                                                                                      cTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacca                                      >  1:1651917/1‑116 (MQ=255)
                                                                                                                      cTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGa                            >  1:621025/1‑126 (MQ=255)
                                                                                                                      cTCATTGAATAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGa                            <  2:621025/126‑1 (MQ=255)
                                                                                                                              aTAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCg                                                                             <  2:597765/69‑1 (MQ=255)
                                                                                                                              aTAAGCGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCg                                                                             >  1:597765/1‑69 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                   >  1:1402256/1‑74 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                   <  2:1402256/74‑1 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGc                                                  <  1:1539120/91‑1 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGc                                                  >  2:1539120/1‑91 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc  >  2:496538/1‑139 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc  <  1:653267/139‑1 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc  >  1:352878/1‑139 (MQ=255)
                                                                                                                                   cGGGACTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc  <  1:496538/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATC  >  NC_000913/3050556‑3050825

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