Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 139,326 C→T G634S (GGC→AGC)  gcd ← glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913139,3260CT97.9% 180.6 / NA 47G634S (GGC→AGC) gcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/24);  total (22/25)
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

CGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCA  >  NC_000913/139153‑139529
                                                                                                                                                                              |                                                                                                                                                                                                            
cggcgTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                                                                                                                                                                       >  1:150446/1‑212 (MQ=255)
cggcgTACCAATACGTTTCTTCCACAACACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                                                                                                                                                                       <  2:150446/212‑1 (MQ=255)
    gTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAg                                                                                                                                        <  3:153958/239‑1 (MQ=255)
     tACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGAt                                                                                                                                                             <  2:153250/217‑1 (MQ=255)
     tACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGAt                                                                                                                                                             >  1:153250/1‑217 (MQ=255)
     tACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGt                                                                                                                                       >  3:341431/1‑239 (MQ=255)
           tACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAtt                                                                                                                                 <  3:141306/239‑1 (MQ=255)
                ttcttcCACACGCATTCATTGGTTTTCAGATCCAGCGCCGAGATATCACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACCCCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGt                                                                                                                             <  4:341431/238‑1 (MQ=255)
                   ttcCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  <  4:117017/182‑1 (MQ=255)
                   ttcCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  >  3:117017/1‑182 (MQ=255)
                                cATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCa                                                                                                            <  1:34381/239‑1 (MQ=255)
                                                             taACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACAACGCGGGATCAGTTTCGCAACAAACGGCAGTGCCATTGGGTTGGGAATCGTCACGTCACGATTTGGATCAACGGAAATCTCGtcc                                                                                 >  4:155009/1‑237 (MQ=255)
                                                                  ccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  <  3:230515/135‑1 (MQ=255)
                                                                  ccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  >  4:230515/1‑135 (MQ=255)
                                                                    aTGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAAc                                                                        >  1:95466/1‑239 (MQ=255)
                                                                          gctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCGGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATTAACGGAAATCCCGCCCCAATCGAAAATtccc                                                                  >  4:317991/1‑239 (MQ=255)
                                                                          gctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGCAATCCCGCCCCATTCGAACATcccc                                                                  >  3:292323/1‑239 (MQ=255)
                                                                               ttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAtt                                                                                                                                 <  2:228688/171‑1 (MQ=255)
                                                                               ttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAtt                                                                                                                                 >  1:228688/1‑171 (MQ=255)
                                                                                          aCCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAg                                                   >  3:162451/1‑238 (MQ=255)
                                                                                          aCCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGACCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTTACGATTTGGATCAACGGACATACAGCCGCATTCGTACAATCCCAGGTTACCCGGGAAgc                                                  >  2:67225/1‑238 (MQ=255)
                                                                                          aCCAAATGGTAGGAAGACGGGGTTGAGCGTAGACACAATCGGTACACCGTACTGTGGATAAATCCCGGTCTCCGTACCCGTCCTTTTGGCCTCTTTCGGCTGCTCCATCGGGTTTCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGa                                                  <  4:292323/239‑1 (MQ=255)
                                                                                           ccAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCGGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGa                                                  <  3:317991/238‑1 (MQ=255)
                                                                                               aTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTTAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGc                                             <  3:155009/239‑1 (MQ=255)
                                                                                                      gagGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGACCATCCCCAGGTTAACCGGGAAGACCAGCggaa                                         >  4:476304/1‑233 (MQ=255)
                                                                                                          aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTAccc                                       <  3:239800/234‑1 (MQ=255)
                                                                                                                 tGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATgg                            <  3:476304/238‑1 (MQ=255)
                                                                                                                     catgtacgccgtACGGTAAGCGGGACTGGGGCTGAATGGGGGATCCGGTACCGGGACTATGGACATCGTTCGGCTGCTCCAGCGGTGTGCCAGGACCACGCGGGATCAGTTGCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTGCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACAGCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGAggg                            <  4:325497/224‑1 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgct                                                                                                                                                                                            <  4:355815/70‑1 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgct                                                                                                                                                                                            >  3:355815/1‑70 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCAc                                                                                                                                                                    >  3:290504/1‑94 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCAc                                                                                                                                                                    <  4:290504/94‑1 (MQ=255)
                                                                                                                                        acCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         >  3:185673/1‑107 (MQ=255)
                                                                                                                                        acCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         <  4:185673/107‑1 (MQ=255)
                                                                                                                                           gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  <  3:283090/239‑1 (MQ=255)
                                                                                                                                            tACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         <  1:224139/103‑1 (MQ=255)
                                                                                                                                            tACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         >  2:224139/1‑103 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  <  3:323377/52‑1 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  >  4:323377/1‑52 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  >  1:5852/1‑228 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  <  2:5852/228‑1 (MQ=255)
                                                                                                                                                                        ccGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGa                                                                          >  3:318874/1‑138 (MQ=255)
                                                                                                                                                                        ccGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGa                                                                          <  4:318874/138‑1 (MQ=255)
                                                                                                                                                                          gTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  <  3:106574/208‑1 (MQ=255)
                                                                                                                                                                          gTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  >  4:106574/1‑208 (MQ=255)
                                                                                                                                                                           tGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATc                                                                     >  4:328242/1‑140 (MQ=255)
                                                                                                                                                                           tGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATc                                                                     <  3:328242/140‑1 (MQ=255)
                                                                                                                                                                              |                                                                                                                                                                                                            
CGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCA  >  NC_000913/139153‑139529

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: