Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	140,015	Δ1 bp	coding (1211/2391 nt)	gcd ←	quinoprotein glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	140,013	0	C	.	100.0%	130.9 / NA	32	coding (1213/2391 nt)	gcd	quinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (15/17);  total (15/17)

CCGGTAGTGGCATTCAGCGCCAGAATCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTGCCCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGATCACGCCAGA  >  NC_000913/139874‑140142                                                                                                                                            |                                                                                                                                  ccGGTAGTGGCATTCAGCGCCAGAATCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTGc                                                                                                                                   <  2:149956/140‑1 (MQ=255)               cAGCGCCAGAATCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGa                                                                                                                    >  1:297790/1‑140 (MQ=255)                    ccAGAATCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGttta                                                                                                               <  1:220876/140‑3 (MQ=255)                      aGAATCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTtaaag                                                                                                             <  1:663256/140‑1 (MQ=255)                        aaTCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGtt                                                                                                                 >  2:447867/1‑134 (MQ=255)                        aaTCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGtt                                                                                                                 <  1:447867/134‑1 (MQ=255)                            gAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGc                                                                                                                      <  1:220805/125‑1 (MQ=255)                            gAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGc                                                                                                                      >  2:220805/1‑125 (MQ=255)                                 ggCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGgtgtg                                                                                                  <  1:136744/140‑1 (MQ=255)                                  gCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGgtgtgt                                                                                                 >  2:809796/1‑140 (MQ=255)                                            ccTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGAc                                                                                        >  2:574271/1‑139 (MQ=255)                                            ccTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGAc                                                                                        <  1:574271/139‑1 (MQ=255)                                               gTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGa                                                                                    <  2:62143/140‑1 (MQ=255)                                               gTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGa                                                                                    <  2:324025/140‑1 (MQ=255)                                                                                 cGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCaa                                                      <  2:53778/136‑1 (MQ=255)                                                                                 cGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCaa                                                      >  1:53778/1‑136 (MQ=255)                                                                                   tGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGccc                                                >  1:519192/1‑140 (MQ=255)                                                                                     gTCACGCCCATTGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCcac                                              >  2:66505/1‑140 (MQ=255)                                                                                     gTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGGAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCcac                                              >  2:549578/1‑140 (MQ=255)                                                                                               tCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTccc                                    >  2:426176/1‑140 (MQ=255)                                                                                                      aTAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGtt                             <  1:108344/140‑1 (MQ=255)                                                                                                           ccAGATCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACAt                        <  1:809796/140‑1 (MQ=255)                                                                                                               aTCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGgtg                               <  2:326383/129‑1 (MQ=255)                                                                                                               aTCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGgtg                               >  1:326383/1‑129 (MQ=255)                                                                                                               aTCCAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCaaa                    >  2:519700/1‑140 (MQ=255)                                                                                                                 ccAGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAAc                  >  1:440322/1‑140 (MQ=255)                                                                                                                   aGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCa                <  2:297790/140‑1 (MQ=255)                                                                                                                   aGCTTCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCa                <  1:426176/140‑1 (MQ=255)                                                                                                                       tCGCGTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGa            <  2:387307/140‑1 (MQ=255)                                                                                                                        cgcgTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGAt           <  2:702797/140‑1 (MQ=255)                                                                                                                         gcgTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGATc          >  1:503588/1‑140 (MQ=255)                                                                                                                          cgTCATAGGCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGATCa         >  1:414137/1‑140 (MQ=255)                                                                                                                                 ggCCGCTGGTG‑CCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGATCACGCCAGa  <  2:519192/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                  CCGGTAGTGGCATTCAGCGCCAGAATCGAGCTGGCATAACGTTCCTGTTCCGGTGTGCGGTTACCGCCCCAGATATCCGGCGTGGTCACGCCCATCGGCAGATAGACCAGATCCAGCTTCGCGTCATAGGCCGCTGGTGCCCAGGAGTTTGGCGAGTTAAAGGTAAAGGTGTGTTCGTCAGACGGGATTGCGTTCGGATCTTTCGCGCCGGGATCAAAAGCCCACAGCAGCTCCCCGGTGTTGACATCAAAACCACGGATCACGCCAGA  >  NC_000913/139874‑140142

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