Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,393,232 A→G 19.7% S89P (TCA→CCA)  moeA ← Molybdopterin molybdenumtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,393,2320AG19.7% 53.2 / 6.7 33S89P (TCA→CCA) moeAMolybdopterin molybdenumtransferase
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (3/3);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

GGTTTACGAATTGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTACAAT  >  USA300TCH1516_ALE/2393096‑2393367
                                                                                                                                        |                                                                                                                                       
ggTTTACGAATTGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCt                                                                                                                                     <  1:23997/141‑1 (MQ=255)
           tGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTCCACCAATATg                                                                                                                          >  1:190881/1‑141 (MQ=255)
                             ttCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGGACCTGCACCCATACGATCAATCACTTTAAACTc                                                                                                        >  2:154794/1‑141 (MQ=255)
                             ttCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTc                                                                                                        >  2:158085/1‑141 (MQ=255)
                             ttCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCACCAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTc                                                                                                        >  2:216739/1‑141 (MQ=255)
                                       tGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACGAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTc                                                                                              <  1:87081/141‑1 (MQ=255)
                                           cAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTTGTGATCCCCAACTAATTTATCAGAAACTGAACCCCCCCCACTATGATAAATCACTTTAAACTAAAAGCGATTCtcaa                                                                                          >  2:96252/1‑138 (MQ=255)
                                                     aacaGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGcccc                                                                                <  2:238005/141‑1 (MQ=255)
                                                            tCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACAAATTAATCCTGAACTGAACAAGCAACAATATGATCAATCACTTTAAACTCAATGCGATACTAACAACTTGCCCCTTGTaaa                                                                         >  2:141017/1‑141 (MQ=255)
                                                              tGCGCCATTCGGTATTTGTGCTACAATGGAAAAACGCACCGCTTCGTGATCCACAACTAATTTATCGGAAACTGAACCTGCACCAATATTATCAATCACTTTAATCTCAATGCGATTCTGACCACTTGCCCCTTGTGAATc                                                                       <  1:56465/141‑1 (MQ=255)
                                                                     ttAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCACCTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTACACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTg                                                                >  2:85696/1‑141 (MQ=255)
                                                                            tttGTGCTCTAGTCATAATACGCACCGCTTCGTGATCCCCTACTAACTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTg                                                         <  2:213916/141‑1 (MQ=255)
                                                                            tttGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTg                                                         <  2:45353/141‑1 (MQ=255)
                                                                              tgtgCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCa                                                       >  1:268770/1‑141 (MQ=255)
                                                                              tgtgCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACACCACCCCCAATTCGATTCCGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCa                                                       >  2:267867/1‑141 (MQ=255)
                                                                               gtgCTCCAGTCATAAGACGCACCGCTTCGCGATCCCCAACTAATTTATCAGAAACTGTACCTGCACCAATATGATCAATCACTTTAAACTCCATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCaa                                                      <  2:260785/141‑1 (MQ=255)
                                                                                           tgataCGGACAGCTTCGTGATTACCAACTAGGTTATCAGAAACTGCACTTGCACCAATATGATCAATCACTTGAAGCTCAGTGCGATTGTGGCGACTTGGCCCTTGTGTTTCAACACTGGGTATTGCAGAGCCataatatg                                          <  2:255494/139‑3 (MQ=255)
                                                                                            aataCGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAgg                                         >  1:209004/1‑141 (MQ=255)
                                                                                                    ccGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTAt                                 >  1:221697/1‑141 (MQ=255)
                                                                                                     tgCTTCTTGATACTCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATc                                <  2:135470/140‑1 (MQ=255)
                                                                                                        ttCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCaaa                             >  1:63765/1‑141 (MQ=255)
                                                                                                        ttCGTGATCCCCAACTAATTTATCAGAAACTGAACCCTCACCAATATGATCAATCACTTTAAACTCAATTCTTTTCTGACCACTTTCCNCTTGTGAATCAACACTCCGAATTGCAAAACCACCATAAGGCCCTTTCTCaaa                             >  2:75429/1‑141 (MQ=255)
                                                                                                               tCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGt                      <  1:269099/141‑1 (MQ=255)
                                                                                                                tcccAACTAATTTATCAGAAACTTAACCTACATCGATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACGTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGta                     <  2:203328/140‑1 (MQ=255)
                                                                                                                   cAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTatat                  <  1:158085/141‑1 (MQ=255)
                                                                                                                   cAACTAATTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTatat                  <  1:216739/141‑1 (MQ=255)
                                                                                                                   cAACTAAGTTATCAGAAACTGGACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTatat                  <  1:154794/141‑1 (MQ=255)
                                                                                                                     aCTAATTCATCAGTCACTGAACCTGCACCAGTATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCa                <  1:131470/141‑1 (MQ=255)
                                                                                                                          tttATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGt           >  2:198593/1‑141 (MQ=255)
                                                                                                                             aTCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTACCAAACCTGGGTGTATCATAAGTAGc        >  1:107652/1‑141 (MQ=255)
                                                                                                                                aGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTAc     <  1:157659/141‑1 (MQ=255)
                                                                                                                                aGAAACTGAACCGGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTAc     <  2:47865/141‑1 (MQ=255)
                                                                                                                                   aaCTGAACCTGCACCAATATGATCAATCCCTTTAAACTCAATGCGATTCTGACCACTTGCCCCCTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTACAat  <  2:209822/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GGTTTACGAATTGTAAATGTATCTTCTAGTTCAATCGTTTGTTCAAACATAACAACAGCATCTGCGCCATTAGGTATTTGTGCTCCAGTCATAATACGCACCGCTTCGTGATCCCCAACTAATTTATCAGAAACTGAACCTGCACCAATATGATCAATCACTTTAAACTCAATGCGATTCTGACCACTTGCCCCTTGTGAATCAACACTGCGAATTGCAAAACCATCATAAGGTGATTTATCAAACCTTGGTATATCATAAGTAGCTACAAT  >  USA300TCH1516_ALE/2393096‑2393367

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