Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,359,587 G→T R215S (CGT→AGT)  yjeH ← L‑methionine/branched‑chain amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,359,5870GT94.3% 123.3 / 0.2 35R215S (CGT→AGT) yjeHL‑methionine/branched‑chain amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (18/15);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

CAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGC  >  NZ_CP009273/4359455‑4359719
                                                                                                                                    |                                                                                                                                    
cAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTg                                                                                                                                 >  1:329600/1‑138 (MQ=255)
                  ctgcCGCCGCCATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAgg                                                                                                               <  2:282023/138‑1 (MQ=255)
                          gccATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATggg                                                                                                      >  2:203221/1‑139 (MQ=255)
                             aTTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCaa                                                                                                   <  2:250097/139‑1 (MQ=255)
                                         tAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGAc                                                                                       <  2:118982/139‑1 (MQ=255)
                                                cGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACaa                                                                                 >  1:140521/1‑138 (MQ=255)
                                                 gAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACaaaa                                                                               >  1:49894/1‑139 (MQ=255)
                                                 gAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACaaaa                                                                               >  1:36484/1‑139 (MQ=255)
                                                    gtgtAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAacac                                                                            <  1:211528/139‑1 (MQ=255)
                                                    gtgtAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAacac                                                                            <  2:208828/139‑1 (MQ=255)
                                                    gtgtAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAacac                                                                            <  2:149193/139‑1 (MQ=255)
                                                       tAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAacac                                                                            >  2:336803/1‑136 (MQ=255)
                                                            cgacTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACAt                                                                    <  1:356702/139‑1 (MQ=255)
                                                               cTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAacac                                                                            >  2:46884/1‑128 (MQ=255)
                                                               cTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAacac                                                                            <  1:46884/128‑1 (MQ=255)
                                                               cTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCAc                                                                 >  1:156926/1‑139 (MQ=255)
                                                                     tACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATaa                                                           >  2:411513/1‑139 (MQ=255)
                                                                        aGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACgc                                                        <  1:148085/139‑1 (MQ=255)
                                                                          ccccAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGAt                                                                                                         <  2:331792/88‑1 (MQ=255)
                                                                          ccccAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGAt                                                                                                         >  1:331792/1‑88 (MQ=255)
                                                                            ccAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACgcagc                                                     >  1:366618/1‑138 (MQ=255)
                                                                               gTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCaaa                                                  >  1:136401/1‑138 (MQ=255)
                                                                                       aTCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCgg                                          >  1:228246/1‑138 (MQ=255)
                                                                                       aTCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGt                                         <  1:203221/139‑1 (MQ=255)
                                                                                        tCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGt                                         >  2:204665/1‑138 (MQ=255)
                                                                                          cTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAg                                      >  2:396300/1‑139 (MQ=255)
                                                                                             ccAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGt                                     <  2:215078/137‑1 (MQ=255)
                                                                                             ccAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTc                                   >  2:31758/1‑139 (MQ=255)
                                                                                                       aCCAATCATCAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAc                           >  1:14170/1‑137 (MQ=255)
                                                                                                           atcatcAAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTg                     >  1:370676/1‑139 (MQ=255)
                                                                                                                   aGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACaa                                                                                 <  2:194361/71‑1 (MQ=255)
                                                                                                                   aGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACaa                                                                                 >  1:194361/1‑71 (MQ=255)
                                                                                                                   aGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTGCGATATTACCAGGTGCCGGAAAg             <  1:67588/139‑1 (MQ=255)
                                                                                                                             aaaTCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCg   <  1:355346/139‑1 (MQ=255)
                                                                                                                              aaTCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGc  <  1:84846/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGC  >  NZ_CP009273/4359455‑4359719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑