Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,595,767 C→A A213S (GCG→TCG)  livH ← branched chain amino acid/phenylalanine ABC transporter membrane subunit LivH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,595,7670CA98.0% 171.5 / NA 49A213S (GCG→TCG) livHbranched chain amino acid/phenylalanine ABC transporter membrane subunit LivH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/24);  total (25/24)
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

GGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGCGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTCCGC  >  NC_000913/3595690‑3595842
                                                                             |                                                                           
ggCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGcc                                                                          <  1:628182/81‑1 (MQ=255)
   cATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAAt                                                                       >  1:628231/1‑81 (MQ=255)
    aTAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATc                                                                      <  2:628182/81‑1 (MQ=255)
      aaaGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATca                                                                     >  2:628229/1‑80 (MQ=255)
             aTGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTc                                                             <  2:628183/81‑1 (MQ=255)
                tAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgc                                                          >  1:628237/1‑81 (MQ=255)
                  ggggTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgcgc                                                        <  2:628185/81‑1 (MQ=255)
                  ggggTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgcgc                                                        <  2:628184/81‑1 (MQ=255)
                     gTTGATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAAt                                                     >  2:628237/1‑81 (MQ=255)
                     gTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAAt                                                     >  1:628235/1‑81 (MQ=255)
                      ttAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAAt                                                     >  1:628242/1‑80 (MQ=255)
                      ttAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATc                                                    <  1:628187/81‑1 (MQ=255)
                       tAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCa                                                   >  2:628256/1‑81 (MQ=255)
                           gACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcg                                               >  1:628256/1‑81 (MQ=255)
                             cGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggt                                             >  2:628253/1‑81 (MQ=255)
                              gCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtc                                            >  1:628250/1‑81 (MQ=255)
                               ccGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcg                                           <  2:628187/81‑1 (MQ=255)
                                cGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcgg                                          >  1:628243/1‑81 (MQ=255)
                                cGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcgg                                          >  1:628239/1‑81 (MQ=255)
                                 gTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcggt                                         >  1:628234/1‑81 (MQ=255)
                                  tAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGTGCAATCACCCGGTCGgtg                                        <  1:628189/81‑1 (MQ=255)
                                     aaCTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTa                                     <  2:628189/81‑1 (MQ=255)
                                        tGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATg                                  >  1:628240/1‑81 (MQ=255)
                                         gACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGc                                 >  2:628250/1‑81 (MQ=255)
                                         gACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGc                                 >  2:628240/1‑81 (MQ=255)
                                           ccGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCa                               <  1:628191/81‑1 (MQ=255)
                                            cGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCaa                              <  1:628192/81‑1 (MQ=255)
                                             gatcagcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAg                             >  1:628266/4‑81 (MQ=255)
                                              atcagcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAg                             >  1:628262/3‑80 (MQ=255)
                                               tcagcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCa                           <  2:628192/80‑1 (MQ=255)
                                                cagcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAg                             >  1:628247/1‑78 (MQ=255)
                                                 agcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAg                          <  1:628195/80‑1 (MQ=255)
                                                 agcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAg                          <  2:628193/80‑1 (MQ=255)
                                                 agcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGa                         >  1:628264/1‑81 (MQ=255)
                                                   cacaCCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACt                       <  2:628196/81‑1 (MQ=255)
                                                   cacaCCCCCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACt                       <  2:628195/81‑1 (MQ=255)
                                                     cacCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCg                     >  2:628258/1‑81 (MQ=255)
                                                       cccGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGcc                   <  1:628198/81‑1 (MQ=255)
                                                          gCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAtt                <  1:628199/81‑1 (MQ=255)
                                                           ccACCGCCGCCATCGCCGCGCCAATCACACAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAttt               >  2:628257/1‑81 (MQ=255)
                                                            cACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAtttt              >  2:628247/1‑81 (MQ=255)
                                                             aCCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTc             <  2:628201/81‑1 (MQ=255)
                                                             aCCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTc             <  2:628200/81‑1 (MQ=255)
                                                             aCCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTc             <  1:628200/81‑1 (MQ=255)
                                                                 ccgccATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGAt         <  1:628203/81‑1 (MQ=255)
                                                                  cgccATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATc        >  2:628270/1‑81 (MQ=255)
                                                                    ccATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCtt      <  2:628203/81‑1 (MQ=255)
                                                                     cATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTc     >  2:628259/1‑81 (MQ=255)
                                                                     cATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTc     <  1:628205/81‑1 (MQ=255)
                                                                         gccgaGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTCcgc  <  2:628205/80‑1 (MQ=255)
                                                                             |                                                                           
GGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGCGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTCCGC  >  NC_000913/3595690‑3595842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: