Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	140,924	C→T	G101D (GGC→GAC)	gcd ←	quinoprotein glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	140,924	0	C	T	96.4%	88.0 / NA	28	G101D (GGC→GAC)	gcd	quinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (14/13);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCGGATCGTTAAATCCGGCCCAGGTCAGGATACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGCCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCAC  >  NC_000913/140793‑141055                                                                                                                                    |                                                                                                                                    ccTGCGGATCGTTAAATCCGGCCCAGGTCAGGATACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCgaag                                                                                                                                <  1:119249/137‑1 (MQ=255) ccTGCGGATCGTTAAATCCGGCCCAGGTCAGGATACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCgaag                                                                                                                                >  2:119249/1‑137 (MQ=255)    gCGGATCGTTAAATCCGGCCCAGGTCAGGATACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGAcc                                                                                                                          >  1:140215/1‑140 (MQ=255)              aaaTCCGGCCCAGGTCAGGATACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGc                                                                                                                <  2:140215/140‑1 (MQ=255)                                  tACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAg                                                                                            <  2:535157/140‑1 (MQ=255)                                     caccGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGATGACCAGACTGTAGCCGCGCGGCGTCAGCGCCCAGAAg                                                                                         >  1:370808/1‑140 (MQ=255)                                     caccGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAAATGTCGAACAAGACAACAATGTCGCAGCGCGGAGTCAGCGCCCAGAAg                                                                                         >  2:386408/1‑140 (MQ=255)                                          cTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGaa                                                                                    >  1:854824/1‑140 (MQ=255)                                             aTCAGCAGTGCGACCACCAGTGGGGCAACTGGGCCGCTGTCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAAcc                                                                                 <  2:370808/140‑1 (MQ=255)                                                  cagTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACtt                                                                            >  1:692603/1‑140 (MQ=255)                                                   agTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTc                                                                           >  2:146033/1‑140 (MQ=255)                                                      gCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCa                                                                        <  1:659013/140‑1 (MQ=255)                                                        gACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTTGAAACCAACTTCCCAGa                                                                      <  2:56188/140‑1 (MQ=255)                                                            accaGTGCGGCAACTGCTCCGCCGGCAGGAATGACCAGGAGACGCCAGGCAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGcc                                                                  <  2:254707/140‑1 (MQ=255)                                                            accaGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGcc                                                                  >  2:511182/1‑140 (MQ=255)                                                                      cAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACGAAGTTCCCAGACGCCCCAAATCATg                                                        >  1:254707/1‑140 (MQ=255)                                                                         cTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTg                                                     <  2:550754/140‑1 (MQ=255)                                                                                 cTGGCAGGAATGACCCGGCGACACCAAACCAACGGCATGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGagcag                                             <  1:321492/140‑1 (MQ=255)                                                                                        gAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGCCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCt                                      >  2:794218/1‑140 (MQ=255)                                                                                              ccAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGtat                                <  1:705714/140‑1 (MQ=255)                                                                                                  gCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGtat                                >  1:827501/1‑136 (MQ=255)                                                                                                  gCGACGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGtat                                <  2:827501/136‑1 (MQ=255)                                                                                                      cGCCAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAACTCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAg                        >  1:373714/1‑140 (MQ=255)                                                                                                         cAGACAAACGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAgcgc                     <  2:854824/140‑1 (MQ=255)                                                                                                                aCGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGt              >  1:492626/1‑140 (MQ=255)                                                                                                                 cGGCAGGATCAGCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGtt             >  2:123672/1‑140 (MQ=255)                                                                                                                           aGCCAGATGTCGAAGAAGACCACAAGGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCca   <  1:317726/140‑1 (MQ=255)                                                                                                                            gCCAGATGTCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGTGCTGCGTATAGCCAAAACGCGACGCGTTTACTGCGCcac  >  2:252516/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                    CCTGCGGATCGTTAAATCCGGCCCAGGTCAGGATACCACCGCTAATCAGCAGTGCGACCACCAGTGCGGCAACTGCGCCGCTGGCAGGAATGACCAGGCGACGCCAGACAAACGGCAGGATCAGCCAGATGCCGAAGAAGACCAGAATGTCGCTGCGCGGAGTCAGCGCCCAGAAGTCGAAACCAACTTCCCAGACGCCCCAAATCATGGTGCCGAGCAGCAGGGCTGCGTATAGCCAAAGCGCGGCGCGTTTACTGCGCCAC  >  NC_000913/140793‑141055

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