Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,028,049 A→G V199A (GTT→GCT)  USA300HOU_RS10230 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,028,0490AG100.0% 96.2 / NA 27V199A (GTT→GCT) USA300HOU_RS10230hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/11);  total (16/11)

GCATTGAACCAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAACTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTA  >  USA300TCH1516_ALE/2027913‑2028183
                                                                                                                                        |                                                                                                                                      
gCATTGAACCAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCttt                                                                                                                                    >  2:162060/1‑141 (MQ=255)
gCATTGAACCAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCtgt                                                                                                                                    >  2:79739/1‑139 (MQ=255)
  aTTGAACCAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTtat                                                                                                                                  <  2:108670/141‑1 (MQ=255)
      aaCCAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGa                                                                                                                              <  1:134788/141‑1 (MQ=255)
         cAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGAttt                                                                                                                           <  2:30138/141‑1 (MQ=255)
         cAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGAttt                                                                                                                           <  2:162794/141‑1 (MQ=255)
                           ggTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACt                                                                                                         >  2:128990/1‑141 (MQ=255)
                           ggTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACt                                                                                                         >  1:219281/1‑141 (MQ=255)
                                     tAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATaaaaa                                                                                               >  1:14092/1‑141 (MQ=255)
                                       aaTCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAagc                                                                                             >  1:176648/1‑141 (MQ=255)
                                         tCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAagcag                                                                                           >  1:65105/1‑141 (MQ=255)
                                            ccTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTa                                                                                        >  2:58758/1‑141 (MQ=255)
                                                    tctcCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTaaa                                                                                <  2:184150/141‑1 (MQ=255)
                                                               cGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTa                                                                     >  1:29286/1‑141 (MQ=255)
                                                               cGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTa                                                                     >  2:162517/1‑141 (MQ=255)
                                                                                 aaCTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCCGCTTTATATGATTTTTCATTCACGTAATTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAAc                                                   <  2:179077/141‑1 (MQ=255)
                                                                                 aaCTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAAc                                                   <  1:111025/141‑1 (MQ=255)
                                                                                         aGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACtt                                           <  1:180980/141‑1 (MQ=255)
                                                                                               ccGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTa                                     <  2:65105/141‑1 (MQ=255)
                                                                                                                 cAATATGAATATTGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAAcaacaa                   >  1:26371/1‑141 (MQ=255)
                                                                                                                            ttGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTAttt        <  1:162517/141‑1 (MQ=255)
                                                                                                                            ttGTTATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTAttt        <  2:29286/141‑1 (MQ=255)
                                                                                                                                tATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCt    >  1:97447/1‑141 (MQ=255)
                                                                                                                                tATTTTCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCt    >  1:118923/1‑141 (MQ=255)
                                                                                                                                  ttttCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTa  >  1:102749/1‑141 (MQ=255)
                                                                                                                                  ttttCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTa  >  1:222367/1‑141 (MQ=255)
                                                                                                                                  gtttCAGCTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTa  >  1:213469/2‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GCATTGAACCAATTGTTCTAATGTAGTGGTTCATCGATAAATCACCTCTACGTCTCCGATAAACGTTGATACGTAGATATTAACTGTATAGTTATCCGGTTTCATCATTTCTTCAATATGAATATTGTTATTTTCAACTTTATATGATTTTTCATTCACGTAAGTACTTCCATAAAAAGCAGCTACATGTAAATTAATATTGTAATTAACCGGCAATATAACCTGCACTTTACCTAAAATGTGTCTAACAACAATGGTATTATTTTCCTTA  >  USA300TCH1516_ALE/2027913‑2028183

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