Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,440,085 G→T intergenic (‑45/+126) cspA_2 ← / ← USA300TCH1516_01338 Cold shock protein CspA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,440,0850GT93.2% 134.8 / 0.9 45intergenic (‑45/+126)cspA_2/USA300TCH1516_01338Cold shock protein CspA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (16/25);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATCAA  >  USA300TCH1516_ALE/1439948‑1440218
                                                                                                                                         |                                                                                                                                     
tGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCat                                                                                                                                    <  1:181947/141‑1 (MQ=255)
   aaaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATAtt                                                                                                                                 <  1:898332/141‑1 (MQ=255)
     aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTAc                                                                                                                               <  1:384887/141‑1 (MQ=255)
     aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAAGCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTAc                                                                                                                               <  2:320876/141‑1 (MQ=255)
        gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaa                                                                                                                            <  2:14605/141‑1 (MQ=255)
           cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaaaaa                                                                                                                         >  1:583595/1‑141 (MQ=255)
           cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAAt                                                                                                                        >  2:174909/1‑141 (MQ=255)
               tACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACa                                                                                                                     <  1:312422/141‑1 (MQ=255)
                       tttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAAtt                                                                                                             >  1:151368/1‑141 (MQ=255)
                           tcCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCAc                                                                                                         >  1:334070/1‑141 (MQ=255)
                                      cGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCt                                                                                              <  1:470146/141‑1 (MQ=255)
                                        aTAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt                                                                                            <  2:223015/141‑1 (MQ=255)
                                             gCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAat                                                                                       >  2:252805/1‑141 (MQ=255)
                                               cGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat                                                                                     <  1:116220/141‑1 (MQ=255)
                                               cGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat                                                                                     <  1:348974/141‑1 (MQ=255)
                                                  aTCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGa                                                                                  >  1:506632/1‑141 (MQ=255)
                                                  aTCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGa                                                                                  >  2:184943/1‑141 (MQ=255)
                                                       ttttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaataat                                                                             <  2:961/141‑1 (MQ=255)
                                                              cGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGc                                                                      <  1:665673/141‑1 (MQ=255)
                                                                    aCCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACt                                                                >  1:867687/1‑141 (MQ=255)
                                                                         ttAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtat                                                           <  1:550830/141‑1 (MQ=255)
                                                                         ttAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtat                                                           <  1:855356/141‑1 (MQ=255)
                                                                            aCTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtattat                                                        >  1:899852/1‑141 (MQ=255)
                                                                               gTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGt                                                     >  1:735107/1‑141 (MQ=255)
                                                                                tACCTTATTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCATTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTg                                                    <  1:497681/141‑1 (MQ=255)
                                                                                      gTTTCATAATCTGAACCCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGtta                                              >  2:417734/1‑141 (MQ=255)
                                                                                          cATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTAc                                          <  2:133153/141‑1 (MQ=255)
                                                                                            tAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCa                                        >  2:277501/1‑141 (MQ=255)
                                                                                                cTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATAGTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAAtt                                    <  2:167496/141‑1 (MQ=255)
                                                                                                  gAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGc                                  <  2:103478/141‑1 (MQ=255)
                                                                                                  gAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGc                                  <  2:192068/141‑1 (MQ=255)
                                                                                                  gAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGc                                  <  2:196343/141‑1 (MQ=255)
                                                                                                   aaaCCTCCAAGACTAAAATTCATTCAATATGCTTATTCCCAGATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGca                                 >  2:414135/1‑141 (MQ=255)
                                                                                                       cTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatac                             <  2:38707/141‑1 (MQ=255)
                                                                                                         ccAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatacaa                           >  1:281179/1‑141 (MQ=255)
                                                                                                         ccAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATACTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatacaa                           >  2:149277/1‑141 (MQ=255)
                                                                                                          cAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAGTAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatacaat                          >  1:61578/1‑141 (MQ=255)
                                                                                                                 aaaaTTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATaaaa                   >  1:146590/1‑141 (MQ=255)
                                                                                                                  aaaTTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATaaaaa                  <  1:718425/141‑1 (MQ=255)
                                                                                                                    attcattcAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAat                <  1:771543/141‑1 (MQ=255)
                                                                                                                       cattcaATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATAtt             >  2:55048/1‑141 (MQ=255)
                                                                                                                             atatGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTa       >  1:255931/1‑141 (MQ=255)
                                                                                                                              tatGCTTCTTATAATATTACACAACATACATGCACCTTACCCTCACAAATAATTTGCAATATGGCATACACATCAGACTAACCTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTaa      <  1:492331/141‑1 (MQ=255)
                                                                                                                               atGCTTATTCGCATATGACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAAt     <  1:32912/141‑1 (MQ=255)
                                                                                                                                  cTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATCaa  <  1:170807/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATCAA  >  USA300TCH1516_ALE/1439948‑1440218

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