Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,416,458 A→T 100% I216N (ATT→AAT)  USA300HOU_1314 ← oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,416,4580AT97.7% 45.8 / ‑6.3 18I216N (ATT→AAT) USA300HOU_1314oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (9/8);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATC  >  CP000730/1416323‑1416576
                                                                                                                                       |                                                                                                                      
ttAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTAc                                                                                                                   >  2:349543/1‑141 (MQ=255)
ttAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTGGACCTTCTTTTAAAACATAATTggat                                                                                                                   >  1:237186/1‑138 (MQ=255)
                aaTCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATa                                                                                                   >  2:314699/1‑141 (MQ=255)
                      tctTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGAc                                                                                             <  1:326446/141‑1 (MQ=255)
                           aaTTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAg                                                                                        >  1:284047/1‑141 (MQ=255)
                                   tcgaggtttGATTGGTCGTTGAGTATACGGGTGCAACTGTTCATACCTGAAATGTTTAATATTACGTCGTTCAATGATTTGACCTTCTTTTAACACATAATTGTACTGACAATATTTCAATACATGGCTTAAGTGAtgtgt                                                                                <  2:302362/135‑1 (MQ=255)
                                       tatTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGata                                                                            <  1:154365/141‑1 (MQ=255)
                                                         gTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                                                          <  2:102629/141‑1 (MQ=255)
                                                          tataCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCAt                                                         <  2:65672/141‑1 (MQ=255)
                                                           ataCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATATGTGATAATAAATAATGTTTGACCAtt                                                        >  2:325125/1‑140 (MQ=255)
                                                                 aTGCAAATGCTCATACTTTAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTaa                                                  <  1:314699/141‑1 (MQ=255)
                                                                  tGCAAATGCTCATACCTGAAATGATTAATATTCCCTCGTCCAATAATTTGACCTTCTTTTAAAACATAATTGTACTGACAATAGTTCCATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAt                                                 <  1:38441/141‑1 (MQ=255)
                                                                                     aaTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAAt                              >  1:306965/1‑141 (MQ=255)
                                                                                        gATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTCATACAATATGCTGGAATAAATCCa                           >  1:232633/1‑141 (MQ=255)
                                                                                            aaTATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGTTAATAAAGAAGGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCat                       >  1:373891/1‑140 (MQ=255)
                                                                                               aatACATCGTTCAATGATTTGACCTTCTTCAAAAACATAAAAGAACAAAAAATAAATCAACACATTACTTAAGTTACGCTTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTg                    <  2:181562/139‑1 (MQ=255)
                                                                                                         tCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCACGTTTTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTc          >  2:27520/1‑141 (MQ=255)
                                                                                                                 tgtacatgtttttAAAACATAATTGTACTGAAAATATTTCAATACATGACTTAATTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATc  <  2:348263/132‑1 (MQ=255)
                                                                                                                                       |                                                                                                                      
TTAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATC  >  CP000730/1416323‑1416576

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