Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 19,009 C→T 100% S222F (TCC→TTC)  yybT → DHH subfamily phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073019,0090CT100.0% 165.2 / NA 49S222F (TCC→TTC) yybTDHH subfamily phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/27);  total (22/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAG  >  CP000730/18871‑19146
                                                                                                                                          |                                                                                                                                         
cGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCg                                                                                                                                         >  2:360560/1‑141 (MQ=255)
  aCATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGAt                                                                                                                                       <  1:86266/141‑1 (MQ=255)
      tattTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAAt                                                                                                                                   >  1:96996/1‑141 (MQ=255)
         ttttAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCg                                                                                                                                >  1:60739/1‑141 (MQ=255)
          attAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGt                                                                                                                               >  1:385087/2‑141 (MQ=255)
            tAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAg                                                                                                                             <  1:500189/141‑1 (MQ=255)
                 aaCTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTAt                                                                                                                        <  1:189900/141‑1 (MQ=255)
                            gaTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTTGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCaaaa                                                                                                             <  1:364421/141‑1 (MQ=255)
                                     aaaTATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGc                                                                                                    >  1:359009/1‑141 (MQ=255)
                                          tgaatgaTACGCAGCGTTCGGAAAGCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACt                                                                                               <  2:412209/141‑1 (MQ=255)
                                           gaatgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCTATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACtt                                                                                              >  2:272611/1‑141 (MQ=255)
                                           gaatgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACtt                                                                                              >  2:325350/1‑141 (MQ=255)
                                             atgaTACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGTCTATTTAAATCAAAAAATATTAGCTGACTTag                                                                                            >  2:65117/1‑141 (MQ=255)
                                                    gCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATc                                                                                     <  1:86206/141‑1 (MQ=255)
                                                     cAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATACGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCt                                                                                    <  1:412623/141‑1 (MQ=255)
                                                          ttCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAAtt                                                                               <  1:169124/141‑1 (MQ=255)
                                                            cGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGACAAATAAAAAAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTg                                                                             <  2:385087/141‑1 (MQ=255)
                                                                      cTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGAATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAg                                                                   >  2:259769/1‑141 (MQ=255)
                                                                       tCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGt                                                                  >  2:438861/1‑141 (MQ=255)
                                                                            ggTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAAtt                                                             >  2:426388/1‑141 (MQ=255)
                                                                             gTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTa                                                            <  1:328945/141‑1 (MQ=255)
                                                                               aaCGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACg                                                          >  1:258287/1‑141 (MQ=255)
                                                                               aaCGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACg                                                          <  2:58309/141‑1 (MQ=255)
                                                                                       tCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGaaaaaa                                                  <  2:359009/141‑1 (MQ=255)
                                                                                          ttAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAgtg                                               <  1:55154/141‑1 (MQ=255)
                                                                                           tAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAgtgt                                              >  2:171305/1‑141 (MQ=255)
                                                                                             gTCGATGGGCAACTGAGTATAATATATTTTTAAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGTTATCTTGAGTAAAAAAAAAAAAAAAAGTGTTg                                            <  2:287723/141‑1 (MQ=255)
                                                                                              tCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTgg                                           <  1:365064/141‑1 (MQ=255)
                                                                                               cGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGt                                          <  1:434352/141‑1 (MQ=255)
                                                                                               cGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCATATATAAATAAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGt                                          <  2:402635/141‑1 (MQ=255)
                                                                                                  tGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTAt                                       >  1:190620/1‑141 (MQ=255)
                                                                                                    ggCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCg                                     <  2:106622/141‑1 (MQ=255)
                                                                                                    ggCAACTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCg                                     <  1:280389/141‑1 (MQ=255)
                                                                                                      cAACTGAGTATAATATATTTTCTAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTGAAGCAAAAAATATTTGCTGACTTAGAAGAATCTAAATTTGATATCTTGTGTCAATTACGTGAAAAAAGTGTTGGTTTTCGTg                                   <  2:273768/141‑1 (MQ=255)
                                                                                                       aaCTGAGTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATGTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTgg                                  <  1:386214/141‑2 (MQ=255)
                                                                                                            agaataataTATTTTTCAAAAGATACAGTTTCGATAAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATGTGATATCTTGAGTCAAAAAAGTAAAAAAAGTGTTGGTTATCGTGCCCAAt                             <  2:451529/138‑1 (MQ=255)
                                                                                                             gTATAATATATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAAtt                            >  2:65570/1‑141 (MQ=255)
                                                                                                                taataTATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAAc                         >  1:384982/1‑141 (MQ=255)
                                                                                                                taataTATTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAAc                         >  1:310044/1‑141 (MQ=255)
                                                                                                                  atatatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACAt                       <  1:20548/141‑1 (MQ=255)
                                                                                                                    atatTTTTCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATATAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTa                     >  2:483725/1‑141 (MQ=255)
                                                                                                                          ttCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATc               >  1:493317/1‑141 (MQ=255)
                                                                                                                           tCAAAAGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCg              >  1:79221/1‑141 (MQ=255)
                                                                                                                             aaaaGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGgt            <  2:388446/141‑1 (MQ=255)
                                                                                                                             aaaaGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGgt            <  1:498317/141‑1 (MQ=255)
                                                                                                                              aaaGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTAAGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGgtg           <  2:96996/141‑1 (MQ=255)
                                                                                                                                aGATACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGtt         >  2:13505/1‑141 (MQ=255)
                                                                                                                                  aTACAGTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTgg       <  2:328497/141‑1 (MQ=255)
                                                                                                                                       gTTTCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAg  >  1:115839/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTAGCCTATTTAAATCAAAAAATATTAGCTGACTTAGAAGAATCTAAATTTGATATCTTGAGTCAATTACGTGAAAAAAGTGTTGGTTATCGTGCCCAATTAACATTAAGTATCGGTGTTGGTGAAG  >  CP000730/18871‑19146

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: ‑