Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,513,930 A→T 100% L413F (TTA→TTT)  USA300HOU_1379 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,513,9300AT100.0% 116.5 / NA 35L413F (TTA→TTT) USA300HOU_1379hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/15);  total (20/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACA  >  CP000730/1513802‑1514065
                                                                                                                                |                                                                                                                                       
gtatgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAgg                                                                                                                              <  2:285299/140‑1 (MQ=255)
  atgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                                                                            <  1:546405/140‑1 (MQ=255)
  atgtacCGTTTTATTAGTAGCCATGGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAGTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                                                                            <  2:434615/140‑1 (MQ=255)
             tattaGTAGGCATAGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTCTATTGTATTAGGATTATTAATGCCa                                                                                                                 <  2:481059/140‑1 (MQ=255)
             tactatcagCAATTGGCATGGGGTTTGTCTCTGGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATGATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                                                                 <  2:505979/133‑1 (MQ=255)
                 agtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                                             <  1:615604/140‑1 (MQ=255)
                  gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCATATCTCGTAATAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                                                                            <  1:44381/140‑1 (MQ=255)
                                                aTAAAGCGTACAAAAAATACTGCAATGAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                                                                              >  2:451659/1‑140 (MQ=255)
                                                aTAAAGCGTACAAAAAATACGGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCCTCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                                                                              <  2:305966/140‑1 (MQ=255)
                                                 tAAAGCGGACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                                                                             >  2:720557/1‑140 (MQ=255)
                                                  aaaGCGTACAAAAAATACTTAAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATCGTCTTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                                                                            <  2:424274/140‑1 (MQ=255)
                                                    aGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGaa                                                                          <  1:388149/140‑1 (MQ=255)
                                                     gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                                                                         >  2:521469/1‑140 (MQ=255)
                                                     gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                                                                         >  2:587395/1‑140 (MQ=255)
                                                       gTACAAAAAATACTGCAATAAAAGGACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTAGATTATCACTTATCCTAGTATTTATAATTGAATTa                                                                       >  1:464064/1‑140 (MQ=255)
                                                           aaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                                                                   >  1:127446/1‑140 (MQ=255)
                                                             aaaaTACTGCAATAAAAGTATCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatta                                                                 <  1:529526/140‑1 (MQ=255)
                                                               aaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTTAATTAATATtata                                                               >  2:642547/1‑140 (MQ=255)
                                                                  agtGCAATAAAAGTAACAGATTCCGTAAAATAACCAGCCTCTATATCACTCATAATATGTGTTATTGTATTAGGATTATTAATGCCATGATTTGGATTATCACTTATCCTTGTATTTATAAGTGAATTAATATTatatat                                                            <  2:137028/138‑1 (MQ=255)
                                                                  aCTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTatatat                                                            <  1:574076/140‑1 (MQ=255)
                                                                         tAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt                                                     >  1:290235/1‑140 (MQ=255)
                                                                             agcaCCACATCGCGTAAAAAAACCAGCATCTATATCACAAAGACAATGTTTTATTGTAATAGGATTATTAATGCAATTATTTGCCTTATCACTCATCCTCGGATTTATAATTCACTTAATATTTTCTACTCAAGAccgta                                                 >  2:484522/4‑139 (MQ=255)
                                                                               tACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATcgtcgt                                               >  2:51115/1‑140 (MQ=255)
                                                                                aCCACATCGCGTAAAAAAACCAGCTTCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTg                                              <  1:415924/140‑1 (MQ=255)
                                                                                aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATATTTGAATTAATATTATATATTAAAGATTGTCGTg                                              >  1:488866/1‑140 (MQ=255)
                                                                                                   ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTa                           <  1:501498/140‑1 (MQ=255)
                                                                                                       cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGt                       >  1:726872/1‑140 (MQ=255)
                                                                                                       cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGt                       >  2:528285/1‑140 (MQ=255)
                                                                                                         tCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGttt                     >  2:45849/1‑140 (MQ=255)
                                                                                                           tataTCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAACTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTg                   >  2:521370/1‑140 (MQ=255)
                                                                                                                   tCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACg           >  2:175224/1‑140 (MQ=255)
                                                                                                                        atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATTTTATAGATCGTCGTGCTAAACAATATTGCACTTATAGTTTTGAACTGATGAGAtt      >  1:506396/1‑140 (MQ=255)
                                                                                                                          atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATGTGTCAGAATAAAGATCGTCGTGCTAAACAATAATGCATTTAAAGTTTTGAACTGACGTAATTTa    >  1:167803/1‑140 (MQ=255)
                                                                                                                           tGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGACCTTAAAGTTTAGAACTGCCGAAATTTAc   >  2:23929/1‑140 (MQ=255)
                                                                                                                            gTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAAATCGTCGTACTAAACAATAATGCACTTAAAGTTTTAAACTGACGAAATTTACa  >  1:202465/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                       
GTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACA  >  CP000730/1513802‑1514065

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