Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,621,188 A→G 8.4% Y85C (TAT→TGT)  zwf → glucose‑6‑phosphate 1‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,621,1880AG 8.4% 123.0 / 4.0 48Y85C (TAT→TGT) zwfglucose‑6‑phosphate 1‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base G (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

AACATATTGCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGCACCACAATTCTTTGGCGTTATTT  >  CP000730/1621056‑1621325
                                                                                                                                    |                                                                                                                                         
aaCATATTGCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGAc                                                                                                                                   >  1:694519/1‑141 (MQ=255)
     atTGCCATCATCGGTCTTGGACGTCGTGACATTCCTAGTGGTGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGTACGTAAACGATACAAAAAAAATTTACGCGTTTATGGAACATGTCTTCTATCATAGACatgat                                                                                                                              <  2:4198/141‑1 (MQ=255)
     atTGCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACatgat                                                                                                                              <  2:148743/141‑1 (MQ=255)
     atTGCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACatgat                                                                                                                              <  1:212177/141‑1 (MQ=255)
        gCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGtt                                                                                                                           >  1:46863/1‑141 (MQ=255)
        gCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGtt                                                                                                                           >  2:201905/1‑141 (MQ=255)
        gCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGtt                                                                                                                           >  1:342698/1‑141 (MQ=255)
             catcGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTaa                                                                                                                      <  2:632943/141‑1 (MQ=255)
               tcGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATg                                                                                                                    <  2:685478/141‑1 (MQ=255)
                  gTATTGGACGTCGTGACATTACTGATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATgaag                                                                                                                 <  2:407862/141‑1 (MQ=255)
                   tATTGGACGTCGGGACATTACTAATGGTGATCTCCGTAATCACGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACCACATTGACGCGGTTATGGAACATGTCTTCTATCATAGACATGCTGTTAGTACTgaaga                                                                                                                >  2:310188/1‑141 (MQ=255)
                    aTTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATgaagaa                                                                                                               <  1:532403/141‑1 (MQ=255)
                    aTTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATgaagaa                                                                                                               <  2:701050/141‑1 (MQ=255)
                        gACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATAAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGCTCATGGATAATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCt                                                                                                           >  2:631698/1‑141 (MQ=255)
                         aCGTCGTGACATTACTAATGATGATTTCCGGAATAAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGCAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTa                                                                                                          >  1:462697/1‑141 (MQ=255)
                            tcgtGACATTACTAATGATGATTTCCGTAACCAAGTAAAATCATCAATTCAAACGCACGTACCAGCTACAACCACTCTTGCCGCGTTCTCGCAGCCAGCCTTATCTCCTCGCCGTGATGTTTGTAATGAAGAAAGCTATCa                                                                                                       <  1:278323/141‑1 (MQ=255)
                               tGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGa                                                                                                    >  2:369512/1‑141 (MQ=255)
                                gaaattagcaaTTATGACTTACGGATTCAGGGAAAATCATAAATTCAAAAGCACGTAAAAGCTACACACAAAATCGCGGCGTTTAAAGAACATGTCGTCTATAATAGACATGAGGTGACGAGTAAAGAAAGGTGTCAGGaa                                                                                                   >  2:128224/10‑141 (MQ=255)
                                 aCATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTGTCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAAt                                                                                                  >  1:263864/1‑141 (MQ=255)
                                 aCATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGCTAGTAATGAAGAAAGCTagtcaacat                                                                                                  >  1:660005/1‑133 (MQ=255)
                                    ttACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTAc                                                                                               <  1:218755/141‑1 (MQ=255)
                                        tAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGa                                                                                           <  2:446325/141‑1 (MQ=255)
                                          atgatgATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGAtt                                                                                         <  1:197035/141‑1 (MQ=255)
                                            gatggtTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAGAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGAtttt                                                                                       >  1:581322/1‑141 (MQ=255)
                                                atTTCCGTAATCAAGTAAAATAATCAATTCAAAAGCGCGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTTTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTa                                                                                   <  2:515667/141‑1 (MQ=255)
                                                 tttCCGTAATCGAGTAAAATCATCAATTCATAAGCACGTACAAGATACATATAAAATGAACGCGTTAATTGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTCTCAAGAATTACTAGATTTTAGTaa                                                                                  >  2:668382/1‑141 (MQ=255)
                                                 tttCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTGTCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTAc                                                                                               <  1:232186/128‑1 (MQ=255)
                                                 tttCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTGTCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTAc                                                                                               >  2:232186/1‑128 (MQ=255)
                                                 tttCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTaa                                                                                  >  2:668494/1‑141 (MQ=255)
                                                   tCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATg                                                                                <  2:419618/141‑1 (MQ=255)
                                                              gTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGc                                                                     <  1:548868/141‑1 (MQ=255)
                                                                      atcaATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGa                                                             >  2:195436/1‑141 (MQ=255)
                                                                      atcaATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGa                                                             >  2:325236/1‑141 (MQ=255)
                                                                        caATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAAt                                                           >  1:451629/1‑141 (MQ=255)
                                                                         aaTTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAAtt                                                          <  1:306892/141‑1 (MQ=255)
                                                                                aaGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGt                                                   >  2:435533/1‑141 (MQ=255)
                                                                                aaGCACGTAAAAGATACAAACAAAATTGACGCGTAAATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGt                                                   >  2:82543/1‑141 (MQ=255)
                                                                                 aGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTa                                                  >  1:417081/1‑141 (MQ=255)
                                                                                    aCGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTGTCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATc                                               <  2:263864/141‑1 (MQ=255)
                                                                                               acaaacaaAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGActattctat                                    <  1:128224/141‑1 (MQ=255)
                                                                                                          tGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGc                         >  1:499608/1‑141 (MQ=255)
                                                                                                              gcgTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGcacca                     >  1:699320/1‑141 (MQ=255)
                                                                                                               cgTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCACTTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGcaccac                    <  1:432801/141‑1 (MQ=255)
                                                                                                               cgTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGcaccac                    <  2:382837/141‑1 (MQ=255)
                                                                                                                 tttATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGCACCACaa                  <  1:98695/141‑1 (MQ=255)
                                                                                                                            tGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGCACCACAATTCTTTGGCGt       <  1:275227/141‑1 (MQ=255)
                                                                                                                               cttctATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGCACCACAATTCTTTGGCGTTAt    >  2:370255/1‑141 (MQ=255)
                                                                                                                                 tctATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGCACCACAATTCTTTGGCGTTAttt  >  2:107145/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                         
AACATATTGCCATCATCGGTATTGGACGTCGTGACATTACTAATGATGATTTCCGTAATCAAGTAAAATCATCAATTCAAAAGCACGTAAAAGATACAAACAAAATTGACGCGTTTATGGAACATGTCTTCTATCATAGACATGATGTTAGTAATGAAGAAAGCTATCAAGAATTACTAGATTTTAGTAATGAATTAGATAGCCAATTTGAATTAAAAGGTAATCGACTATTCTATTTAGCAATGGCACCACAATTCTTTGGCGTTATTT  >  CP000730/1621056‑1621325

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