Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 134,393 G→A L12L (TTG→TTA sbnA → putative siderophore biosynthesis protein SbnA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE134,3930GA96.7% 102.8 / NA 30L12L (TTG→TTAsbnAputative siderophore biosynthesis protein SbnA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/19);  minor base C (1/0);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTGTTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCAAGTAC  >  USA300TCH1516_ALE/134257‑134516
                                                                                                                                        |                                                                                                                           
aaTGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAg                                                                                                                         <  1:382348/141‑1 (MQ=255)
aaTGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAg                                                                                                                         <  1:309188/141‑1 (MQ=255)
  tGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGAt                                                                                                                       <  1:348586/141‑1 (MQ=255)
    cATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTc                                                                                                                     <  1:406284/141‑1 (MQ=255)
      tAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTAATTATTAGATTCTg                                                                                                                   >  1:308107/1‑141 (MQ=255)
                     tAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTa                                                                                                    <  1:127360/141‑1 (MQ=255)
                               aTAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACt                                                                                          <  1:128825/141‑1 (MQ=255)
                                    gagaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatc                                                                                     <  1:434736/141‑1 (MQ=255)
                                    gagaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatc                                                                                     <  1:230252/141‑1 (MQ=255)
                                     agaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatca                                                                                    <  1:413277/141‑1 (MQ=255)
                                     agaATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTtcatca                                                                                    <  1:383465/141‑1 (MQ=255)
                                       aaTCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAAc                                                                                  <  1:96994/141‑1 (MQ=255)
                                                  aTTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGaaa                                                                       >  1:353461/1‑141 (MQ=255)
                                                     aGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACAt                                                                    <  1:348855/141‑1 (MQ=255)
                                                        tAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATGGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGaa                                                                 >  1:255281/1‑141 (MQ=255)
                                                                   attatAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGAGTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCaaa                                                      >  1:208269/1‑141 (MQ=255)
                                                                          acttttAAATTTGAGGAGGAAGCGCTGTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTagag                                               <  1:173583/139‑1 (MQ=255)
                                                                               tAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTCTTAGATTCTGTAGGGCAAACACCTATGGTGTTACGTCATCAACTATTTCCGAAACATGAAGGGGTTGCAAAGTTAGAGTatat                                          >  1:388278/1‑141 (MQ=255)
                                                                                aaaTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATg                                         <  1:391959/141‑1 (MQ=255)
                                                                                aaaTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATg                                         >  1:377012/1‑141 (MQ=255)
                                                                                aaaTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCACCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATg                                         <  1:141228/141‑1 (MQ=255)
                                                                                  aTTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCGAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGaa                                       >  1:317659/1‑141 (MQ=255)
                                                                                               gcgcTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGgcagc                          >  1:158187/1‑141 (MQ=255)
                                                                                               gcgcTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACAACTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGgcagc                          >  1:442796/1‑141 (MQ=255)
                                                                                                  cTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATg                       <  1:306178/141‑1 (MQ=255)
                                                                                                       gattgaAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGa                  <  1:258163/141‑1 (MQ=255)
                                                                                                         ttgaAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAGCTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGGAAGATc                >  1:85129/1‑141 (MQ=255)
                                                                                                         ttgaAAAAAGTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATc                >  1:172943/1‑141 (MQ=255)
                                                                                                                  gTCAAGCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCa       <  1:483917/141‑1 (MQ=255)
                                                                                                                       gCATGTCACGATTCATTATTAGATTCTGTAGGGCAAACACCTATGGTTCAACTAAAATACCCATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCAAGTAc  <  1:492015/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                           
AATGCATAAAAAATATAATCCTAATTACTTGATAGTGAGAATCATTATCAATTAGGTAACACACAATATTATAGAATTTTAAATTTGAGGAGGAAGCGCTTTTGATTGAAAAAAGTCAAGCATGTCACGATTCATTGTTAGATTCTGTAGGGCAAACACCTATGGTTCAACTTCATCAACTATTTCCGAAACATGAAGTGTTTGCAAAGTTAGAGTATATGAATCCTGGAGGCAGCATGAAAGATCGACCTGCCAAGTAC  >  USA300TCH1516_ALE/134257‑134516

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