Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,615,188 A→C S69R (AGT→AGG eamA ← O‑acetylserine/cysteine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,615,1880AC92.3% 82.4 / ‑0.2 26S69R (AGT→AGGeamAO‑acetylserine/cysteine exporter
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (8/16);  total (8/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 = 3.97e-01

AAACGCCTGTGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGAC  >  NZ_CP009273/1615101‑1615268
                                                                                       |                                                                                
aaaCGCCTGTGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCt                                                                                >  2:909509/1‑90 (MQ=255)
  aCGCCTGTGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGa                                                                              <  2:339882/90‑1 (MQ=255)
  aCGCCTGTGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGa                                                                              <  1:95598/90‑1 (MQ=255)
  aCGCCTGTGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGa                                                                              <  1:906206/90‑1 (MQ=255)
  aCGCCTGAGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGa                                                                              <  1:67471/90‑1 (MQ=255)
       tgtgcctgtAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGtt                                                                         <  2:365748/90‑1 (MQ=255)
       tgtgcctgtAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGtt                                                                         <  1:109644/90‑1 (MQ=255)
       tgtgcctgtAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGtt                                                                         <  1:475989/90‑1 (MQ=255)
        gtgcctgtAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTa                                                                        <  2:771868/90‑1 (MQ=255)
         tgcctgtAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTaa                                                                       <  1:909509/90‑1 (MQ=255)
         tgcctgtAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTaa                                                                       <  2:645692/90‑1 (MQ=255)
                          aaGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAgcagca                                                      <  2:303535/90‑1 (MQ=255)
                           aGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCaa                                                     <  1:337373/90‑1 (MQ=255)
                           aGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCaa                                                     >  1:195882/1‑90 (MQ=255)
                            gCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGa                                                                              >  1:70346/1‑64 (MQ=255)
                            gCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGa                                                                              <  2:70346/64‑1 (MQ=255)
                             ccAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAAt                                                   <  1:526444/90‑1 (MQ=255)
                                gTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCa                                                >  2:94712/1‑90 (MQ=255)
                                    agcagGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTgg                                            <  2:724465/90‑1 (MQ=255)
                                                  aGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTg                              >  2:700759/1‑90 (MQ=255)
                                                           cacaAAAAAGAAAAGCAAACTGCGCAAACGTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACaaaaa                     >  1:732855/1‑90 (MQ=255)
                                                           cacaAAAAAGAAAAGCAAACTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACaaaaa                     >  2:183986/1‑90 (MQ=255)
                                                                             aCTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGa   >  1:838462/1‑90 (MQ=255)
                                                                              cTGCGCAAGCCGGAAGGTTAATCCATCCCCCAGCAGCAAATTCAGTCGTCCTTTCGGTCGTGCGACAAAAAAGGTAGCCGGAAAAGCGAc  <  2:13266/90‑1 (MQ=255)
                                                                              cTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGAc  <  2:195882/90‑1 (MQ=255)
                                                                              cTGCGCAAACCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGAc  <  2:838462/90‑1 (MQ=255)
                                                                                       |                                                                                
AAACGCCTGTGCCTGTAACACCAGCGAAGCCAGTCCAGCAGGCATACCGAAGTTAATGGCACAAAAAAGAAAAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGAC  >  NZ_CP009273/1615101‑1615268

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