Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,064,126 G→C 100% Y121* (TAC→TAG USA300HOU_1000 ← LytR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,064,1260GC100.0% 68.1 / NA 21Y121* (TAC→TAGUSA300HOU_1000LytR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (11/10);  total (11/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTC  >  CP000730/1064005‑1064253
                                                                                                                         |                                                                                                                               
ccTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGaa                                                                                                              <  2:256032/141‑1 (MQ=255)
   cGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACTACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGt                                                                                                           <  2:87535/141‑1 (MQ=255)
     aCAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAg                                                                                                         >  2:7630/1‑141 (MQ=255)
     aCAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAg                                                                                                         >  2:178722/1‑141 (MQ=255)
          tGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGAt                                                                                                    >  1:187681/1‑141 (MQ=255)
                cATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATc                                                                                              >  1:282073/1‑141 (MQ=255)
                  tatCGACTCTCACATAATAAGCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCAc                                                                                            <  1:40432/141‑1 (MQ=255)
                             acaTAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTACTAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTa                                                                                 >  1:62956/1‑141 (MQ=255)
                                             gTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGAt                                                                 <  1:192917/141‑1 (MQ=255)
                                                     ctGGTAGCTTCTATCGAATCCATTGCGGCCATTGGACCACCATATGCATGTGCATGGGTTATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTttgat                                                         <  2:285841/140‑1 (MQ=255)
                                                       ggTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTttgattg                                                       >  1:120736/1‑141 (MQ=255)
                                                                    gAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATCCTAATTTGATGTTTTGATTGATTAAAAGTAGtt                                          <  2:265420/141‑1 (MQ=255)
                                                                           ttGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCa                                   >  1:169120/1‑141 (MQ=255)
                                                                                           accaTATGCATGTGCATGGGTAATCTTATCCTAATAGCGATCTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCGACCTAATTTTTTGTTTTGATTGATTGAAAGTAGCTAAATTCACAGCGTCTGATCTaga                   <  2:116459/141‑1 (MQ=255)
                                                                                               tatGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATcacgaatgc               >  2:213045/1‑133 (MQ=255)
                                                                                                atGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGtt              >  1:81472/1‑141 (MQ=255)
                                                                                                      gtgcatgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCAt        >  1:244193/1‑141 (MQ=255)
                                                                                                         catgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGGGTGTTCAGCATCCt     >  2:225338/1‑141 (MQ=255)
                                                                                                           tgGGTAATCTTATCCTAATAGCCAACTTTATGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTACATTGTTCAGCATCCTGt   <  1:225003/141‑1 (MQ=255)
                                                                                                           tgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGt   <  1:166330/141‑1 (MQ=255)
                                                                                                            gggTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTc  <  1:152436/141‑1 (MQ=255)
                                                                                                                         |                                                                                                                               
CCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTC  >  CP000730/1064005‑1064253

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