Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,042,574 T→C Y192C (TAC→TGC)  yqfA ← hemolysin‑III family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,5740TC97.6% 136.8 / ‑4.9 42Y192C (TAC→TGC) yqfAhemolysin‑III family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/13);  total (29/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGC  >  NC_000913/3042449‑3042710
                                                                                                                             |                                                                                                                                         
gACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTg                                                                                                                            <  2:199184/140‑1 (MQ=255)
       gcgcGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGt                                                                                                                     >  1:151615/1‑140 (MQ=255)
        cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTa                                                                                                                    <  1:529230/140‑1 (MQ=255)
            tCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGaaa                                                                                                                >  1:618309/1‑140 (MQ=255)
                 tCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCAc                                                                                                           >  2:705207/1‑140 (MQ=255)
                      ccTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGa                                                                                                      <  2:618309/140‑1 (MQ=255)
                       cTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAg                                                                                                     >  1:734397/1‑140 (MQ=255)
                       cTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAg                                                                                                     >  1:188527/1‑140 (MQ=255)
                         aCAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCg                                                                                                   >  1:555688/1‑140 (MQ=255)
                               gCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATaaa                                                                                             >  1:110557/1‑140 (MQ=255)
                               gCATACCCATTACGCCTGCCCAATATACAAAAGAACCCCAGAAAGTGGCAAAAACCACCGCCCAGCACGAAACCGTGGCAGAATGCAAGGTTGCATGGGATGGGGTTGCAAACGTAAAAAAACAACCCGAGCGAATAAAc                                                                                            >  2:327551/1‑140 (MQ=255)
                                cATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAAc                                                                                            >  1:350645/1‑140 (MQ=255)
                                  tACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCa                                                                                          >  1:593496/1‑140 (MQ=255)
                                   aCCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGa                                                                                                  >  1:250029/1‑131 (MQ=255)
                                   aCCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGa                                                                                                  <  2:250029/131‑1 (MQ=255)
                                   aCCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAc                                                                                         <  1:612979/140‑1 (MQ=255)
                                    cccATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcg                                                                                        <  2:614479/140‑1 (MQ=255)
                                     ccATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgc                                                                                       <  2:151615/140‑1 (MQ=255)
                                        ttACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcgccgc                                                                                    <  2:734397/140‑1 (MQ=255)
                                         tACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcg                                                                                        >  2:983212/1‑135 (MQ=255)
                                         tACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCAcg                                                                                        <  1:983212/135‑1 (MQ=255)
                                         tACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACgccgcc                                                                                   >  1:457320/1‑140 (MQ=255)
                                            cccTCCCAAATATAAAAATAGATCGCAAGAAGGTGGCCCCCCCTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCCCCCCGAGCGAATAACCCACGCCGCCTAc                                                                                <  2:188527/139‑1 (MQ=255)
                                             ccTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCCCCCCGAGCGAATAAACCACGCCGCCTAcc                                                                               <  1:382024/140‑1 (MQ=255)
                                                 cccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCa                                                                           >  1:88899/1‑140 (MQ=255)
                                                 cccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCa                                                                           >  2:670360/1‑140 (MQ=255)
                                                  ccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGACAg                                                                          >  1:945530/1‑140 (MQ=255)
                                                  ccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGACAg                                                                          <  2:945530/140‑1 (MQ=255)
                                                  caaaTATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAg                                                                          <  1:670360/138‑1 (MQ=255)
                                                    aaTATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTa                                                                        <  1:705207/140‑1 (MQ=255)
                                                                gATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGc                                                            >  2:235612/1‑140 (MQ=255)
                                                                 aTCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGcc                                                           >  2:966109/1‑140 (MQ=255)
                                                                  tCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGccc                                                          >  1:730161/1‑140 (MQ=255)
                                                                  tCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGccc                                                          >  1:857730/1‑140 (MQ=255)
                                                                   cGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCccg                                                         >  2:24690/1‑140 (MQ=255)
                                                                                       tACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCAtt                                     >  2:339980/1‑140 (MQ=255)
                                                                                         ccgccgAGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTc                                   >  1:819070/1‑140 (MQ=255)
                                                                                               aGCACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAt                             >  1:304768/1‑140 (MQ=255)
                                                                                                 cACGAAGCCGTGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAACCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTa                           >  2:187142/1‑140 (MQ=255)
                                                                                                           tGCCAGATGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGa                 >  2:814214/1‑140 (MQ=255)
                                                                                                                  tGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagccag          >  1:512012/1‑140 (MQ=255)
                                                                                                                  tGGCATGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagccag          >  1:337299/1‑140 (MQ=255)
                                                                                                                     gcgggTTGGAAGGAATTGCTTTTGCAGAGGAAAAAAACCCCCCCGAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACagccagcc        <  2:88899/137‑1 (MQ=255)
                                                                                                                      aTGGTTGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCa      >  1:572273/1‑140 (MQ=255)
                                                                                                                          ttGCATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGc  >  2:727903/1‑140 (MQ=255)
                                                                                                                             |                                                                                                                                         
GACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGC  >  NC_000913/3042449‑3042710

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