Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,433,237 C→T 16.7% W602* (TGG→TAG)  pfo ← pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,433,2370CT16.7% 50.7 / 4.7 24W602* (TGG→TAG) pfopyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (8/12);  new base T (2/2);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

ACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCG  >  CP009273/1433089‑1433378
                                                                                                                                                    |                                                                                                                                             
aCCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGcc                                                                                                                                               <  1:56007/149‑1 (MQ=255)
        cGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGt                                                                                                                                       >  1:137294/1‑149 (MQ=255)
                            aCGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCCGTTGCGTTCCACCAGATCCTGGCCctt                                                                                                                   >  1:81000/1‑149 (MQ=255)
                               gTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCTGAGCCTGCCAGTTGCGTGCCACCAGATCCTGGCCTTTGCt                                                                                                                >  1:48175/1‑149 (MQ=255)
                                ttttCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTCACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTa                                                                                                               <  1:174505/149‑1 (MQ=255)
                                    cACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGt                                                                                                           <  1:247592/149‑1 (MQ=255)
                                                             aCCACTGGCGGTCGGTTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATGCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCAcc                                                                                        <  2:103702/143‑1 (MQ=255)
                                                             aCCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCAcc                                                                                        <  2:277781/143‑1 (MQ=255)
                                                             aCCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCAcc                                                                                        <  2:228984/143‑1 (MQ=255)
                                                                 ccggCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGc                                                                                    <  2:246/141‑1 (MQ=255)
                                                                     cGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGa                                                                          <  1:94069/149‑1 (MQ=255)
                                                                                     tgCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATAGCCCCCTGCAATTCTGCCAGGACGCTATCCCCCgg                                                          >  1:120158/1‑149 (MQ=255)
                                                                                       cgatttgaCCGGTTGCAGCGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAAGCGCACCCTGCAATTCTGCGAGGGCGCTATCGc                                                              <  2:131775/136‑1 (MQ=255)
                                                                                           tttACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTACGAGGGCGCTATCGCCAGGCAGact                                                    >  1:224096/1‑147 (MQ=255)
                                                                                             tACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAAttt                                                  >  1:309099/1‑149 (MQ=255)
                                                                                              aCCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCCATTCTACGGGGGCGGTATCGCCCGGCAGGATTTg                                                 >  1:141032/1‑149 (MQ=255)
                                                                                              aCCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTg                                                 <  1:249836/149‑1 (MQ=255)
                                                                                                               ttcttcTACGGATTCACGCGCCAGCGCCAGAGCCTGCTAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGa                                      >  2:31641/1‑143 (MQ=255)
                                                                                                               ttcttcTACGGATTCACGCGCCAGCGCCAGAGCCTGCTAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGa                                      >  2:31691/1‑143 (MQ=255)
                                                                                                                 cttctACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCGGGCAGAATTTGCGTCAGATGGcaa                                    >  2:311615/1‑143 (MQ=255)
                                                                                                                                 cgccagcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGa                    <  2:255743/143‑1 (MQ=255)
                                                                                                                                  gccagcgccagAGCCTGCCAGTTGCGTTCCACCAGGTCCTGGCCTTTGCTACTGTAACTTTTGGCAAGCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGGCAGCTGctttac                   <  2:121504/143‑6 (MQ=255)
                                                                                                                                   ccagcgccagAGCCTGCTAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTa            <  1:35579/149‑1 (MQ=255)
                                                                                                                                   ccagcgccagAGCCTGCTAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTa            <  1:35527/149‑1 (MQ=255)
                                                                                                                                             agCCTGCCAGTTGCCTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGgccg  <  1:331090/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ACCAAGCCCGGCGAGCATCGCAGCGGTTACGGTTTTCACGAAATCAGGGGCGGCATCGGAAACCACTGGCGGTCGATTGGCGCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCG  >  CP009273/1433089‑1433378

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