Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,241,539 T→C 57.6% intergenic (‑182/+284) murA2 ← / ← fba UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase 2/Fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,241,5390TC57.6% 1.0 / 35.2 31intergenic (‑182/+284)murA2/fbaUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase 2/Fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base C (9/9);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-02

TAAATTAATTTTACTTTATTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATATCACAATCAATTT  >  USA300TCH1516_ALE/2241411‑2241670
                                                                                                                                |                                                                                                                                   
taatttaattTTATTTTATTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTc                                                                                                                          <  1:103354/140‑1 (MQ=255)
      aattTTACTTTATTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATAATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGatta                                                                                                                    >  1:233532/1‑140 (MQ=255)
            aCTTTATTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAAc                                                                                                              >  1:66800/1‑140 (MQ=255)
            aCTTTATTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAAc                                                                                                              >  1:601399/1‑140 (MQ=255)
             cTTTTTTTGTAAAATAGTTTGTTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACt                                                                                                             <  1:171934/140‑1 (MQ=255)
                 aTTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAat                                                                                                         <  1:183672/140‑1 (MQ=255)
                     gTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTg                                                                                                     <  1:504570/140‑1 (MQ=255)
                        taaTAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCGTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACg                                                                                                  <  1:596112/139‑1 (MQ=255)
                            aGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTGCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACgttgt                                                                                              <  1:163534/140‑1 (MQ=255)
                                    ttAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTGTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAGCTAAATATTGACGTTGTTTTGACaa                                                                                      >  1:470000/1‑140 (MQ=255)
                                       aTCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTAACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAatta                                                                                   >  1:69080/1‑140 (MQ=255)
                                                   ttGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAAt                                                                       >  1:164805/1‑140 (MQ=255)
                                                   ttGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAAt                                                                       >  1:644440/1‑140 (MQ=255)
                                                      tCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTATTATATGCTGTTCTTCCCCAATTGACCACGGTCTAAATTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCt                                                                    <  1:407784/140‑1 (MQ=255)
                                                      tCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCt                                                                    <  1:85352/140‑1 (MQ=255)
                                                                   gattgattAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGt                                                       <  1:39302/140‑1 (MQ=255)
                                                                      tgattAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTTACGTaaa                                                    >  1:728925/1‑140 (MQ=255)
                                                                             taaCAAATCATTCATCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACcttttt                                             <  1:25852/140‑1 (MQ=255)
                                                                                    tCATTCACCGTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTGGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATaact                                      >  1:561507/1‑138 (MQ=255)
                                                                                     cATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAATTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAAtt                                     >  1:99515/1‑140 (MQ=255)
                                                                                       ttCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGCTTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGc                                   >  1:118303/1‑140 (MQ=255)
                                                                                            cTTTACTAATATATGCTGTTCTTCCCAACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTGAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCAtt                              <  1:636315/140‑1 (MQ=255)
                                                                                               tACTAATATATGCTGTTCTTCCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTc                           <  1:290191/140‑1 (MQ=255)
                                                                                                           cTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTACTTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATtata               >  1:289918/1‑140 (MQ=255)
                                                                                                           cTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATATTCCTTATAATtata               >  1:158402/1‑140 (MQ=255)
                                                                                                             gTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGTCTATTTTCCTTATAATTATATc             >  1:216524/1‑140 (MQ=255)
                                                                                                                cctcCCCACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATATcaca          <  1:449301/138‑1 (MQ=255)
                                                                                                                 ttcCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATATCAcaa         <  1:577866/140‑1 (MQ=255)
                                                                                                                   ccccACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCAAGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATATCAcaatc       >  1:399942/1‑140 (MQ=255)
                                                                                                                     ccACTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGCCGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATATCAcaatcaa     >  1:222571/1‑140 (MQ=255)
                                                                                                                        cTTGACCACGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATGTCACAATCAAttt  >  1:546432/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                   
TAAATTAATTTTACTTTATTTGTAAAATAGTTTATTTTAATCGTTAAAATCTTGTCTTTCTATTTTCGATTGATTAATAACAAATCATTCACCTTTACTAATATATGCTGTTCTTCCCCACTTGACCATGGTCTAAAGTCTGATTATTCAACTAAATATTGACGTTGTTTTGACAATTATTTAGTTTTAATTCTACGGTTTTCACGTAAACTTTTTTCATAAATTGCTCATTTTCCTTATAATTATATCACAATCAATTT  >  USA300TCH1516_ALE/2241411‑2241670

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