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)

TTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGT  >  NC_000913_3_pae_tpiA/4110599‑4110892
                                                                                                                                                     |                                                                                                                                                 
ttGCTGGCGACACAACGAAATCTACTGACACTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGCTAAGGCGAGGAGTTGAGGAAAGTAAGTGCCGGATATGA‑AACCCGGCACCTGTCAGACTTCGggcg                                                                                                                                                   >  2:297146/1‑149 (MQ=255)
            ccatgAAATCTTCTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCCGTGTGAGATTTTGCGTTAAGGCGAAGAGCCAAGGCAAGTAAGTGCCGGATATCA‑AATCCGGCTTCTGTCAGACTTAAGGCTGCCCCGCCGCGa                                                                                                                                       <  2:420238/145‑1 (MQ=255)
             aaCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                                      >  2:255993/1‑149 (MQ=255)
             aaCGAAATCTACGGACGGTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                                      <  2:158542/149‑1 (MQ=255)
              aCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGc                                                                                                                                     >  1:375218/1‑149 (MQ=255)
               cGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGGGTGAGATTTTGCGTTAAGGCGAAGAGAAACGCGAAGGAAG‑ACTGCATAAGAGATTCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCg                                                                                                                                    <  1:246550/149‑1 (MQ=255)
                  aaTCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGtt                                                                                                                                 >  1:498975/1‑149 (MQ=255)
                        cTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCa                                                                                                                           >  2:35847/1‑149 (MQ=255)
                           aCGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGAGCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATCTGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCAAGCGGTGTAGCCCAtt                                                                                                                        >  2:310899/1‑149 (MQ=255)
                               ttAACCGTGCTTCAGCCGAGGAAATTCAGGATCATAAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTaa                                                                                                                    <  1:223986/149‑1 (MQ=255)
                                                            gATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCCCCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAAGTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCa                                                                                       >  1:299650/1‑149 (MQ=255)
                                                              tcatcaGTGTGAGATTTTGCGTTAAGGCGACGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGACCCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAAATCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACg                                                                                     >  1:423803/1‑149 (MQ=255)
                                                                       tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                            <  1:322164/133‑1 (MQ=255)
                                                                       tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                            >  2:322164/1‑133 (MQ=255)
                                                                                 cGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAaca                                                                  <  1:255993/149‑1 (MQ=255)
                                                                                     aaGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTcc                                                              >  1:88799/1‑149 (MQ=255)
                                                                                                  aaGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGt                                                 <  1:198711/149‑1 (MQ=255)
                                                                                                     gAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGGCTTTa                                              >  2:460606/1‑149 (MQ=255)
                                                                                                                   gATATGA‑AATCCGGAACCTGTCAGACTTAAGGCTGCCCCGCCGCTAGCGGTTTCGCAAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTTGTGCCGAGTTTAAGCGCGGCTTCAAc                                <  2:370611/149‑1 (MQ=255)
                                                                                                                    atatGA‑AATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACg                               <  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)
                                                                                                                                                     |                                                                                                                                                 
TTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGA‑AATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCCCAGT  >  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: ‑