Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	139,187	T→C	N680S (AAT→AGT)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	139,187	0	T	C	96.8%	204.9 / ‑2.0	62	N680S (AAT→AGT)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (30/30);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAC  >  NC_000913/139090‑139282                                                                                                  |                                                                                                cAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTa                                                                                               <  1:848567/100‑1 (MQ=255)   gCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGt                                                                                             <  2:1263761/100‑1 (MQ=255)    cATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGtt                                                                                            >  1:647297/1‑100 (MQ=255)       cGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTc                                                                                         <  2:401264/100‑1 (MQ=255)        ggCATACCCATATTGAACGGCACCGGAACTGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCATTACGTTTCTTCCACACCACTTCACTAGTTTTCa                                                                                        <  2:782241/100‑1 (MQ=255)         gCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAg                                                                                       >  1:200673/1‑100 (MQ=255)         gCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAg                                                                                       <  2:690728/100‑1 (MQ=255)            tACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATc                                                                                    >  1:121675/1‑100 (MQ=255)            tACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATc                                                                                    <  2:149197/100‑1 (MQ=255)              cccATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCa                                                                                  >  1:821635/1‑100 (MQ=255)                 atatTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAgcg                                                                               <  1:737191/100‑1 (MQ=255)                   atTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGcc                                                                             >  2:741131/1‑100 (MQ=255)                      gAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCATGCGGCTTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCga                                                                           >  2:300133/1‑99 (MQ=255)                       accGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCgaga                                                                         >  1:885646/3‑100 (MQ=255)                       aaCGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCgaga                                                                         <  1:1242301/100‑1 (MQ=255)                         cGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGata                                                                       >  2:686287/1‑100 (MQ=255)                         cGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGata                                                                       >  1:1263301/1‑100 (MQ=255)                             aCCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAAc                                                                   >  1:761738/1‑100 (MQ=255)                                ggAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAAcccc                                                                >  2:824700/1‑100 (MQ=255)                                 gAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCa                                                               <  2:1032988/100‑1 (MQ=255)                                      ggCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctg                                                          >  1:1211296/1‑100 (MQ=255)                                        cATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctggc                                                        <  1:631755/100‑1 (MQ=255)                                           ctggAATGTCATACTGTCCTGCGGTATACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGt                                                     <  2:1190725/98‑1 (MQ=255)                                             ggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGttt                                                   >  1:258381/1‑100 (MQ=255)                                              gAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTa                                                  <  2:463241/100‑1 (MQ=255)                                               aaCGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTAc                                                 <  1:775237/100‑1 (MQ=255)                                                aCGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACa                                                <  1:271664/100‑1 (MQ=255)                                                 cGGCATACTGTCTTGCGGCGTACCAATACGTTTCTTCCACACCACTTCTCTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACa                                                <  2:254813/99‑1 (MQ=255)                                                  ggCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATg                                              <  2:1124841/100‑1 (MQ=255)                                                  ggCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATg                                              <  2:553044/100‑1 (MQ=255)                                                  ggCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATg                                              >  1:1237620/1‑100 (MQ=255)                                                   gCATACTGTCCTGCTGCGTACCAATACGTTTCTTCGACGCAACGGCTCAAGTATTCTGATTGAGCGCCGGGATAGACCCCCATGCTGGCTGTTTACATg                                              <  2:931876/99‑1 (MQ=255)                                                    cATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGc                                            <  2:1168700/100‑1 (MQ=255)                                                    cATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGc                                            <  2:669217/100‑1 (MQ=255)                                                    cATACTGTCCTGCGGCGTACCAATACGTGTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGc                                            <  2:473150/100‑1 (MQ=255)                                                         tGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACGTGGCGGAcc                                       >  2:1106652/1‑100 (MQ=255)                                                         tGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGAcc                                       >  1:860150/1‑100 (MQ=255)                                                          gTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCa                                      >  1:692793/1‑100 (MQ=255)                                                           tCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCaa                                     >  2:960140/1‑100 (MQ=255)                                                            ccTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCaaa                                    >  1:1233893/1‑100 (MQ=255)                                                             cTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCaa                                     >  2:383364/1‑98 (MQ=255)                                                                 ggcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTg                               <  1:470503/100‑1 (MQ=255)                                                                  gcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              <  1:50104/100‑1 (MQ=255)                                                                   cgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTgag                             >  1:219066/1‑100 (MQ=255)                                                                   cgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              <  2:885639/99‑1 (MQ=255)                                                                     tACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGagag                           >  1:355661/1‑100 (MQ=255)                                                                      aCCAATACGTTTCTTCGACACCTCTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTgag                             <  2:541176/97‑1 (MQ=255)                                                                       ccAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGa                         >  2:447284/1‑100 (MQ=255)                                                                             cGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGt                   >  2:126185/1‑100 (MQ=255)                                                                              gTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGtt                  >  2:734350/1‑100 (MQ=255)                                                                                   ttcCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCg             >  1:1174649/1‑100 (MQ=255)                                                                                       acacCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGac         <  2:469617/100‑1 (MQ=255)                                                                                       acacCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGac         <  2:652651/100‑1 (MQ=255)                                                                                       acacCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGGTGAGCGTGac         >  1:911193/1‑100 (MQ=255)                                                                                          ccacTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACAcc      >  1:1134828/1‑100 (MQ=255)                                                                                          ccacTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACAcc      >  1:236035/1‑100 (MQ=255)                                                                                           cacTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCa     >  1:592151/1‑100 (MQ=255)                                                                                            acTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCAt    >  2:393797/1‑100 (MQ=255)                                                                                             cTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCa     <  2:560296/98‑1 (MQ=255)                                                                                              ttCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAc  <  2:502241/100‑1 (MQ=255)                                                                                              ttCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAc  <  1:741136/100‑1 (MQ=255)                                                                                              ttCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAc  <  2:355658/100‑1 (MQ=255)                                                                                                  |                                                                                                CAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAC  >  NC_000913/139090‑139282

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