Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,289,324 G→T 100% Q56K (CAG→AAG)  USA300HOU_2153 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,289,3240GT100.0% 75.6 / NA 26Q56K (CAG→AAG) USA300HOU_2153hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/13);  total (13/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAATTTTTA  >  CP000730/2289192‑2289452
                                                                                                                                    |                                                                                                                                
gaCCGCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCaa                                                                                                                          <  2:268992/140‑1 (MQ=255)
                ttCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTACCTTCTGACCAATAACAATTACTAatat                                                                                                          >  1:86122/1‑141 (MQ=255)
                         aCGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAg                                                                                                 <  2:139503/141‑1 (MQ=255)
                         aCGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAg                                                                                                 >  1:266628/1‑141 (MQ=255)
                         aCGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAg                                                                                                 <  2:198175/141‑1 (MQ=255)
                          agTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATTTTCTGACCACTAACAATTATTAATATCCCTTTAAga                                                                                                >  2:93399/2‑140 (MQ=255)
                           gTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTa                                                                                               <  2:231546/141‑1 (MQ=255)
                             cTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAAt                                                                                             >  1:84384/1‑141 (MQ=255)
                             cTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTGGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAAt                                                                                             >  2:136192/1‑141 (MQ=255)
                              tAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAAtt                                                                                            >  2:216114/1‑141 (MQ=255)
                                              tatTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTcc                                                                            >  1:107221/1‑141 (MQ=255)
                                                               aTACCCCCCATTGAATAACGATTTCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGtgatga                                                           <  1:230532/141‑1 (MQ=255)
                                                               aTACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGtgatca                                                           >  1:249852/1‑139 (MQ=255)
                                                                    cctcATTGAATTAAGATATCGAATAATTTAGTTGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCTATAACAATTACTTGTATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGagaaga                                                      <  1:129741/138‑1 (MQ=255)
                                                                    ccccATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGagaaga                                                      <  1:202181/141‑1 (MQ=255)
                                                                       cATTGAATAACTATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATACCAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATaa                                                   <  2:10782/141‑1 (MQ=255)
                                                                            aaTAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAATTaa                                              >  1:247098/1‑140 (MQ=255)
                                                                                    tatCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCaaa                                      <  1:216114/141‑1 (MQ=255)
                                                                                        gAATAATTTAGATGTTGCAGTCAAGATCAAGAACATAATTATCTTCTGACCAAGAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTACCTTAAAAACCGTGATGAGAAGATAAAGAACATAACaaaaaaa                                  >  2:164834/1‑137 (MQ=255)
                                                                                          ataatTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTTACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCa                                <  1:71873/141‑1 (MQ=255)
                                                                                           taatTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAg                               <  2:86122/141‑1 (MQ=255)
                                                                                                  gATGTTGCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAAAAATTACTAATATCACTTTAAGTAATTGTATTGATTTAGTTCCATTAAAGACCGTGATGAGAAGATAAAGTACATACAAAACTATCCGTAAAtaa                        >  2:222958/1‑139 (MQ=255)
                                                                                                        gCAGTCAAGTTCAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATa                  >  1:110518/1‑141 (MQ=255)
                                                                                                                 ttCAATATCATACTTATCTTCTAAACCATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCCGTAAATCAAGGGTACTCGTTAct         <  2:122241/141‑2 (MQ=255)
                                                                                                                   cAATATCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAAt       <  2:266628/141‑1 (MQ=255)
                                                                                                                        tCATACTTATCTTCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAATTTTTa  >  1:45713/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                
TACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAGATAAAGTACATACCAAACTATCAGTAAATCAAGGATACTCGTTACAATTTTTA  >  CP000730/2289192‑2289452

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