Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 915,857 C→T 9.0% A5V (GCA→GTA)  sufB_2 → FeS cluster assembly protein SufB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE915,8570CT 9.0% 135.2 / 3.2 44A5V (GCA→GTA) sufB_2FeS cluster assembly protein SufB
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base T (2/2);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GATATTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGATTTCCGT  >  USA300TCH1516_ALE/915719‑915996
                                                                                                                                          |                                                                                                                                           
gATATTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCAc                                                                                                                                           <  2:258849/141‑1 (MQ=255)
  tatTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCt                                                                                                                                         <  2:64515/141‑1 (MQ=255)
            ataAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTgggg                                                                                                                               <  1:203833/141‑1 (MQ=255)
                     taAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAat                                                                                                                      <  1:267580/141‑1 (MQ=255)
                          aTTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGa                                                                                                                 <  1:238822/141‑1 (MQ=255)
                               tataAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCa                                                                                                            >  2:398100/1‑141 (MQ=255)
                               tataAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCa                                                                                                            >  2:180584/1‑141 (MQ=255)
                                 taAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCacg                                                                                                          >  2:446917/1‑141 (MQ=255)
                                  aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         >  1:201020/1‑141 (MQ=255)
                                  aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         >  1:293493/1‑141 (MQ=255)
                                  aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         <  2:315014/141‑1 (MQ=255)
                                  aaGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCacga                                                                                                         <  2:408871/141‑1 (MQ=255)
                                            gtcatTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGAtgtg                                                                                               <  1:7580/141‑1 (MQ=255)
                                                 tGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCAt                                                                                          <  1:185976/141‑1 (MQ=255)
                                                 tGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCAt                                                                                          <  1:107454/141‑1 (MQ=255)
                                                  gTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCAtt                                                                                         <  1:6942/141‑1 (MQ=255)
                                                     ataataTGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtc                                                                                      >  1:286996/1‑141 (MQ=255)
                                                        atatGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcaga                                                                                   >  1:219818/1‑141 (MQ=255)
                                                          atGGTTTACATCATGAATTAAAAACTTACGCACGTCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatc                                                                                 <  2:136581/141‑1 (MQ=255)
                                                          atGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatc                                                                                 <  2:6513/141‑1 (MQ=255)
                                                          atGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatc                                                                                 <  1:309131/141‑1 (MQ=255)
                                                            ggTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTtcagatcag                                                                               >  1:301083/1‑141 (MQ=255)
                                                                tACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACg                                                                           >  2:266811/1‑141 (MQ=255)
                                                                tACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACg                                                                           >  2:425658/1‑141 (MQ=255)
                                                                 aCATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGt                                                                          >  1:32604/1‑141 (MQ=255)
                                                                    tcatGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGt                                                                       >  1:141458/1‑141 (MQ=255)
                                                                          aTTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACt                                                                 <  2:227123/141‑1 (MQ=255)
                                                                              aaaaCTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTgaga                                                             >  2:279169/1‑141 (MQ=255)
                                                                              aaaaCTTACGCACGCCGTGGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTgaga                                                             <  2:113991/141‑1 (MQ=255)
                                                                                     acgcacgcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGt                                                      <  1:446917/141‑1 (MQ=255)
                                                                                        cacgcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTag                                                   <  2:6188/141‑1 (MQ=255)
                                                                                        cacgcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTag                                                   <  2:398844/141‑1 (MQ=255)
                                                                                           gcCGTTGTAAATATATTTTTAAGGGGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAgaga                                                <  2:201020/141‑1 (MQ=255)
                                                                                           gcCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGTACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAgaga                                                <  2:293493/141‑1 (MQ=255)
                                                                                                tGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAAttt                                           >  2:248148/1‑141 (MQ=255)
                                                                                                          tttttAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTCGGGGACGACAAAGCTGGATTCCACGACGATGATGTGTCCATTTTCAGACCCGAACGGGGGTCAACTGAGAATATCGTTAGAGAAACTTCTAAcgcgaa                                 >  2:69249/1‑136 (MQ=255)
                                                                                                               aaGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATg                            <  1:40312/141‑1 (MQ=255)
                                                                                                                       aTTGAAATGACTAAAAAAGCAACTGATGTTGGGGATTATAAATATGGATTCCACGACAATGATGTGTCCATTTTCAGATCAGAACGTGGTTTATCTGTGAAGATCGTTAGAGAAATAAAAAAAAAAAAAAATGAGCCGGaa                    <  2:376426/141‑1 (MQ=255)
                                                                                                                           aaaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGa                >  2:182734/1‑141 (MQ=255)
                                                                                                                            aaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGAt               >  2:25825/1‑141 (MQ=255)
                                                                                                                            aaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGAt               >  1:323173/1‑141 (MQ=255)
                                                                                                                            aaTGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGAt               >  1:237713/1‑141 (MQ=255)
                                                                                                                                     aaaaGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGAttt      <  2:296909/141‑1 (MQ=255)
                                                                                                                                         gCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGATTTCCGt  <  2:219818/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GATATTAGACATATAAAAGTATAAAAATTTTTATAAGATGTCATGTCATTGTTATAATATGGTTTACATCATGAATTAAAAACTTACGCACGCCGTTGTAAATATATTTTTAAGGAGTGATTGAAATGGCTAAAAAAGCACCTGATGTTGGGGATTATAAATATGGATTCCACGACGATGATGTGTCCATTTTCAGATCAGAACGTGGTTTAACTGAGAATATCGTTAGAGAAATTTCTAACATGAAAAATGAGCCGGAATGGATGTTAGATTTCCGT  >  USA300TCH1516_ALE/915719‑915996

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