Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 847,926 T→C 14.8% A289A (GCT→GCC)  trxB1 → thioredoxin‑disulfide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730847,9260TC14.8% 60.0 / 6.0 27A289A (GCT→GCC) trxB1thioredoxin‑disulfide reductase
Reads supporting (aligned to +/- strand):  ref base T (8/15);  new base C (2/2);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

ACAGCGCCATTTAAAGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGATATTTATTTTAATAGTGTCATCACAG  >  CP000730/847789‑848062
                                                                                                                                         |                                                                                                                                        
aCAGCGCCATTTAAAGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACt                                                                                                                                       <  2:181208/141‑1 (MQ=255)
        aTTTAAAGACTTAGGTATTACAAATGATGTTGGTTATATTGTCACCAAAGATGATATGACAACATCAGTCCCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATgg                                                                                                                               <  1:92239/141‑1 (MQ=255)
        aTTTAAAGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATgg                                                                                                                               <  1:26465/141‑1 (MQ=255)
         tttAAAGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGgt                                                                                                                              >  2:77099/1‑141 (MQ=255)
              aGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTTCATCAGTTGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCCACTGGCGATGGTAGTAt                                                                                                                         >  2:178375/1‑141 (MQ=255)
              aGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCCACTGGCGATGGTAGTAt                                                                                                                         >  1:34550/1‑141 (MQ=255)
                   tAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTTCGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAg                                                                                                                    <  2:227882/141‑1 (MQ=255)
                      gTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAgcgc                                                                                                                 <  1:152447/141‑1 (MQ=255)
                           aCAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTTTTTTTTTTTTTTTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGt                                                                                                            >  2:40715/1‑141 (MQ=255)
                           aCAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTTTATTTTGTTTTGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGt                                                                                                            >  2:247673/1‑141 (MQ=255)
                           aCAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGt                                                                                                            <  2:246248/141‑1 (MQ=255)
                            cAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTg                                                                                                           <  2:159898/141‑1 (MQ=255)
                            cAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTg                                                                                                           <  1:213474/141‑1 (MQ=255)
                            cAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCACAAAGTg                                                                                                           <  1:187795/141‑1 (MQ=255)
                            aaaaTGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTg                                                                                                           <  2:78301/140‑1 (MQ=255)
                             aaaTGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGc                                                                                                          <  2:142972/141‑1 (MQ=255)
                                                          atgatATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATc                                                                             >  1:117726/1‑141 (MQ=255)
                                                                                                           cgACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGa                            >  2:58255/1‑141 (MQ=255)
                                                                                                            gACAAAGGTTTACGCCAAATTGTCACTGCCACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGat                           <  1:178375/141‑1 (MQ=255)
                                                                                                            gACAAAGGTTTACGCCAAATTGTCACTGCCACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGat                           <  2:34550/141‑1 (MQ=255)
                                                                                                              cAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGatat                         >  1:181871/1‑141 (MQ=255)
                                                                                                              cAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGatat                         >  2:253100/1‑141 (MQ=255)
                                                                                                                   gTTTACGCCACATTGTCACCGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGATCATTTACACGATCAAGCTTAATGCGAAGCCGAATTAAGATGTTGAGCTGTAAAGTATTTGGAtatttatt                    <  2:103976/141‑1 (MQ=255)
                                                                                                                      tACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGATATTTATTTTa                 <  1:63862/141‑1 (MQ=255)
                                                                                                                            aaaTTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGATATTTATTTTAATAgtg           <  1:163655/141‑1 (MQ=255)
                                                                                                                                    aCTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGATATTTATTTTAATAGTGTCATcaca   >  1:82331/1‑141 (MQ=255)
                                                                                                                                     cTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGATATTTATTTTAATAGTGTCATCACAg  <  2:168343/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
ACAGCGCCATTTAAAGACTTAGGTATTACAAATGATGTTGGTTATATTGTAACAAAAGATGATATGACAACATCAGTACCAGGTATTTTTGCAGCAGGAGATGTTCGCGACAAAGGTTTACGCCAAATTGTCACTGCTACTGGCGATGGTAGTATTGCAGCGCAAAGTGCAGCGGAATATATTGAACATTTAAACGATCAAGCTTAATTCGAAGTCGAATTAAGATGTTGAGCTGTAAATTATTTGGATATTTATTTTAATAGTGTCATCACAG  >  CP000730/847789‑848062

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