Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,113,851 G→A H208Y (CAT→TAT)  fpr ← flavodoxin/ferredoxin‑NADP(+) reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,8510GA95.0% 49.3 / NA 20H208Y (CAT→TAT) fprflavodoxin/ferredoxin‑NADP(+) reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/9);  minor base C (1/0);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACC  >  NC_000913/4113620‑4114091
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ggTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATc                                                                                                                                                                                                                               <  1:280942/250‑1 (MQ=255)
                             gATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATAGATAAGTAACCACTTCCCAGTAATGCTCAGCTGTCATCTGGCCCGGTCGGCGCCGTAAATGTTTAGTCATCTGACGGGTCTCTTTCAGCAACTGTTGTGTCTCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTcccc                                                                                                                                                                                                  >  1:335899/1‑251 (MQ=255)
                                                                                             gATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCttt                                                                                                                                                                                                            >  1:35891/1‑177 (MQ=255)
                                                                                                              cAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGttt                                                                                                                 >  1:22385/1‑251 (MQ=255)
                                                                                                                                      cccGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACt                                                                                                                                         >  1:160308/1‑203 (MQ=255)
                                                                                                                                              ncGACGTNNATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGtt                                                                                                                                                                                                                                           <  2:40065/96‑1 (MQ=255)
                                                                                                                                                    tAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTc                                                                                                                                               <  1:82237/183‑1 (MQ=255)
                                                                                                                                                     aaaTGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTCTGTATCGCGCACCATCTGAGGATCGCAGCACAGCATCACATAGGAGGTTTATTTATTCATCCGCAGGCCAATCGTGCTTTACAGTTCCCCACTTTCCATTAATGCCGGTATCCGTCCGGTGAGCGACCACGCTGCCGTTTCCCGAc                                                                                                                                          >  1:283590/1‑187 (MQ=255)
                                                                                                                                                     aaaTGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGATGTGTATCGCGAACCATCTGTGGCTTGACGCACAGCATCACATCGCTGGTTTCTTTATTCATAGGCAGGCCAATCGTGCTGTCCACGTCCCCACTTTCAATTAATGCCGGTATCCGTTCGGTGAGCCCCCCCACTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCtt                                                                                                            >  1:224276/1‑217 (MQ=255)
                                                                                                                                                                                  cAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAGCCGGCCGGTGAGCGACACCGCTGCCGTGTCCCGACTGACCACCGTCTGAATGCGCAGGTTGCCGTCGGAGCGGTTTGCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTgcg                                              <  1:370553/250‑1 (MQ=255)
                                                                                                                                                                                      aaCTGTTGTGTATCGCTCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGGGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGAcc                                                                                                                                                            <  1:334261/136‑1 (MQ=255)
                                                                                                                                                                                         cgatccgtATCGCGCACCATCTGTGGATGGCCGCACAGCAGCACATAGCTGGGTTCTTTATTCATCGGCAGTCCAATCGTGCTCTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCTACCCCTCTGCTGTTGCCCGACTTACCACCGTCTTAATGCGCAGTTTTCCTTCGGAGCCTGTTTCCAGTTCCCGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCAc                                      <  1:364606/245‑1 (MQ=255)
                                                                                                                                                                                          gttgtGTAGCGCTCACCATCTGTGGATTGCCGCACATCATCACATAGCTGGTTTCTTTATTCATCTGCAGGCCAATCTTGCTTTCCAGTTCCCCACTGTCAAATAATGCCGGTATCCTGCCGGTGATCGACCCCGCGTCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTcc                                                                                                              <  1:337070/178‑1 (MQ=255)
                                                                                                                                                                                                    gcgcACCATCTGTGGATTGCCGCACAGCATCACCTAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCtt                                                                                                                                                                                                             >  1:159751/1‑73 (MQ=255)
                                                                                                                                                                                                                nggATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTTATCGGCAGGCCACCCGGGCTTTCCAGTTCCCCAATTTAACTTAATGCCGGTATCCATCCGGTGGGCGACCCCCCTGCCGTTTCCCGACTGACCCCAGTCTCAATGCGCCGTTTTCCTTCGTAGAGATTTTCCAGTTCATGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTg                                                >  1:59418/2‑218 (MQ=255)
                                                                                                                                                                                                                       ccGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCGCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCTACCCCGCTGCCGTTGCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTGCGTAGCGTTTTTc                                                                                               <  1:322755/164‑1 (MQ=255)
                                                                                                                                                                                                                       ccGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCCGTTTTCCTTCGTAGCGTTTTTCCAGtt                                                                                          >  1:323465/1‑169 (MQ=255)
                                                                                                                                                                                                                          cacaGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTACCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCttt     >  1:433941/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCATCACATAGCTGTTTTCTTTATTCATCGGCAGGCCAATCGTGCTGTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGGTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTGTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACTTGCGGCTTGCACCAGGCCCAGATTTTTGAAGCGATCTAAACCTTTatc  <  1:297053/251‑3 (MQ=255)
                                                                                                                                                                                                                                 tCACATAGCTGGTTTCTTTATTCATCGGCCGTCCAATCGCGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCAGTCCGGTGAGCCACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGACTGCGCAGTTTTCCTTCGTAGCGTTTTTc                                                                                               >  1:220658/1‑154 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACC  >  NC_000913/4113620‑4114091

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