Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,492,250 T→C 7.8% G266G (GGT→GGC)  yhfP → possible zinc (Zn2+) dependent dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,492,2500TC 7.8% 133.6 / 3.3 51G266G (GGT→GGC) yhfPpossible zinc (Zn2+) dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (2/2);  total (26/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 = 8.08e-01

AAGCTTGTGTAGACCCTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTCCAGA  >  CP000730/2492111‑2492387
                                                                                                                                           |                                                                                                                                         
aaGCTTGTGTAGACCCTGTTGGTGGGGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATCCTAATTCTGTATTCCCTCACATTTTAAGAGgtg                                                                                                                                          <  2:132334/141‑1 (MQ=255)
  gCTTGTGTAGACCCTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTa                                                                                                                                        <  1:453375/141‑1 (MQ=255)
  gCTTGTGTAGACCCTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTa                                                                                                                                        <  1:265164/141‑1 (MQ=255)
   cTTGTGTAGACCCTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTaa                                                                                                                                       <  1:162807/141‑1 (MQ=255)
     tgtgTAGACACTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAAc                                                                                                                                     <  2:329893/141‑1 (MQ=255)
             cccTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAgg                                                                                                                             >  2:95873/1‑141 (MQ=255)
                     gtggCGAAGGTATTAATTATGTTACAAAACGTTTAAATCATAGTGGATCAATTGCAGTTATTGGTATGACTGCCGGCAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGCGTAAACATTTTAGGAATTGACt                                                                                                                     >  2:46178/1‑141 (MQ=255)
                       ggCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCg                                                                                                                   >  1:222734/1‑141 (MQ=255)
                         cGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGt                                                                                                                 >  1:263541/1‑141 (MQ=255)
                         cGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGt                                                                                                                 >  1:421094/1‑141 (MQ=255)
                                        tGTTACAAAACGTTTAAATCATAGTGGATCAATTGCAGTTATTGGTATGACTGCCGGCAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGCGTAAACATTTTAGGAATTGACTCAGTATTTACTGCTATGaa                                                                                                  <  1:118908/141‑1 (MQ=255)
                                         gTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGaaa                                                                                                 >  2:450429/1‑141 (MQ=255)
                                         gTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGaaa                                                                                                 >  1:344975/1‑141 (MQ=255)
                                           tACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAAtt                                                                                               >  2:137348/1‑141 (MQ=255)
                                            aCAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTa                                                                                              >  2:130067/1‑141 (MQ=255)
                                              aaaGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAg                                                                                            >  2:297928/1‑141 (MQ=255)
                                                   ttttcactCATAGTGACTCAATTGCAGTTCTTGGTATGACTGCCAGTAATACTTATACTAATTCTGTATGCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAAGTAAGACAgc                                                                                       <  2:70707/134‑1 (MQ=255)
                                                     ttAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAgcgc                                                                                     <  2:93814/141‑1 (MQ=255)
                                                         aTCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGttt                                                                                 >  2:68446/1‑141 (MQ=255)
                                                           cATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTgg                                                                               <  1:162055/141‑1 (MQ=255)
                                                                tGGGTAAATAGCAGTTATTGGTATGACTGCCGGTAATACTTATAATAATTCTGTATACCCTAAAATTTTAAGAGGTGTAAACATTTTAGGAATTGACACGGTATTTACTGCTATAAAATTAAGACAGCGCGTTTggaaaaa                                                                          <  2:178007/141‑6 (MQ=255)
                                                                     cAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCg                                                                     <  1:110474/141‑1 (MQ=255)
                                                                     cAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTCTAAGAGGTTGAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCg                                                                     <  1:200543/141‑1 (MQ=255)
                                                                      aaTTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGc                                                                    <  2:24751/141‑1 (MQ=255)
                                                                        ttGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCaa                                                                  <  2:418908/141‑1 (MQ=255)
                                                                        ttGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCaa                                                                  <  1:143311/141‑1 (MQ=255)
                                                                        ttGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCaa                                                                  <  2:253240/141‑1 (MQ=255)
                                                                           cAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTAGTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAg                                                               <  2:427576/141‑1 (MQ=255)
                                                                                aTTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTa                                                          >  1:310237/1‑141 (MQ=255)
                                                                                    gTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATCTTAGGAATTGACTCGGTATTTACTGCTATGACATTAAGACAGCGCGTTTGGCGTCGCCTCGCAAAAGATTTAGTGc                                                      >  2:342672/1‑141 (MQ=255)
                                                                                      aTGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCt                                                    >  1:274901/1‑141 (MQ=255)
                                                                                           tGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGaaaa                                               <  1:74118/141‑1 (MQ=255)
                                                                                             ccGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAAtt                                             >  2:36218/1‑141 (MQ=255)
                                                                                              cGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAAttt                                            >  1:451420/1‑141 (MQ=255)
                                                                                                              aactCTGTCTTCCCTCACATGTTAAGAGGTGTAAACATGTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATcaagcaa                            <  1:384692/138‑1 (MQ=255)
                                                                                                                ttCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGGTTTAATGCCTGAAAATTTACATGAGATCAAGCAAGt                          >  1:131797/1‑141 (MQ=255)
                                                                                                                  cTGTATTCCCTCACATTTTAAGAGGCGTAAACATTTTAGGAATTGACTCAGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAAAGCCTGAAAATTTACATGAGATCAAGCAAGtta                        <  1:59960/141‑1 (MQ=255)
                                                                                                                    gTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGttatt                      >  2:288681/1‑141 (MQ=255)
                                                                                                                    gTATTCCCTCACATTTTAAGAGGCGTAAACATTTTAGGAATTGACTCAGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAAAGCCTGAAAATTTACATGAGATCAAGCAAGttatt                      >  2:196514/1‑141 (MQ=255)
                                                                                                                      aTTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTAc                    >  1:156967/1‑141 (MQ=255)
                                                                                                                        tCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACAt                  <  1:450429/141‑1 (MQ=255)
                                                                                                                         cccTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATAAAGCAAGTTATTACAtt                 >  1:131900/1‑141 (MQ=255)
                                                                                                                           cTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTtg               <  1:243905/141‑1 (MQ=255)
                                                                                                                                   ttAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACtt       <  1:172331/141‑1 (MQ=255)
                                                                                                                                    tAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAACTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAAGTTc      <  2:390229/141‑1 (MQ=255)
                                                                                                                                     aaGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATcaag                               >  1:77253/1‑115 (MQ=255)
                                                                                                                                     aaGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATcaag                               <  2:77253/115‑1 (MQ=255)
                                                                                                                                     aaGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTcc     >  2:400036/1‑141 (MQ=255)
                                                                                                                                     aaGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTcc     >  2:410484/1‑141 (MQ=255)
                                                                                                                                      agagGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTCCa    >  2:323354/1‑141 (MQ=255)
                                                                                                                                      agagGTGTAAACATTTTAGGAATTGACTAGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTCCa    >  2:111112/1‑141 (MQ=255)
                                                                                                                                        agGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTCCAGa  <  1:165420/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
AAGCTTGTGTAGACCCTGTTGGTGGCGAAGGTATTAATTATGTTACAAAGCGTTTAAATCATAGTGGGTCAATTGCAGTTATTGGTATGACTGCCGGTAATACTTATACTAATTCTGTATTCCCTCACATTTTAAGAGGTGTAAACATTTTAGGAATTGACTCGGTATTTACTGCTATGAAATTAAGACAGCGCGTTTGGCGTCGTCTCGCAAAAGATTTAATGCCTGAAAATTTACATGAGATCAAGCAAGTTATTACATTTGATGAACTTCCAGA  >  CP000730/2492111‑2492387

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