Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,519,951 T→G 11.3% P227P (CCA→CCC narG ← Respiratory nitrate reductase 1 alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,519,9510TG11.3% 76.0 / 2.5 38P227P (CCA→CCCnarGRespiratory nitrate reductase 1 alpha chain
Reads supporting (aligned to +/- strand):  ref base T (21/11);  new base G (3/2);  total (24/14)
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

CTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATC  >  USA300TCH1516_ALE/2519815‑2520088
                                                                                                                                        |                                                                                                                                         
cTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGACACATCAGTTTGCTCTCCCCAAATTTGTggca                                                                                                                                       >  2:143229/1‑139 (MQ=255)
      aGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCCTTATTCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAg                                                                                                                                 <  1:307350/141‑1 (MQ=255)
         ggACTTCAGTCATAAAAGGTCTATCCGGAGTACTTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCATGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggtg                                                                                                                              <  2:62574/141‑1 (MQ=255)
                aGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAACTTTGTGGAGAGGCAGCTGGTAAATc                                                                                                                       >  2:3785/1‑141 (MQ=255)
                        aaTGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATAcc                                                                                                               <  1:7338/141‑1 (MQ=255)
                              cATCCGGAGTACGTGTTGAAGGTACATTATAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                                                                         <  2:287614/141‑1 (MQ=255)
                              aatacGGAGTACTTTGTAAAGGTACAGTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                                                                         <  1:243271/137‑1 (MQ=255)
                                      gTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGGCATAAAAACTa                                                                                                 >  1:195518/1‑141 (MQ=255)
                                         cGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTAAGGCACATCTGTTTGCTCTCCCCGAATTTCTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAACAACTAAGc                                                                                              >  2:93236/1‑141 (MQ=255)
                                                aaGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTcac                                                                                       >  1:220965/1‑141 (MQ=255)
                                                       aTTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAg                                                                                <  2:63703/141‑1 (MQ=255)
                                                         tAGAGCTCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTCGCTCTCCCCAAATGGGTGGAGAGGCAGGTAGTAGATCTGCATACCAGTCATAAAAACTAAGGATTTCAGCACCAAGCa                                                                              >  2:265229/1‑141 (MQ=255)
                                                             gCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGGGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCGGCAGCAAGCAAAtt                                                                          >  1:116473/1‑141 (MQ=255)
                                                                ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTTTCCCCAAATTTGTGGACAGGCAGGTGGTAAATCTGCATACCAGTCACAAAAACTTACAATTTCACAACCAAGCCAATtgat                                                                       >  2:147176/1‑141 (MQ=255)
                                                                    ataataATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAAt                                                                   >  1:25306/1‑141 (MQ=255)
                                                                       ataatGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGa                                                                >  1:181636/1‑141 (MQ=255)
                                                                        taatGTATGATGCGTTATAACAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGCATCGaa                                                               >  2:337141/1‑140 (MQ=255)
                                                                              atggtgCGTTATACCAGTCAGTTGATTCGGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctg                                                         <  2:55655/141‑1 (MQ=255)
                                                                                gatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgct                                                       >  1:308220/1‑141 (MQ=255)
                                                                                 atgCGTTATACCAGTCACTTGATTCAGGCACATCTTTTTTTTTTTTTTTTTTTTATGCAGAGGCAGGTGCTAAATCTTCATACCAGTCATAAAAACTAAGCATTTCTCCACCACGCATCTTGATGAATCGAGCActtcctc                                                      >  2:195092/1‑135 (MQ=255)
                                                                                   gCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCa                                                    >  2:30988/1‑141 (MQ=255)
                                                                                                  cTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGAc                                     >  1:74393/1‑141 (MQ=255)
                                                                                                    tGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACAt                                   >  1:290826/1‑141 (MQ=255)
                                                                                                    tGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACAt                                   >  1:289131/1‑141 (MQ=255)
                                                                                                      aTTCAGGCACAGCGGGTGGCTCTCAGCATTTTTGGGGAGAGGCAGGTGGTAAATCTGCAGACCAGTCATAAAAACTAAGCATTTCACCACCAAGGAAATTGATGAACCGAGCAACTGCTGCATAACTAATAATTGACATCg                                 <  2:167480/141‑1 (MQ=255)
                                                                                                        tCAGGCACATCTGTTTGCTCTCCCCCAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCt                               <  1:105889/141‑1 (MQ=255)
                                                                                                        acaagCACATCTGTTTGATCTCCCCAAATTTGTGAAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCt                               <  2:250009/137‑1 (MQ=255)
                                                                                                             cacaTCTGTTTGTTCTCACCACATCTACGCAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAAt                          <  1:264362/141‑1 (MQ=255)
                                                                                                                  cTGTTTTCTCTCCCCAAAGTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGgtg                     <  1:291761/141‑1 (MQ=255)
                                                                                                                       tGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCCTTCACAGCGCCGTAACAGgtgtccat                >  2:274361/1‑137 (MQ=255)
                                                                                                                        gCTCTCCCCAAAGTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGCATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATc               <  1:95784/141‑1 (MQ=255)
                                                                                                                              cccAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGa         >  1:171296/1‑141 (MQ=255)
                                                                                                                              cacAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAAATAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCCTTACCATCGCTGGAATAGGTGTAAATCCTGCGa         >  1:322/3‑141 (MQ=255)
                                                                                                                                 aaaTTTGGCGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTACGCATTTCACCAACAAGGAAATTGATGAATCGAGCACCTGCTGCATATCTAATCATTGACATCGCTCGAATAGGTGTAAATCCTGCGATTc      <  2:220544/141‑1 (MQ=255)
                                                                                                                                   aTTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGa    >  2:97201/1‑141 (MQ=255)
                                                                                                                                   aTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGa    >  1:179472/1‑141 (MQ=255)
                                                                                                                                   aTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAACATTTCTCCACAAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGa    >  1:178997/1‑141 (MQ=255)
                                                                                                                                     ttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTACTCATTGCCATCGCTGGAATCGGTTTTAATACTGCGATTCGATc  >  2:294882/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                         
CTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATC  >  USA300TCH1516_ALE/2519815‑2520088

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: