Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,595,057 A→C Y142D (TAC→GAC)  livM ← branched‑chain amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,595,0570AC98.2% 206.7 / ‑3.6 57Y142D (TAC→GAC) livMbranched‑chain amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (24/32);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACC  >  NC_000913/3594962‑3595148
                                                                                               |                                                                                           
tAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCcccc                                                                                         <  2:1315617/100‑1 (MQ=255)
tAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCcccc                                                                                         <  1:255358/100‑1 (MQ=255)
  aTCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccag                                                                                       >  2:1494603/1‑100 (MQ=255)
   tCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagc                                                                                      >  2:106656/1‑100 (MQ=255)
    ccAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagca                                                                                     <  2:1207414/100‑1 (MQ=255)
     cAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagcac                                                                                    >  2:1341120/1‑100 (MQ=255)
      aGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagc                                                                                      <  2:764304/97‑1 (MQ=255)
        cAATCGGCAGGCAGGTCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagcaccagc                                                                                <  2:423195/100‑1 (MQ=255)
         aaTCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagcaccagc                                                                                <  1:127253/100‑1 (MQ=255)
         aaTCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCccagcaccagc                                                                                >  1:97833/1‑100 (MQ=255)
           tCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACcagcag                                                                              <  2:636273/100‑1 (MQ=255)
             ggcaggcaGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGAc                                                                            <  2:423917/100‑1 (MQ=255)
                   caggTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGaaa                                                                      >  1:1278090/1‑100 (MQ=255)
                    aggTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAg                                                                     <  1:106656/100‑1 (MQ=255)
                    aggTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAg                                                                     >  1:1040497/1‑100 (MQ=255)
                        ccAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACca                                                                 >  2:883832/1‑100 (MQ=255)
                         cAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACcaa                                                                >  1:320147/1‑100 (MQ=255)
                         cAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACcaa                                                                >  2:1194425/1‑100 (MQ=255)
                               cccAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGGTGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAAcc                                                           <  1:1365029/99‑1 (MQ=255)
                                ccAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCAc                                                         >  1:1563733/1‑100 (MQ=255)
                                 cAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACg                                                        <  2:1040490/100‑1 (MQ=255)
                                  aaGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGt                                                       <  1:1277168/100‑1 (MQ=255)
                                    gCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTc                                                     <  2:1147680/100‑1 (MQ=255)
                                    gCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTc                                                     >  1:389233/1‑100 (MQ=255)
                                     ccGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCa                                                    >  1:1139432/1‑100 (MQ=255)
                                     ccGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCa                                                    >  1:1198246/1‑100 (MQ=255)
                                     acGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCa                                                    >  2:385547/2‑100 (MQ=255)
                                        tAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTCAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGcc                                                 <  2:1004283/100‑1 (MQ=255)
                                        tAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGcc                                                 >  2:1268543/1‑100 (MQ=255)
                                          aTAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCg                                               <  2:1139425/100‑1 (MQ=255)
                                           tAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCga                                              <  1:1164905/100‑1 (MQ=255)
                                             gTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCgaga                                            >  1:172237/1‑100 (MQ=255)
                                               gATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGAcc                                          >  2:74138/1‑100 (MQ=255)
                                               gATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGAcc                                          <  2:4071/100‑1 (MQ=255)
                                                aTTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCg                                         <  1:1092449/100‑1 (MQ=255)
                                                  tGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAg                                       <  2:1413250/100‑1 (MQ=255)
                                                    agcagcGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGa                                     >  2:425807/1‑100 (MQ=255)
                                                    agcagcGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGa                                     >  2:52842/1‑100 (MQ=255)
                                                      cagcGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATa                                   <  2:154751/100‑1 (MQ=255)
                                                          gcAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGt                               <  1:1313212/100‑1 (MQ=255)
                                                            aaaaGTGTAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAg                             <  2:325543/100‑1 (MQ=255)
                                                                gtgtAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCa                         <  2:654431/100‑1 (MQ=255)
                                                                 tgtAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCAt                        <  1:74138/100‑1 (MQ=255)
                                                                  gtAAGCGCCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATg                       <  1:1140888/100‑1 (MQ=255)
                                                                    aaGCGCTGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATgg                      <  2:320145/99‑1 (MQ=255)
                                                                      gcgcCGATGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCa                   >  2:584837/1‑100 (MQ=255)
                                                                             tGGCGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGc            <  2:1077042/100‑1 (MQ=255)
                                                                                cGTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAat         >  2:1422492/1‑100 (MQ=255)
                                                                                 gTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAata        >  1:1393252/1‑100 (MQ=255)
                                                                                 gTAAAAACCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAata        >  1:786808/1‑100 (MQ=255)
                                                                                   aaaaaCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATc      <  2:1393242/100‑1 (MQ=255)
                                                                                    aaaaCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATcc     >  1:318742/1‑100 (MQ=255)
                                                                                     aaaCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCa    >  2:21392/1‑100 (MQ=255)
                                                                                      aaCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCAc   <  1:52842/100‑1 (MQ=255)
                                                                                       aCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCAcc  <  2:1278081/100‑1 (MQ=255)
                                                                                       aCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCAcc  <  1:650663/100‑1 (MQ=255)
                                                                                       aCCGCCGTCCCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCAcc  <  2:8013/100‑1 (MQ=255)
                                                                                               |                                                                                           
TAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACC  >  NC_000913/3594962‑3595148

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