Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,036,314 A→G F65S (TTC→TCC)  yqfA ← hemolysin III family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,036,3140AG96.2% 82.7 / ‑6.1 26F65S (TTC→TCC) yqfAhemolysin III family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (12/13);  total (12/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

AATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCC  >  NZ_CP009273/3036183‑3036438
                                                                                                                                   |                                                                                                                            
aaTCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGGAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATAGCGTGATAGAGCGTCGAAGCGAGGGAcagcag                                                                                                                       >  1:497762/1‑139 (MQ=255)
                  aaaaCGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATa                                                                                                     <  2:218997/139‑1 (MQ=255)
                         cGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTg                                                                                              >  1:50142/1‑139 (MQ=255)
                                  tCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAg                                                                                           >  2:579725/1‑133 (MQ=255)
                                  tCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAg                                                                                           <  1:579725/133‑1 (MQ=255)
                                              cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                         >  2:475169/1‑139 (MQ=255)
                                                           cAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCAtt                                                                 >  1:101524/1‑134 (MQ=255)
                                                                          ttcAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGAt                                                                                                                    <  1:292797/68‑1 (MQ=255)
                                                                          ttcAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGAt                                                                                                                    >  2:292797/1‑68 (MQ=255)
                                                                            cAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCa                                             <  1:30000/137‑1 (MQ=255)
                                                                                    tttttGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAacac                                   <  2:50142/139‑1 (MQ=255)
                                                                                            cGTTGATGAGGAATGCGGTGATAGAGCGTGAAAGGGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACCCCAGCCCaa                           <  2:512231/139‑1 (MQ=255)
                                                                                             gTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTaa                                                               >  2:97504/1‑102 (MQ=255)
                                                                                             gTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTaa                                                               <  1:97504/102‑1 (MQ=255)
                                                                                             gTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAAc                          >  2:551254/1‑139 (MQ=255)
                                                                                              ttGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAAc                          >  2:497020/1‑138 (MQ=255)
                                                                                              ttGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTACCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACg                         <  2:101524/139‑1 (MQ=255)
                                                                                                gatgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGa                        <  1:184913/138‑1 (MQ=255)
                                                                                                gatgaGGAATGGCGTGATAGAGCAACGAAGCGAGGGAGAGCGGGATAATACTGCCGCCATAGAGGCTGTAGCTGGGTATCGCTGGGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGa                        <  2:400455/138‑1 (MQ=255)
                                                                                                    aGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                   <  1:497020/139‑1 (MQ=255)
                                                                                                         tGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAacac              >  1:210425/1‑139 (MQ=255)
                                                                                                                aTAGAGCGTCTAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAAccc        >  2:186143/1‑138 (MQ=255)
                                                                                                                aTAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCa       <  2:179167/139‑1 (MQ=255)
                                                                                                                     gCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGcc  <  1:73226/139‑1 (MQ=255)
                                                                                                                          gAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                   <  2:154179/117‑1 (MQ=255)
                                                                                                                          gAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                   >  1:154179/1‑117 (MQ=255)
                                                                                                                                   |                                                                                                                            
AATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCC  >  NZ_CP009273/3036183‑3036438

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