Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,934,854 C→A 34.2% W487C (TGG→TGT)  zwf ← NADP(+)‑dependent glucose‑6‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,934,8540CA34.2% 29.3 / 34.7 49W487C (TGG→TGT) zwfNADP(+)‑dependent glucose‑6‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (16/15);  new base A (8/8);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

TTCAGTCAGTGTAATAAAAAAAGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGTCTACCCATTTCCAGGCTTCTTCCA  >  NC_000913/1934712‑1934990
                                                                                                                                               |                                                                                                                                        
ttCAGTCAGGGTAATAAAAAAAGGCTCGTGGGGGGATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCCACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTc                                                                                                                                          <  2:721510/143‑1 (MQ=255)
     tcagtGTAATAAAAAAAGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGCAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAAGATAAGCGCAGATATTACCCAAACCC‑ATTcacgaa                                                                                                                                     >  2:764310/1‑138 (MQ=255)
          gtAATAAAAAACG‑CGCGTTGGTGAGCCCGCACGAGGCCCGAAAGTGTAAAAATTGTACAAAAATCTGCGCAAGAGCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGAcc                                                                                                                               <  2:128603/143‑1 (MQ=255)
          gtAATAAAAAAAGCCTCGTGGGTGAATCCGCAGGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGAc                                                                                                                                <  2:13057/143‑1 (MQ=255)
           gcttaaAAATAGTCTCGTGGGTCAATCCGCACGAGGCCTGAACGTGTAACGATGGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAACAACCACAAAGGCGAAGCGCAGATATTACTCAAACTC‑ATTTCAGGAACGAcc                                                                                                                               <  2:114898/140‑1 (MQ=255)
            aataaAAAGAGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCa                                                                                                                              <  2:361970/143‑1 (MQ=255)
             ataaAAAAAGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATTATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATGACAGGaagactcag                                                                                                                             >  1:799345/1‑136 (MQ=255)
               aaaaaaaGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAGGGGTAAGCGCAGATATTAGTCAAACTC‑ATTCCAGGAACGAcggtaa                                                                                                                           >  2:514509/1‑138 (MQ=255)
                aaaaaaGCCTCGTGGGTGAATCCGCACGACGCCTGAAAGTGTAAAAATTGTCCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGACAAGCGCAGATAATACCCAAACTC‑ATTCCAGGACCGCCCATCAc                                                                                                                          >  2:200887/1‑143 (MQ=255)
                 aaaaaGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATACTACTCAAACCC‑ATCCCACGAACAGCCATCACg                                                                                                                         >  2:973121/1‑143 (MQ=255)
                         tCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAatca                                                                                                                 >  2:189592/1‑143 (MQ=255)
                             gtgTGAATCCGTACGAGGCCTGAAAGTGTAAAGAGTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAGTAGCCATATAGTATAAGCGCAGATATTACTCCAACTC‑ATTCCAGGAACGACCATCACGTGTAATCATCGc                                                                                                             <  2:334968/141‑1 (MQ=255)
                                 gAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCAcc                                                                                                         >  2:448094/1‑143 (MQ=255)
                                    tCCGGACGAGGCATGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAg                                                                                                      <  1:136626/143‑1 (MQ=255)
                                         aCGAGGCCTGAAAGTGTAAAAATTGTTCAACAAGCTGCGCAAGACCATGTTACCGGTAAAATAACCATAAAGGAGAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAAc                                                                                                 <  1:371948/143‑1 (MQ=255)
                                              gCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTCCCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGCTAAACATCGCCCCCGAGGCAACgcgac                                                                                            >  2:49382/1‑139 (MQ=255)
                                               ccTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCAAGTTACCGGGAAAATAACCATAAAGGATAAGCGCAGATAGTACTCAAACCC‑ATTACAGGAACGACCATCACGCGCAGTCATCGCCACACAGGCAAcagcaca                                                                                           >  2:406255/1‑137 (MQ=255)
                                                  gAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATATACGGAACGACCATCACGGGTAATCATCGCAACCAAGGCAACGGGTCCCCa                                                                                        >  1:863184/1‑143 (MQ=255)
                                                   aacgtgtAAAAATTGTTCTACAATCTTCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATGCCAACAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAg                                                                                       <  2:477608/140‑1 (MQ=255)
                                                   aaaGTGTATGAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAg                                                                                       <  2:1094912/143‑1 (MQ=255)
                                                   aaaGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAg                                                                                       >  2:701192/1‑143 (MQ=255)
                                                   aaaGTCTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAg                                                                                       <  2:911166/143‑1 (MQ=255)
                                                         tAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCg                                                                                 <  1:701192/143‑1 (MQ=255)
                                                           aaagTTGTACTACAATCTGCGCCAGAGCATGTTACCGGTAGAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑AGTCCAGGAACGACCATCACGGGTAAGCGTCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGc                                                                               <  2:672705/143‑1 (MQ=255)
                                                           aaaaTTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGc                                                                               <  1:189592/143‑1 (MQ=255)
                                                            aaaTTTTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGcc                                                                              <  2:879563/143‑1 (MQ=255)
                                                                  ttCTACAATCTGCGCGAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGCTCATCATCGCGACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATAc                                                                        <  1:572403/143‑1 (MQ=255)
                                                                        aaTCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCCCCGCTGTTAGACAAACCCAAAACCCTGTTTG‑CTGTCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTc                                                                  >  1:1102338/1‑143 (MQ=255)
                                                                              cgcAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGca                                                            <  2:868129/143‑1 (MQ=255)
                                                                              cgcAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGca                                                            >  1:687841/1‑143 (MQ=255)
                                                                               gcAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACACAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTGTGgcgcct                                                           >  1:469295/1‑141 (MQ=255)
                                                                                   gATCATGTTACCGGTAAAATAACCATAAAGGATCAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCAtt                                                       <  1:354451/143‑1 (MQ=255)
                                                                                     tCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGCTCCGGCCTGATCCGCATCCGGCGcatccctgc                                                     >  2:926979/1‑138 (MQ=255)
                                                                                                    aaTAACCATAAAGGATAAGCGCAGATATTACTCTATTTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCt                                      <  1:817289/143‑1 (MQ=255)
                                                                                                           aaagaGGACAAGCGCAGATATTACTCGAACGC‑ATTCCAGGAACGACCATCCCGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAAt                               <  1:780033/139‑1 (MQ=255)
                                                                                                           aTAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGCACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGCTCCGGCCTGCTCCGGGTGAGCCGCCCCATTGTCCATCGCCCACGCCCCAgcacc                               >  2:778009/1‑139 (MQ=255)
                                                                                                            tAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATg                              >  1:693529/1‑143 (MQ=255)
                                                                                                              aagtacaaGCGCAGAGATTACTCAAACTC‑ATTCCAGGAACGACCAGCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGa                            <  2:946493/137‑1 (MQ=255)
                                                                                                              aagcacaaGCGCAGCTAGCACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCCTCATTGTCCATCGCCCACGCCTCAGTAATGGa                            <  1:89565/137‑1 (MQ=255)
                                                                                                                    aaGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGTCTAc                      <  1:764310/143‑1 (MQ=255)
                                                                                                                    aaGCGCAGATATTACTCAAACTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAACGGAGTCTAc                      >  1:903980/1‑143 (MQ=255)
                                                                                                                      gcgcAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGGCTAccc                    >  1:68473/1‑143 (MQ=255)
                                                                                                                       cgcAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGAGCATCATTGTCCATCGCCCACGCCTCCGTTATGGAGTCTACCCa                   >  1:485898/1‑143 (MQ=255)
                                                                                                                       cgcAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCAACCGAGGCACCGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCGTCGCCCACGCCTCAGGAATGGAGTCTAcccc                   >  1:17605/1‑142 (MQ=255)
                                                                                                                               ttACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGCTACGGTTTCGGCGCATCCTTGTCCATCGCCCACGCCTCAGTCATGGAGTCTACCCATTTCCAgg           >  1:925033/1‑143 (MQ=255)
                                                                                                                                   tCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTACTGGAGTCTACCCATTTCAAGGcttg       >  1:504582/1‑142 (MQ=255)
                                                                                                                                     taaCTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGTCTACCCATTTCCAGGcttctt     <  1:873023/142‑1 (MQ=255)
                                                                                                                                     aaaCTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGCCTACCCATTTCCAGGcttctt     >  1:317025/1‑143 (MQ=255)
                                                                                                                                       aCTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGCTACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGGAATGGAGTCTACCCATTTCCAGGcttccgga   >  2:873023/1‑139 (MQ=255)
                                                                                                                                        cTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGTCTACCCATTTCCAGGCTTCTTCCa  >  1:725541/1‑143 (MQ=255)
                                                                                                                                        cTC‑ATTACAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGTCTACCCATTTCCAGGCTTCTTCCa  <  2:908244/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
TTCAGTCAGTGTAATAAAAAAAGCCTCGTGGGTGAATCCGCACGAGGCCTGAAAGTGTAAAAATTGTTCTACAATCTGCGCAAGATCATGTTACCGGTAAAATAACCATAAAGGATAAGCGCAGATATTACTCAAACTC‑ATTCCAGGAACGACCATCACGGGTAATCATCGCCACCGAGGCAACGGGTCCCCAGGTTCCGGCCTGATACGGTTTCGGCGCATCATTGTCCATCGCCCACGCCTCAGTAATGGAGTCTACCCATTTCCAGGCTTCTTCCA  >  NC_000913/1934712‑1934990

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