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_tpiA	4,110,747	0	G	T	100.0%	105.4 / NA	33	intergenic (+62/+784)	cdh/yiiQ	CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/16);  total (17/16)

GCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACC  >  NC_000913_3_pae_tpiA/4110605‑4110883                                                                                                                                               |                                                                                                                                         gCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCccg                                                                                                                                    >  1:364229/1‑149 (MQ=255)   gACACAACGAAATCTACTGTCGCTGAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGGAAGTGCCGGATATGAAATCCGGCGCCTGTCAGACTTTAGGCTGCCccgcc                                                                                                                                  <  1:570032/149‑1 (MQ=255)        aaCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                             >  1:458166/1‑149 (MQ=255)         aCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATCTGAAAAAAGGCGCCGGTCAGACTTAAGGCTGCCCCGCCGCGAGc                                                                                                                            <  1:341800/149‑1 (MQ=255)                      aCGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGCGCGGTTTCGCCAAtt                                                                                                               >  2:19144/1‑149 (MQ=255)                      aCGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAAtt                                                                                                               >  2:235586/1‑149 (MQ=255)                        gCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGGGCGGTTTCGCCAATTCt                                                                                                             >  1:393473/1‑149 (MQ=255)                          ttAACCGTGCCTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTaa                                                                                                           <  1:393779/149‑1 (MQ=255)                               cGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCt                                                                                                      >  1:152018/1‑149 (MQ=255)                                gTGCTTCAGACAAGAAAATGTAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAAGTCTAAGATCtt                                                                                                     <  1:437796/149‑1 (MQ=255)                                 tGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTa                                                                                                    >  1:41814/1‑149 (MQ=255)                                       aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGGTTCGCCAATTCTAAGATCTTAAGatat                                                                                              >  1:558337/1‑149 (MQ=255)                                        ggcGAGGATATTCTGGCTCAGCAGTTGGAGATTTGTCGTTAAGGGGAAGAGTTAATGAAAGTAAGTGCCGGATATGAAAACGCGCACCTGTCAGACTGAAGTCTGTCCCGCCGCCGGCGGTTTGGCGAATTCTGAGATCTTAAGATATg                                                                                             <  2:568613/147‑1 (MQ=255)                                        gCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATg                                                                                             <  1:235586/149‑1 (MQ=255)                                         ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTACGATATgg                                                                                            >  1:562753/1‑149 (MQ=255)                                         ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATgg                                                                                            <  1:411177/149‑1 (MQ=255)                                         ccGAGGAAAGTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATgg                                                                                            <  1:316808/149‑1 (MQ=255)                                              gAAATTCAGGATCATCAGTGTGAGATTTGACGTTCGGGGGGAGGGGAAAGGGAAGTAAGCGCCGGGTAGGTCTTGGGGGACGTCGCAGACTTAAGACTGCCCCGACGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGAAATAGGTCtt                                                                                       >  2:37572/1‑149 (MQ=255)                                                  ttCAGGATAATCAGTGTGAGATTTTGCGTTAAGGTGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAACCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                   <  1:615441/149‑1 (MQ=255)                                                                  tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCtt                                                                   >  1:453994/1‑149 (MQ=255)                                                                  tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCtt                                                                   <  1:243989/149‑1 (MQ=255)                                                                    agaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGc                                                                 <  2:453994/149‑1 (MQ=255)                                                                          tGCGGTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCaa                                                           <  2:445988/149‑1 (MQ=255)                                                                                  ggCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTcccc                                                   <  1:40252/149‑1 (MQ=255)                                                                                   gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCt                                                  >  1:37546/1‑149 (MQ=255)                                                                                              aGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGtt                                       >  1:105881/1‑149 (MQ=255)                                                                                                gAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTa                                     >  1:409781/1‑149 (MQ=255)                                                                                                           ccGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTc                          >  2:261965/1‑149 (MQ=255)                                                                                                              gATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAAc                       <  1:18306/149‑1 (MQ=255)                                                                                                                tatGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCGGCGCTTGCCAGCAACACTCCCCTCGGGCGGAGTTTAAGCGCGGCTGCacggt                     >  2:616411/1‑145 (MQ=255)                                                                                                                        ccGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCTCCAATTCTAAGATCTTAAGATATGGGTCTTATTCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCg             <  2:424723/149‑1 (MQ=255)                                                                                                                            cACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTCCACGGCTGCGCTTGCCAGCAACACTCCCCCCGTGCCGAGTTTAAGCGCGCCTTCAACGTCCTCTCCCctct         <  2:268873/149‑1 (MQ=255)                                                                                                                                   cAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATAcc  >  1:650308/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                         GCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACC  >  NC_000913_3_pae_tpiA/4110605‑4110883

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