Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,935,635 T→G D227A (GAT→GCT)  zwf ← glucose‑6‑phosphate 1‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,935,6350TG98.2% 211.9 / ‑4.1 56D227A (GAT→GCT) zwfglucose‑6‑phosphate 1‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (30/25);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTCGATGCGGCGCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAA  >  NC_000913/1935490‑1935779
                                                                                                                                                 |                                                                                                                                                
gTCGATGCGGCGCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCaaaa                                                                                                                                             <  2:721437/151‑1 (MQ=255)
 tCGATGCGGCGCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAAt                                                                                                                                            >  1:847891/1‑151 (MQ=255)
  cGATGCGGCGCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATa                                                                                                                                           <  1:74709/151‑1 (MQ=255)
    aTGCGGCGCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGc                                                                                                                                         >  1:241873/1‑151 (MQ=255)
      gcggcgCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGccc                                                                                                                                       <  2:401913/151‑1 (MQ=255)
       cggcgCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGcccc                                                                                                                                      >  2:464003/1‑151 (MQ=255)
        ggcgCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCa                                                                                                                                     <  2:889068/151‑1 (MQ=255)
                            ttCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCAc                                                                                                                   <  2:194737/149‑1 (MQ=255)
                             tCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCActt                                                                                                                 <  2:177310/150‑1 (MQ=255)
                                  ttttCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACttcttccg                                                                                                           >  1:297335/1‑149 (MQ=255)
                                      cATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCAc                                                                                                       >  2:462220/1‑151 (MQ=255)
                                      cATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCAc                                                                                                       >  2:799329/1‑151 (MQ=255)
                                      cATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCAc                                                                                                       >  1:52019/1‑151 (MQ=255)
                                              aTGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAAtct                                                                                               >  1:196675/1‑151 (MQ=255)
                                                gCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAAtctc                                                                                              >  2:403850/1‑150 (MQ=255)
                                                   gTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAAc                                                                                           >  2:507298/1‑150 (MQ=255)
                                                       gcgcTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATg                                                                                        <  2:347778/149‑1 (MQ=255)
                                                        cgcTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATc                                                                                     >  1:599963/1‑151 (MQ=255)
                                                         gcTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATc                                                                                     >  1:846736/1‑150 (MQ=255)
                                                           tCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAAt                                                                                  >  2:597637/1‑151 (MQ=255)
                                                            cAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAAt                                                                                  <  2:771592/150‑1 (MQ=255)
                                                                      ggcggAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGt                                                                       >  1:748819/1‑151 (MQ=255)
                                                                         ggAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTccc                                                                    >  1:760754/1‑151 (MQ=255)
                                                                          gagaCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCa                                                                   <  2:259095/151‑1 (MQ=255)
                                                                          gagaCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCAATTGTCCCa                                                                   <  1:210837/151‑1 (MQ=255)
                                                                           agaCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAg                                                                  <  1:541940/151‑1 (MQ=255)
                                                                            gaCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGt                                                                 <  2:513199/151‑1 (MQ=255)
                                                                             aCATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGtt                                                                >  2:326225/1‑151 (MQ=255)
                                                                                 cGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTc                                                            <  1:292218/151‑1 (MQ=255)
                                                                                    aaTCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTcac                                                          >  2:846783/1‑150 (MQ=255)
                                                                                    aaTCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTcac                                                          >  1:319277/1‑150 (MQ=255)
                                                                                        aTGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAAc                                                      >  2:408506/1‑150 (MQ=255)
                                                                                         tGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAAc                                                      >  2:310746/1‑149 (MQ=255)
                                                                                           cAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAggg                                                  <  1:234015/151‑1 (MQ=255)
                                                                                             aaGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCAGTTGACCTGCTTTATCAAAATATCCCCAGCGCCCTTCGAGCCCCACTTCTTCTTCCACGGTAAACTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGa                                                 <  2:732557/150‑1 (MQ=255)
                                                                                                aaTTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTa                                             <  1:647670/151‑1 (MQ=255)
                                                                                                   ttGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGc                                           <  1:260113/150‑1 (MQ=255)
                                                                                                            gtggtTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAgc                                 <  2:89922/151‑1 (MQ=255)
                                                                                                              ggtTCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAgcg                                <  2:241868/150‑1 (MQ=255)
                                                                                                                ttCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGcc                              >  2:479032/1‑150 (MQ=255)
                                                                                                                 tCTGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCa                             >  2:805211/1‑150 (MQ=255)
                                                                                                                   tGGATCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAAc                           >  1:403117/1‑150 (MQ=255)
                                                                                                                      aTCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGt                       <  1:513026/151‑1 (MQ=255)
                                                                                                                      aTCATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGt                       >  2:839995/1‑151 (MQ=255)
                                                                                                                        cATGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTc                     <  1:635846/151‑1 (MQ=255)
                                                                                                                         aTGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTc                     <  2:133905/150‑1 (MQ=255)
                                                                                                                          tGTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTc                     >  1:443751/1‑149 (MQ=255)
                                                                                                                           gTCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGc                  <  1:669888/151‑1 (MQ=255)
                                                                                                                            tCGCGCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGc                  >  1:643547/1‑150 (MQ=255)
                                                                                                                             cgcgCATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCa                 <  2:245430/150‑1 (MQ=255)
                                                                                                                                gcATCTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGt             >  2:739120/1‑151 (MQ=255)
                                                                                                                                    cTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCt         >  2:546530/1‑151 (MQ=255)
                                                                                                                                    cTGACCGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCt         >  2:415712/1‑151 (MQ=255)
                                                                                                                                        ccGGCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTAc     <  2:826768/151‑1 (MQ=255)
                                                                                                                                           gCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCaa  <  1:689584/151‑1 (MQ=255)
                                                                                                                                           gCTTTAGCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCa   >  1:638040/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTCGATGCGGCGCAGAGACTTCAGTACTTTCACTTTTTCATCGCGGATGCTGTCTGCGCTCAGGTCAGACGGCGGAGACATCGCAATCATGCAAAGAATTTGCAGCAGGTGGTTCTGGATCATGTCGCGCATCTGACCGGCTTTATCAAAATAGCCCCAGCGCCCTTCGATCCCCACTTCTTCTGCCACGGTAATCTCAACATGATCAATGGTGCGATTGTCCCAGTTATTCACAAACAGGGAGTTAGCAAAACGCAGCGCCAACAGGTTCAGCACCGTTTCTTTACCAA  >  NC_000913/1935490‑1935779

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