Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 996,768 G→A S59S (TCC→TCT)  ssuA ← aliphatic sulfonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913996,7680GA97.1% 125.4 / NA 35S59S (TCC→TCT) ssuAaliphatic sulfonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/16);  minor base C (1/0);  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 = 1.00e+00

GGGCTTCGGTGGCTCGACGCCCACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAGGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGCTAAGCAATCCC  >  NC_000913/996621‑996915
                                                                                                                                                   |                                                                                                                                                   
gggCTTCGGTGGCTCGACGCCCACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAagag                                                                                                                                                  <  2:569285/151‑1 (MQ=255)
     tCGGTGGCTCGACGCCCACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGAtt                                                                                                                                               <  2:649100/149‑1 (MQ=255)
                aCGCCCACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGAt                                                                                                                                  <  2:357692/151‑1 (MQ=255)
                 cGCCCACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATa                                                                                                                                 <  2:349139/151‑1 (MQ=255)
                    ccACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCg                                                                                                                              >  1:662411/1‑151 (MQ=255)
                        gTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGtttt                                                                                                                          >  2:105721/1‑151 (MQ=255)
                            acCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTcc                                                                                                                       >  2:387614/1‑150 (MQ=255)
                                          ggCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCtt                                                                                                        >  2:668746/1‑151 (MQ=255)
                                          ggCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCtt                                                                                                        >  2:315608/1‑151 (MQ=255)
                                                             ggCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCaa                                                                                     <  2:663755/151‑1 (MQ=255)
                                                             ggCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCACGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCa                                                                                      >  2:240737/1‑150 (MQ=255)
                                                                ggAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATAc                                                                                  <  2:350214/151‑1 (MQ=255)
                                                                        cccGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTc                                                                          >  1:19151/1‑151 (MQ=255)
                                                                            gTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGAt                                                                      >  1:236997/1‑151 (MQ=255)
                                                                              ccTGACGAGATCATTACTCCCTACGTTTAACGCTGCGAACATGAACGGACCCGCGGGGGGCTCCACCCACGAGCTTTGTGATTCCGGATCACGTTGGTCAAGCAACTTTCGGCATTTTGCCAGCACCATACCAATACGCCCGTTATGATa                                                                     >  1:154838/2‑150 (MQ=255)
                                                                                          aaTACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAAc                                                        >  2:617369/1‑151 (MQ=255)
                                                                                              cTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCtt                                                    >  2:112926/1‑151 (MQ=255)
                                                                                                 ccAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAg                                                 <  2:472916/151‑1 (MQ=255)
                                                                                                   aaCGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGc                                               >  1:8468/1‑151 (MQ=255)
                                                                                                      gTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAg                                            <  1:337238/151‑1 (MQ=255)
                                                                                                         tAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGAt                                         >  1:334177/1‑151 (MQ=255)
                                                                                                          aaCGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGAt                                         <  2:178121/150‑1 (MQ=255)
                                                                                                            cGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCt                                      <  1:5596/151‑1 (MQ=255)
                                                                                                            cGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCt                                      >  2:344944/1‑151 (MQ=255)
                                                                                                            cGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCt                                      >  1:673786/1‑151 (MQ=255)
                                                                                                              cTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTgc                                    >  1:394535/1‑151 (MQ=255)
                                                                                                                tCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTgca                                   <  2:489099/150‑1 (MQ=255)
                                                                                                                tCCAACATTTGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTgca                                   <  2:118014/150‑1 (MQ=255)
                                                                                                                        ttGCGGACCCGCGGGGAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGc                           <  2:546706/150‑1 (MQ=255)
                                                                                                                          gCGGACCCGCGGGGAACTCCACCCAAGAGATTTGTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGTCTTTTTGCCAGCACCATACCAATACTGCCTGTCTGATAGCCTATGCGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCaa                         <  1:204424/150‑1 (MQ=255)
                                                                                                                                    ggggAACTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACg              <  2:141181/151‑1 (MQ=255)
                                                                                                                                        aaCTCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGCTaa          <  2:41969/151‑1 (MQ=255)
                                                                                                                                           tCCACCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGCTAAGc        >  2:625384/1‑150 (MQ=255)
                                                                                                                                              aCCCAAGAGATTTTTGATTCCGGATAGCGTTTTTCAAGTAACTGGTGGCTTTTTGCCAGTACCATACCAATAGTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGATAAGCAAt     >  1:254599/1‑150 (MQ=255)
                                                                                                                                               cccAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGCTAAGCAATc    >  1:8362/1‑150 (MQ=255)
                                                                                                                                                ccAAGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGCTAAGCAATccc  <  1:543792/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGCTTCGGTGGCTCGACGCCCACGTACACCAAATCAGCCCCGGCAGCCTGGGCAAAGATTGGCGGAATATCCCCGGTACTGCCGAGATCAATACTGCCAACGTTTAACGCTTCCAACATTTGCGGACCCGCGGGGAACTCCACCCAGGAGATTTTTGATTCCGGATAGCGTTTTTCCAGTAACTGGTGGCTTTTTGCCAGCACCATACCAATACTGCCTTTCTGATAGCCTATACGTAACGCTTCAGGCGAGGATTCTGCAGCAACCGCAAAAGTAGAGACGCTAAGCAATCCC  >  NC_000913/996621‑996915

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