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

TTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGAC  >  CP000730/1513795‑1514055
                                                                                                                                       |                                                                                                                             
tttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATAGGTTTTATTGt                                                                                                                          <  2:275371/141‑1 (MQ=255)
    tttgtaatgatCGCTTTATTAGTAGCCAGTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAttt                                                                                                                      <  2:189392/132‑2 (MQ=255)
       gtatgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGa                                                                                                                   >  1:248903/1‑141 (MQ=255)
                                 tGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCAAATCGCGTAAAAAAACCAGCATCTATATCACTAATAATTTGTTTTATTTTATTAGGATTATTAATGCCATTATTTGGATTAta                                                                                         >  2:7116/1‑140 (MQ=255)
                                            ttGTCTCTTGGAAAAAAAAAAAAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTg                                                                              <  2:255690/141‑1 (MQ=255)
                                                      gATAAAGCATACAAAAAATACTGCAATAAAAGTACCACCTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGAGTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAt                                                                    <  1:243065/141‑1 (MQ=255)
                                                         aaaGCGTACAAAAAATACTGCAATAAAAGTCCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGa                                                                 <  2:180243/141‑1 (MQ=255)
                                                                          cTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATAGAATCAATATTATATATTa                                                <  1:156733/141‑1 (MQ=255)
                                                                            gCAATAAACGTATCACATCGAGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTTTATTTATAATTGAATTAATATTATATATTaaa                                              <  1:213931/141‑1 (MQ=255)
                                                                               aaaaaaTTACTACATCGCGTAAAAAAACCAGCATCTATATCACTCGTAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAt                                           <  2:129292/139‑1 (MQ=255)
                                                                                    aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgtcg                                      >  1:256652/1‑141 (MQ=255)
                                                                                     gTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATcgtcgt                                     >  2:114663/1‑141 (MQ=255)
                                                                                       aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGc                                   >  1:88753/1‑141 (MQ=255)
                                                                                         cacaTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTa                                 <  1:50161/141‑1 (MQ=255)
                                                                                           caTCGCGTAAAAAAACCACCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATCATTCCATCAAGATTATATATTAAAGATCGTCGTGCTaaa                               >  2:58439/1‑141 (MQ=255)
                                                                                                   aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATg                       >  2:81916/1‑141 (MQ=255)
                                                                                                                        aCTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGAc  >  1:107057/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                             
TTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGAC  >  CP000730/1513795‑1514055

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: