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% 96.7 / NA 34L413F (TTA→TTT) USA300HOU_1379hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (19/15);  total (19/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACAAA  >  CP000730/1513795‑1514067
                                                                                                                                       |                                                                                                                                         
tttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATACAAGTACCACATCGCGTAAAAAACCCAGCATCTATATCACTCATAATATGTTTTATTGt                                                                                                                                      >  2:90801/1‑141 (MQ=255)
  tatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAt                                                                                                                                    <  1:204995/141‑1 (MQ=255)
      ttaatgtagCGTTTTATTAGTAGTCATTTGCATGGGGTTTGTCTCTGGGATAAAGCGTACAAAAAATACTGCAATAAAAGTTCCACATCTCGAATAAAAACCTGCTTCTATATCACTCATAATATGTTTTTTTTTATTAgg                                                                                                                                <  2:337538/138‑1 (MQ=255)
         atgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTAGATAGAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatt                                                                                                                             >  1:392231/1‑141 (MQ=255)
              tCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTaa                                                                                                                        >  1:93208/1‑141 (MQ=255)
                   ttattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCa                                                                                                                   >  2:483876/1‑141 (MQ=255)
                       tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                                               >  1:256633/1‑141 (MQ=255)
                        agttgCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTAAAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                                                                              <  2:291871/141‑1 (MQ=255)
                        agtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattatt                                                                                                              <  2:510016/141‑1 (MQ=255)
                           agCCAGTGTCATGAGGTTTGTCGTTTGAATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCTTCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTgg                                                                                                           <  2:328967/141‑1 (MQ=255)
                              cATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAtt                                                                                                        <  1:147/141‑1 (MQ=255)
                                  ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                                                                    <  2:98459/141‑1 (MQ=255)
                                  ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCAATTATTTGGAtgttca                                                                                                    >  2:303462/1‑136 (MQ=255)
                                      tGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCCTTATTTGGATTATCACTTa                                                                                                >  1:139041/1‑141 (MQ=255)
                                          gtctgttgctTGGATAAAGCGTACAAAAAATACTGCAATAACAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACGCATTATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCt                                                                                            <  2:298476/133‑1 (MQ=255)
                                            ttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTg                                                                                          >  1:486601/1‑141 (MQ=255)
                                              gTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCCCTTATCCTTGTa                                                                                        >  1:68618/1‑141 (MQ=255)
                                                ctctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAtt                                                                                      >  1:211376/1‑141 (MQ=255)
                                                       aTAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                                                                               <  1:331643/141‑1 (MQ=255)
                                                            gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGT‑AAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTa                                                                         >  2:534980/1‑141 (MQ=255)
                                                                     aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATATTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtata                                                                 >  1:457217/1‑141 (MQ=255)
                                                                           tGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATCTTATATATTaa                                                           >  2:302144/1‑141 (MQ=255)
                                                                             cAATAAAAGTACCAAAACGCGTAAAAAAACCAACATCTATATCACTCATAATATGTTGTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAg                                                         <  1:100924/141‑1 (MQ=255)
                                                                                   aaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgtc                                                   <  2:486601/141‑1 (MQ=255)
                                                                                       aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTgt                                               >  1:109602/1‑140 (MQ=255)
                                                                                        ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATAATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCt                                              >  2:204393/1‑141 (MQ=255)
                                                                                             tCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTGAAGATCGTCGTGCTAAACa                                         >  2:364515/1‑141 (MQ=255)
                                                                                                   aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATg                                   >  2:430883/1‑141 (MQ=255)
                                                                                                       aaaCCAGCCTCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt                               <  1:67256/141‑1 (MQ=255)
                                                                                                       aaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt                               <  1:546682/141‑1 (MQ=255)
                                                                                                             gCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGAATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCAATTAAAGt                         >  2:413028/1‑141 (MQ=255)
                                                                                                                 cTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTg                     <  2:449627/141‑1 (MQ=255)
                                                                                                                              aataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAtt        >  2:90476/1‑141 (MQ=255)
                                                                                                                              aataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAtt        <  2:434577/141‑1 (MQ=255)
                                                                                                                                    ttttATTGTTTTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCCAAACCATGATGCACTTAAAGTTTTGAACTGACGAAATTTACaaa  <  2:429786/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
TTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACAAA  >  CP000730/1513795‑1514067

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