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% 125.0 / NA 39E175* (GAG→TAG) USA300HOU_0930possible HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/17);  total (22/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTACCCA  >  CP000730/973087‑973350
                                                                                                                              |                                                                                                                                         
cGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCAt                                                                                                                             >  1:221890/1‑141 (MQ=255)
 gAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCATCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATa                                                                                                                            >  1:274445/1‑141 (MQ=255)
   aTGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGAGCGGCATaaa                                                                                                                          >  2:274823/1‑141 (MQ=255)
    tGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCGCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaaa                                                                                                                         >  2:136648/1‑141 (MQ=255)
      cAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaaaag                                                                                                                       >  1:23112/1‑140 (MQ=255)
      cAATAATATTATTTCGATCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAAt                                                                                                                       >  1:185557/1‑141 (MQ=255)
          aataTTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGt                                                                                                                   >  1:329141/1‑141 (MQ=255)
           atatTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGAAATGGTGCGCCCCAGAGTCGACACCAAAGATGATCGGCATAAAAATGAGTa                                                                                                                  >  2:217433/1‑141 (MQ=255)
                       tCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATCATTTTg                                                                                                      >  1:155805/1‑141 (MQ=255)
                       tCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTg                                                                                                      >  2:309133/1‑141 (MQ=255)
                           tatTTACTAGTTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAttt                                                                                                  <  1:240609/141‑1 (MQ=255)
                               tAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAgg                                                                                              <  1:305393/141‑1 (MQ=255)
                                aaaaaTTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATTTGATCGGCATCAAAATGAGTAAGCATAATTTTGATTTCAGGt                                                                                             >  2:75754/1‑141 (MQ=255)
                                aaaaaTTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGt                                                                                             >  1:13526/1‑141 (MQ=255)
                                aaaaaTTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTAAGGt                                                                                             >  2:44604/1‑141 (MQ=255)
                                        tCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACt                                                                                     >  1:131390/1‑141 (MQ=255)
                                         cTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACtt                                                                                    <  1:82446/141‑1 (MQ=255)
                                         cTAACTTGCTCAATGCCTCTTGCTTTAGTAATACCAAGTGTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACtt                                                                                    <  2:157024/141‑1 (MQ=255)
                                             cTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttctt                                                                                <  1:58701/141‑1 (MQ=255)
                                                gCTCAATGCAGCTTGCTTTATTAATACCACGTCTTTAGTGACCAAAGACAGGGGATTTGGCGCCCCAGCGTCGCTGCTAAATATGATCGGCATAAAAGTGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCgg                                                                             >  1:244276/1‑141 (MQ=255)
                                                  tCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                                                                           <  2:247321/141‑1 (MQ=255)
                                                    aaTGCCTCTTGCTTTATTAATACCAAGATTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTATACTTTCTTCGGCTTc                                                                         >  1:317856/1‑141 (MQ=255)
                                                      tGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCaa                                                                       <  1:88145/141‑1 (MQ=255)
                                                      tGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCaa                                                                       <  1:217433/141‑1 (MQ=255)
                                                       gCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                                                                      <  1:220605/141‑1 (MQ=255)
                                                        ccgctTGCGTTCTTAATCCCAATTGTAACAATTTCAATGACACGGAATGGTGCGCCCCAGCGTCGACGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                                                                     <  1:91415/138‑1 (MQ=255)
                                                        ccTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                                                                     <  2:248679/141‑1 (MQ=255)
                                                            ttGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaaa                                                                 >  2:169857/1‑141 (MQ=255)
                                                             tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAt                                                                <  2:329141/141‑1 (MQ=255)
                                                                       aTACCTAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCAGCATAAAAATGAGTAAGCATATTTTTCATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAggg                                                      <  2:286322/141‑1 (MQ=255)
                                                                         aCCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGa                                                    >  1:219611/1‑141 (MQ=255)
                                                                          ccAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGAATATCTGTGTTTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGTTAGGGGAt                                                   >  1:63068/1‑141 (MQ=255)
                                                                              attttACAATTTCAAAGACAAGGAATTGTGCGCCCCAGCTTCTATGCTAAAAAAAAAATAAAAAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCtt                                               <  2:219288/140‑1 (MQ=255)
                                                                                  tACAATTGCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCaa                                           <  2:116483/141‑1 (MQ=255)
                                                                                              gACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTaa                               >  2:96844/1‑141 (MQ=255)
                                                                                                                ccAGCGTCATTGCTAAATTTGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCt             <  2:302792/141‑1 (MQ=255)
                                                                                                                   gCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAGTTACCAGTTTGAATTCTTgg          >  1:329075/1‑141 (MQ=255)
                                                                                                                        gATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTAc     >  1:74042/1‑141 (MQ=255)
                                                                                                                           gCTAAATAAGATAGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTACCCa  <  2:19582/141‑1 (MQ=255)
                                                                                                                              |                                                                                                                                         
CGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTACCCA  >  CP000730/973087‑973350

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