Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	2,679,525	G→T	A1074D (GCT→GAT)	yphG ←	DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	2,679,525	0	G	T	100.0%	152.4 / NA	47	A1074D (GCT→GAT)	yphG	DUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (30/17);  total (30/17)

GGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGG  >  NC_000913_3_pae_tpiA/2679382‑2679667                                                                                                                                                |                                                                                                                                               ggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                                                                                                           <  1:251136/149‑1 (MQ=255) ggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                                                                                                           <  1:707533/149‑1 (MQ=255)     tttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                                                                                                       >  2:529872/1‑149 (MQ=255)       ccttaACAACTGCTTATTCTAATTAGTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATAACAGTTAATAGGGTTTAGTGAACGCGCTTAAAAATAAAAGACATCACGGAGGTGAATCATGCCGCTTTCGCTCAATCCTGACTAATc                                                                                                                                     >  2:141034/3‑149 (MQ=255)           aacaacAGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                                                                                                 <  1:202998/149‑1 (MQ=255)            acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATGTTGCAGATAATGGGGTTTAGTGAACGCGGTTCAAAAGAAAAGGCATCATGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTGATCTCCGt                                                                                                                                >  1:964563/1‑149 (MQ=255)                cGGCTTAGTGTAATTGTTTTTTGATGAGCGGTATGGCGTGCCGCTTAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAt                                                                                                                            <  2:828071/149‑1 (MQ=255)                 ggCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAtt                                                                                                                           >  2:281705/1‑149 (MQ=255)                          tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGac                                                                                                                  >  2:837921/1‑148 (MQ=255)                           aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGc                                                                                                                 <  1:529872/149‑1 (MQ=255)                             ttattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTACAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                                                                                                               >  1:1014443/1‑149 (MQ=255)                             ttattTTGTGATGAGCGGCAACGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                                                                                                               >  2:39813/1‑149 (MQ=255)                              gatgttGTGATGGGCGGCAGCGCGTGCCGGTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGTTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGGGCAGa                                                                                                              <  1:734633/145‑1 (MQ=255)                                     tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATGGCGATGGCGCTAATCGGTTTACACTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCATACAATTCAGCGCGAGCAGACTCatat                                                                                                       >  1:1123991/1‑149 (MQ=255)                                      gatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt                                                                                                      <  1:674220/149‑1 (MQ=255)                                      gatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTGAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt                                                                                                      >  1:892199/1‑149 (MQ=255)                                         gaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcg                                                                                                   >  1:298493/1‑149 (MQ=255)                                         gaGCGGCAGCGAGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcg                                                                                                   >  1:699643/1‑149 (MQ=255)                                             ggCAGCGCGTGCCGCTCAAATATTACAGGTAATAGGGTTTAATGTACTTCGTGATTACTAAGAGTGTGCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                               >  1:268851/1‑149 (MQ=255)                                              gCAGCGCGTGCCGCTCAAATAATACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATCAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAcc                                                                                              >  2:1007748/1‑149 (MQ=255)                                                        ccGCTCAAATATTACAGTTAATAGGGTTGAGTGAACGCGGGTAAAAATAAAAGGCATCACGGTGGTGAAGAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTc                                                                                    <  1:139562/149‑1 (MQ=255)                                                          gCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTccc                                                                                  >  1:766098/1‑149 (MQ=255)                                                              aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTTGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCgg                                                                              >  1:357541/1‑149 (MQ=255)                                                              aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCgg                                                                              >  1:459944/1‑149 (MQ=255)                                                              aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCgg                                                                              <  1:814467/149‑1 (MQ=255)                                                                   ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc                                                                         >  2:802027/1‑149 (MQ=255)                                                                    tACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCa                                                                        <  1:929685/149‑1 (MQ=255)                                                                        gTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAAGCCAGCAt                                                                    >  2:974573/1‑149 (MQ=255)                                                                         ttAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATc                                                                   >  1:288589/1‑149 (MQ=255)                                                                               gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTACGCCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGCCGCCAGCATCGtcacc                                                             >  1:420737/1‑149 (MQ=255)                                                                                ggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGccacca                                                            <  1:929720/149‑1 (MQ=255)                                                                                             cgGTGAAAAAGCACAGGCACCCCGGTGGTGAATAATGCCGCTTTCGGCAAATGCGGATTAATCGCCGTCAATTCAGCACAAATAAACTGATATTGCGCAGCCGCCCCGGGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTc                                               <  1:503635/149‑1 (MQ=255)                                                                                                     aaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgct                                       <  1:281705/149‑1 (MQ=255)                                                                                                        aaaaGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTCATCGCCGTCAATTAAGCGCGAGCAGACTCATATTGGGCAACCTCCCAGAGTCCCCCGGACACCAGCATAGCCAACATCAGGCATGTTTCGTTATgctgct                                    >  1:725982/1‑149 (MQ=255)                                                                                                           aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgtt                                 >  1:44123/1‑149 (MQ=255)                                                                                                               aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCCGCATCGCCACCATCAGGCATTTTCATTTATGCCGCTGTTGTaa                             >  2:517768/1‑149 (MQ=255)                                                                                                                           aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCCTTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAg                 >  2:1045903/1‑149 (MQ=255)                                                                                                                           aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCACGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAg                 >  1:892861/1‑149 (MQ=255)                                                                                                                           aataatGCCGCTTTCGGCCAATCAGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAg                 >  1:848831/1‑149 (MQ=255)                                                                                                                             taatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                               <  1:1057050/69‑1 (MQ=255)                                                                                                                             taatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                               >  2:1057050/1‑69 (MQ=255)                                                                                                                               atGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGac             <  1:968235/149‑1 (MQ=255)                                                                                                                                tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCgcaa            >  2:714466/1‑146 (MQ=255)                                                                                                                                tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGaca            >  1:467608/1‑149 (MQ=255)                                                                                                                                    gCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa        >  1:197562/1‑149 (MQ=255)                                                                                                                                      tttCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGt      <  1:483648/149‑1 (MQ=255)                                                                                                                                          ggCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAgg  <  2:197562/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               GGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGG  >  NC_000913_3_pae_tpiA/2679382‑2679667

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: ‑