Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 18,890 C→G N182K (AAC→AAG)  yybT → DHH subfamily phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073018,8900CG85.7% 80.5 / 7.4 28N182K (AAC→AAG) yybTDHH subfamily phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base G (14/10);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TAATCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATT  >  CP000730/18755‑19029
                                                                                                                                       |                                                                                                                                           
taAGCAGTATCATTTCCACGTACGTTATTCTGAGAATGATCACTTCCTCTATATCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATACTAAAACAAACAGTGCCTCATTATTTTTAGATAAGTACGa                                                                                                                                        >  2:159302/1‑141 (MQ=255)
        aTCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTa                                                                                                                                >  1:74660/1‑141 (MQ=255)
           aTTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGc                                                                                                                             <  1:7546/141‑1 (MQ=255)
                  aGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAatat                                                                                                                      >  2:363770/1‑141 (MQ=255)
                             cTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGc                                                                                                           <  1:65069/141‑1 (MQ=255)
                                agaATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGc                                                                                                        >  1:295647/1‑141 (MQ=255)
                                agaATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGc                                                                                                        >  2:316217/1‑141 (MQ=255)
                                      aTCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCgg                                                                                                  >  1:389752/1‑141 (MQ=255)
                                          cTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAAt                                                                                              >  2:19146/1‑141 (MQ=255)
                                          cTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAAGGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAAt                                                                                              <  1:217317/141‑1 (MQ=255)
                                                  aTTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCaa                                                                                      <  2:151968/141‑1 (MQ=255)
                                                           atatAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAAcgc                                                                             <  1:384836/141‑1 (MQ=255)
                                                                aCTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTc                                                                        >  1:153971/1‑141 (MQ=255)
                                                                  tGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCAt                                                                      <  1:152873/141‑1 (MQ=255)
                                                                   gAACAAGTACAAACACATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCAtt                                                                     <  1:102434/141‑1 (MQ=255)
                                                                     aCAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAg                                                                   >  1:46289/1‑141 (MQ=255)
                                                                     aCAAGTACAAAAAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAg                                                                   <  1:252405/141‑1 (MQ=255)
                                                                                  aaTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACt                                                      >  2:350585/1‑141 (MQ=255)
                                                                                          atatTAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTataa                                              <  2:371012/141‑1 (MQ=255)
                                                                                          atatGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTataa                                              <  1:19146/141‑1 (MQ=255)
                                                                                                       aaaCCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCaaa                                 >  2:404247/1‑141 (MQ=255)
                                                                                                          ccAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGa                              >  1:1009/1‑141 (MQ=255)
                                                                                                                aTTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGt                        >  1:48009/1‑141 (MQ=255)
                                                                                                                aTTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGGATAATATATTTTTCAAAAGATACAGt                        >  2:207251/1‑141 (MQ=255)
                                                                                                                 ttGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGtt                       >  1:44164/1‑141 (MQ=255)
                                                                                                                              tttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGt          >  2:318763/1‑141 (MQ=255)
                                                                                                                                tAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAg        <  1:363770/141‑1 (MQ=255)
                                                                                                                                      aCTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCAGATCAATTCGTAGCCTAtt  >  1:265064/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                           
TAATCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATT  >  CP000730/18755‑19029

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