Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 870,807 C→T A466V (GCG→GTG)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913870,8070CT97.4% 127.7 / ‑4.2 38A466V (GCG→GTG) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/15);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGGCACTATCGCCGCTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAAC  >  NC_000913/870663‑870954
                                                                                                                                                |                                                                                                                                                   
gggCACTATCGCCGCTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGgcgc                                                                                                                                                >  2:324304/1‑150 (MQ=255)
 ggCACTATCGCCGCTGTTTGCCTCGCAGAACTGGCCCCCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAg                                                                                                                                              <  1:70817/151‑1 (MQ=255)
  gCACTATCGCCGCTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAgg                                                                                                                                             <  1:390321/151‑1 (MQ=255)
    aCTATCGCCGCTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGat                                                                                                                                           >  1:669901/1‑151 (MQ=255)
         cgccgcTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATatcat                                                                                                                                      <  2:619370/151‑1 (MQ=255)
         cgccgcTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATatca                                                                                                                                       <  1:374264/150‑1 (MQ=255)
            cgcTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTg                                                                                                                                   <  1:340206/151‑1 (MQ=255)
                   tGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGa                                                                                                                             >  2:385782/1‑150 (MQ=255)
                                tGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATcc                                                                                                               >  1:630926/1‑151 (MQ=255)
                                             tGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTggtggt                                                                                                     >  1:556656/1‑148 (MQ=255)
                                                  aaTACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAAc                                                                                             >  2:598102/1‑151 (MQ=255)
                                                           ttttACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCgg                                                                                    >  1:208889/1‑151 (MQ=255)
                                                           ttttACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCgg                                                                                    >  2:437981/1‑151 (MQ=255)
                                                              tACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGca                                                                                  <  1:9846/150‑1 (MQ=255)
                                                                        aaCAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTaaaa                                                                        >  2:221598/1‑150 (MQ=255)
                                                                           aGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCt                                                                    <  2:669821/151‑1 (MQ=255)
                                                                               ggATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGAcc                                                                >  2:499481/1‑151 (MQ=255)
                                                                                    aCTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttt                                                           >  1:151910/1‑151 (MQ=255)
                                                                                      ttCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                         >  1:657427/1‑151 (MQ=255)
                                                                                      ttCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                         >  2:518453/1‑151 (MQ=255)
                                                                                       tCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                         <  1:536572/150‑1 (MQ=255)
                                                                                        ccTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                       <  2:198536/151‑1 (MQ=255)
                                                                                               cAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGc                                                >  2:480350/1‑151 (MQ=255)
                                                                                               cAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCAt                                                  <  2:530861/149‑1 (MQ=255)
                                                                                                aGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGcc                                               >  2:177834/1‑151 (MQ=255)
                                                                                                  gCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGcc                                               >  2:328371/1‑149 (MQ=255)
                                                                                                         aaaCTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCgg                                      >  2:146463/1‑151 (MQ=255)
                                                                                                           aCTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCt                                    >  1:120816/1‑151 (MQ=255)
                                                                                                                    ccGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACGCCGGCTTCAGCttt                            <  1:635797/150‑1 (MQ=255)
                                                                                                                    ccGGCGGAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCttt                            <  2:444248/150‑1 (MQ=255)
                                                                                                                      ggcggAAAAGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGa                          >  2:567117/1‑150 (MQ=255)
                                                                                                                              aGACCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCgg                  >  2:510009/1‑150 (MQ=255)
                                                                                                                                aCCCGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAtt               >  1:54678/1‑151 (MQ=255)
                                                                                                                                 cccGCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAttt              <  1:71196/151‑1 (MQ=255)
                                                                                                                                    gCTTATATAAGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCTGCTTCAGCTTTGAAGACGCGGATTTGCa           <  1:163296/151‑1 (MQ=255)
                                                                                                                                            aaGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCaa   <  1:408530/151‑1 (MQ=255)
                                                                                                                                            aaGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCCGATTTGCAATa        >  1:289298/1‑146 (MQ=255)
                                                                                                                                             aGGTGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAAc  <  2:128752/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGGCACTATCGCCGCTGTTTGCCTCGCAGAACTGGCCACCGACGCTGTTTAATACCGCGTTTTACAGCAATAAACAGGTGGATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAAC  >  NC_000913/870663‑870954

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