Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 973,213 C→A 100% E175* (GAG→TAG)  USA300HOU_0930 ← possible HAD superfamily hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730973,2130CA100.0% 79.4 / NA 25E175* (GAG→TAG) USA300HOU_0930possible HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/10);  total (15/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCTTCATCACCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATT  >  CP000730/973076‑973345
                                                                                                                                         |                                                                                                                                    
aTCTTCATCACCGACTGCCTTCATTTTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTCATACCAAGTTTTACAATTTCAATGACAGGAAATGGTGCGCCCCAGCGTCGATGCTAAat                                                                                                                                   >  2:30544/1‑141 (MQ=255)
        caCCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCgg                                                                                                                           <  1:135380/141‑1 (MQ=255)
           cGAAGGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGGTAAATATGATCGGCAt                                                                                                                        <  1:134607/141‑1 (MQ=255)
                                        tttAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCa                                                                                           <  1:22126/141‑1 (MQ=255)
                                                            cTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGc                                                                       >  1:36845/1‑141 (MQ=255)
                                                                aTGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCa                                                                   <  1:115537/141‑1 (MQ=255)
                                                                         gCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGGTTCAATTAAAAtt                                                          >  2:116467/1‑141 (MQ=255)
                                                                            ttattaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAg                                                       >  2:234572/1‑141 (MQ=255)
                                                                             tattaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGAGTTGATGCTAAATATGATCGCCATAAAAATGAGTAATCATCTTTTTGATTTCAGGTCTTTTACTTTCTTCGGCTTCAATTAAAATTGacg                                                      >  2:132766/1‑139 (MQ=255)
                                                                             tatgaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCTCCGATGCTAAATATGATCGGCATAAAAATTAGTAAGCATATTTTTATTATCAGGTATTTTACTTTCTTCGTTTTCAATTACAATTGCgg                                                      >  2:132776/1‑141 (MQ=255)
                                                                               ttaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTa                                                    >  1:67153/1‑141 (MQ=255)
                                                                                 aagaCCAAGTTTTACAAGTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCACATCTTTGATTTCAGGTATCTCACTTTCTTCGACTTCAACTAAAATTGAgggtgc                                                  <  2:225864/138‑5 (MQ=255)
                                                                                      cAAGTCTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGAtt                                             <  1:237834/141‑1 (MQ=255)
                                                                                         gTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCtt                                          >  2:249798/1‑141 (MQ=255)
                                                                                           tttACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTc                                        >  2:79137/1‑141 (MQ=255)
                                                                                              aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAg                                     >  2:78155/1‑141 (MQ=255)
                                                                                                aaTTGCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAACGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTg                                   <  1:39982/141‑1 (MQ=255)
                                                                                                           cAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAAGGGATTCTTTCAAGTGGACAAGTAAAt                        >  1:192683/1‑141 (MQ=255)
                                                                                                                  tGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGt                 <  2:36845/141‑1 (MQ=255)
                                                                                                                     tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGCTTTCCGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAACTGACCAAGTAAATTACCAGTTTg              >  2:235357/1‑141 (MQ=255)
                                                                                                                     tGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTg              <  2:67153/141‑1 (MQ=255)
                                                                                                                             aGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTg      <  1:79137/141‑1 (MQ=255)
                                                                                                                             aGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTGGAATTCTTg      >  2:107761/1‑141 (MQ=255)
                                                                                                                                gTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGAt   >  1:68034/1‑141 (MQ=255)
                                                                                                                                 tCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTAATTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGAtt  >  1:144279/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                    
ATCTTCATCACCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATT  >  CP000730/973076‑973345

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