Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,856,087 T→A 100% P80P (CCT→CCA USA300HOU_2698 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,856,0870TA100.0% 129.7 / NA 40P80P (CCT→CCAUSA300HOU_2698hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (19/21);  total (19/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTA  >  CP000730/2855951‑2856219
                                                                                                                                        |                                                                                                                                    
gggtGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAAc                                                                                                                                  <  2:89016/138‑1 (MQ=255)
 gaTGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACg                                                                                                                                 >  2:91639/1‑141 (MQ=255)
  gcccATGGCGTTCGTCCTCACGCCATCGCAGTATTGGAACAAGCTTGATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGt                                                                                                                                <  2:66676/138‑1 (MQ=255)
       tGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTcaaca                                                                                                                           >  2:301824/1‑141 (MQ=255)
            ttGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGa                                                                                                                      <  1:310422/141‑1 (MQ=255)
              ggTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATa                                                                                                                    <  1:12185/141‑1 (MQ=255)
                 ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                 <  1:342666/141‑1 (MQ=255)
                           tatCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttt                                                                                                       >  2:193833/1‑140 (MQ=255)
                            tctaaTATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGCAAGTTCAGTCAGCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAgttgt                                                                                                      <  2:84255/138‑1 (MQ=255)
                              cAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTa                                                                                                    >  1:272056/1‑141 (MQ=255)
                                  atTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGACAGGATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATACTTGTTAcagc                                                                                                <  1:250562/141‑1 (MQ=255)
                                     ggAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAGCAAATGATAGATTTTATAGTTGTTACagcagc                                                                                             >  1:144320/1‑141 (MQ=255)
                                        aCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAtt                                                                                          >  1:301808/1‑141 (MQ=255)
                                                  ttacAGCACCAACTACATTTTAAAACGGACTTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTAGTTACAGCAGCATTAGATATTGaa                                                                                >  2:3319/4‑141 (MQ=255)
                                                        cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACCGCAGCATTAGATATTGActtattc                                                                          >  1:282980/1‑140 (MQ=255)
                                                        cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAAttatta                                                                          >  1:108375/1‑141 (MQ=255)
                                                         aCCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAGTGATAGATTTTATATTTGTTACAGCAGCATTAGATATTGAATGATTAc                                                                         >  2:160570/1‑141 (MQ=255)
                                                         aCCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTCTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTAc                                                                         <  1:91639/141‑1 (MQ=255)
                                                            aaCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGtt                                                                      >  1:272358/1‑141 (MQ=255)
                                                            aaCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGtt                                                                      <  2:110833/141‑1 (MQ=255)
                                                             aCTACATTTTAACACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGCGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTc                                                                     >  2:245619/1‑141 (MQ=255)
                                                               tACATTTTAAACCGTAAGTTCAGTCATCATTCCCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                   <  1:160156/141‑1 (MQ=255)
                                                               tACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCAc                                                                   <  1:34370/141‑1 (MQ=255)
                                                                aCATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTCAATTATTACGTTCACt                                                                  >  2:13857/1‑141 (MQ=255)
                                                                  aTTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCGATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGTATTATTACGTTccgtgg                                                                >  2:230883/1‑136 (MQ=255)
                                                                       aaacctgAAGTTCAGTCATCATTACCTTTTAAGTTAACATGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAAGGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTaca                                                           <  1:308626/135‑1 (MQ=255)
                                                                             gAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGaa                                                     <  1:332281/141‑1 (MQ=255)
                                                                                 ttCAGTCATCATTACCTTTTAAGTAAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAattta                                                 >  1:220535/1‑141 (MQ=255)
                                                                                  tCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttat                                                <  1:301786/141‑1 (MQ=255)
                                                                                     gTCATCATTACCTTTTCAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTg                                             >  1:204504/1‑141 (MQ=255)
                                                                                      tcatcaTTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAACCGTCAACAAATGATAGATTTTATAGTGGTTACAGAAGGCTTAGATCTTGAATTATTACGTTCACTGGGTAAAACGGATTGTATTTgg                                            >  1:43907/1‑141 (MQ=255)
                                                                                      gcatcatTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTACAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTgg                                            <  1:242509/140‑1 (MQ=255)
                                                                                                     aaGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGc                             <  2:61521/141‑1 (MQ=255)
                                                                                                                  aTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACg                >  1:105981/1‑141 (MQ=255)
                                                                                                                         aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCa         <  2:323538/141‑1 (MQ=255)
                                                                                                                              aTCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAatat    <  1:158579/141‑1 (MQ=255)
                                                                                                                               ccATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATAtt   <  1:301824/140‑1 (MQ=255)
                                                                                                                                cATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTa  >  1:4261/1‑141 (MQ=255)
                                                                                                                                cATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTa  <  1:297718/141‑1 (MQ=255)
                                                                                                                                cATGCGCAAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTa  <  2:129950/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AGATGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTA  >  CP000730/2855951‑2856219

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