Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,110,747 G→T intergenic (+62/+784) cdh → / ← yiiQ CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,110,7470GT100.0% 67.8 / NA 22intergenic (+62/+784)cdh/yiiQCDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/9);  total (13/9)

TTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGT  >  NC_000913_3_pae_tpiA/4110599‑4110892
                                                                                                                                                     |                                                                                                                                                 
ttGCTGGCGACACAACGAAATCTACTGACACTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGCTAAGGCGAGGAGTTGAGGAAAGTAAGTGCCGGATATGAAACCCGGCACCTGTCAGACTTCGggcg                                                                                                                                                   >  2:297146/1‑149 (MQ=255)
            ccatgAAATCTTCTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCCGTGTGAGATTTTGCGTTAAGGCGAAGAGCCAAGGCAAGTAAGTGCCGGATATCAAATCCGGCTTCTGTCAGACTTAAGGCTGCCCCGCCGCGa                                                                                                                                       <  2:420238/145‑1 (MQ=255)
             aaCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                                      >  2:255993/1‑149 (MQ=255)
             aaCGAAATCTACGGACGGTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                                      <  2:158542/149‑1 (MQ=255)
              aCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGc                                                                                                                                     >  1:375218/1‑149 (MQ=255)
               cGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGGGTGAGATTTTGCGTTAAGGCGAAGAGAAACGCGAAGGAAGACTGCATAAGAGATTCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCg                                                                                                                                    <  1:246550/149‑1 (MQ=255)
                  aaTCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGtt                                                                                                                                 >  1:498975/1‑149 (MQ=255)
                        cTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCa                                                                                                                           >  2:35847/1‑149 (MQ=255)
                           aCGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGAGCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATCTGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCAAGCGGTGTAGCCCAtt                                                                                                                        >  2:310899/1‑149 (MQ=255)
                               ttAACCGTGCTTCAGCCGAGGAAATTCAGGATCATAAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTaa                                                                                                                    <  1:223986/149‑1 (MQ=255)
                                                            gATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCCCCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAAGTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCa                                                                                       >  1:299650/1‑149 (MQ=255)
                                                              tcatcaGTGTGAGATTTTGCGTTAAGGCGACGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGACCCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAAATCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACg                                                                                     >  1:423803/1‑149 (MQ=255)
                                                                       tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                            <  1:322164/133‑1 (MQ=255)
                                                                       tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                            >  2:322164/1‑133 (MQ=255)
                                                                                 cGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAaca                                                                  <  1:255993/149‑1 (MQ=255)
                                                                                     aaGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTcc                                                              >  1:88799/1‑149 (MQ=255)
                                                                                                  aaGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGt                                                 <  1:198711/149‑1 (MQ=255)
                                                                                                     gAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGGCTTTa                                              >  2:460606/1‑149 (MQ=255)
                                                                                                                   gATATGAAATCCGGAACCTGTCAGACTTAAGGCTGCCCCGCCGCTAGCGGTTTCGCAAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTTGTGCCGAGTTTAAGCGCGGCTTCAAc                                <  2:370611/149‑1 (MQ=255)
                                                                                                                    atatGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACg                               <  1:35847/149‑1 (MQ=255)
                                                                                                                                          aGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGCCTTCTACGTCCTCTCCGCTCCCACTAccc          >  2:87975/1‑149 (MQ=255)
                                                                                                                                            aCTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGc        >  2:325397/1‑149 (MQ=255)
                                                                                                                                              ttAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCCCACTACCCGccc      >  2:20216/1‑149 (MQ=255)
                                                                                                                                                  ggctgcCCCGCCGCGGGCGGTTGCGCAAATTCTAAGCTCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCGCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATCCCCGCCCCAGt  <  2:79539/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGT  >  NC_000913_3_pae_tpiA/4110599‑4110892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: