Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,506,757 C→G 8.6% F151L (TTC→TTG)  gqr → glutathionyl‑hydroquinone reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,506,7570CG 8.6% 171.7 / 4.8 56F151L (TTC→TTG) gqrglutathionyl‑hydroquinone reductase
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (3/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GTACCTGCACCAGCGCTACACTCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGCAAGACGCT  >  NC_002947/4506619‑4506892
                                                                                                                                          |                                                                                                                                        
gTACCTGCACCAGCGCTACACTCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGa                                                                                                                                >  8:128580/1‑149 (MQ=255)
            gcgcTACACTCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCgg                                                                                                                            <  1:279482/141‑1 (MQ=255)
            gcgcTACACTCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCgg                                                                                                                            >  2:279482/1‑141 (MQ=255)
              gcTACACTCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTg                                                                                                                  <  7:128580/149‑1 (MQ=255)
                     tCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATAcc                                                                                                                     >  1:152133/1‑139 (MQ=255)
                     tCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATAcc                                                                                                                     <  2:152133/139‑1 (MQ=255)
                       aGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTc                                                                                                                                          <  6:148562/116‑1 (MQ=255)
                       aGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTc                                                                                                                                          >  5:148562/1‑116 (MQ=255)
                              cccGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGaa                                                                                                                               >  1:256873/1‑119 (MQ=255)
                              cccGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGAc                                                                                                                               <  2:256873/120‑1 (MQ=255)
                               ccGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATc                                                                                                 >  3:23130/1‑149 (MQ=255)
                                        acaGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCaa                                                                                                                                        <  2:12282/101‑1 (MQ=255)
                                        acaGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCaa                                                                                                                                        >  1:12282/1‑101 (MQ=255)
                                             gcgcGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCa                                                                                   >  4:43192/1‑149 (MQ=255)
                                              cgcgTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCAt                                                                                  <  4:23130/149‑1 (MQ=255)
                                                 gTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCg                                                                                >  2:294290/1‑148 (MQ=255)
                                                      ggTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCTCGTTGAACGCGCTgt                                                                                                                                <  1:289858/95‑2 (MQ=255)
                                                      ggTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCTCGTTGAACGCGCTgt                                                                                                                                >  2:289858/1‑94 (MQ=255)
                                                      ggTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGAGCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCAc                                                                          <  8:1134/149‑1 (MQ=255)
                                                      ggTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCAc                                                                          >  7:1134/1‑149 (MQ=255)
                                                            tGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAAc                                                                     >  3:215459/1‑148 (MQ=255)
                                                                   tgGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAgc                                                                 >  6:116813/1‑145 (MQ=255)
                                                                   tgGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAgc                                                                 <  5:116813/145‑1 (MQ=255)
                                                                           gaaagaGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAgcgc                                                               <  7:129618/139‑1 (MQ=255)
                                                                           gaaagaGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAgcgc                                                               >  8:129618/1‑139 (MQ=255)
                                                                           gaaagaGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCcg‑g                                                    >  1:326355/1‑149 (MQ=255)
                                                                                agaAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGt                                                >  7:169771/1‑148 (MQ=255)
                                                                                 gaAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGt                                                >  2:234406/1‑147 (MQ=255)
                                                                                 gaAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGt                                                <  1:234406/147‑1 (MQ=255)
                                                                                            tCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCAt                                                                                  <  6:253735/103‑1 (MQ=255)
                                                                                            tCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCAt                                                                                  >  5:253735/1‑103 (MQ=255)
                                                                                              aacaaTGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTgcgc                                                                                          >  8:12299/1‑93 (MQ=255)
                                                                                              aacaaTGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTgcgc                                                                                          <  7:12299/93‑1 (MQ=255)
                                                                                              aacaaTGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGa                                                                   <  8:113010/116‑1 (MQ=255)
                                                                                              aacaaTGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGa                                                                   >  7:113010/1‑116 (MQ=255)
                                                                                                  aTGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTg                                               >  5:5866/1‑131 (MQ=255)
                                                                                                  aTGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTg                                               <  6:5866/131‑1 (MQ=255)
                                                                                                      gTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGCCCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAg                         >  2:243275/1‑149 (MQ=255)
                                                                                                       tCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGCCATTTACCCCGCATCGCTGCGCCCGACAATCAAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGAAgg                        >  5:26914/1‑149 (MQ=255)
                                                                                                              aaaTCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCt                 <  4:215459/149‑1 (MQ=255)
                                                                                                               aaTCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCt                 <  3:43192/148‑1 (MQ=255)
                                                                                                                   atcCGCATTATCAACAGCGCGTTGAACGATCTGAACGGCAATACCCTTGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCa             <  6:26914/148‑1 (MQ=255)
                                                                                                                   atcCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCAc            <  8:169771/149‑1 (MQ=255)
                                                                                                                   atcCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCAc            <  2:326355/149‑1 (MQ=255)
                                                                                                                     ccGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACg           >  6:100171/1‑148 (MQ=255)
                                                                                                                      ctcATATTGAACAGCGCGGTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGGTGGGCCCGTCCATCGAAGCACTGAACGCGCGCATCTACCCGCGCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACg           <  1:243275/146‑1 (MQ=255)
                                                                                                                       gCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCttt                    <  2:19547/137‑1 (MQ=255)
                                                                                                                       gCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCttt                    >  1:19547/1‑137 (MQ=255)
                                                                                                                       gCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGc          <  6:94487/147‑1 (MQ=255)
                                                                                                                       gCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGc          >  5:94487/1‑147 (MQ=255)
                                                                                                                       gCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCACTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGc          >  1:99810/1‑147 (MQ=255)
                                                                                                                          tatTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGCAAGAc     <  5:100171/149‑1 (MQ=255)
                                                                                                                            ccggtcAGCTCGTTGAACGCGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGCAAGACGc   >  6:170936/6‑149 (MQ=255)
                                                                                                                              ggtcAGCTCGTTGAACGCGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGCAAGACGCt  >  3:333845/4‑148 (MQ=255)
                                                                                                                               aaCAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGc                  <  4:268459/131‑1 (MQ=255)
                                                                                                                               aaCAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGc                  >  3:268459/1‑131 (MQ=255)
                                                                                                                                 cAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCAc            >  7:155416/1‑135 (MQ=255)
                                                                                                                                 cAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCAc            <  8:155416/135‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GTACCTGCACCAGCGCTACACTCAGGATGACCCGCATTACACAGGGCGCGTGACGGTGCCTGTACTGTGGGACAAGAAAGAGAAGCGCATCGTCAACAATGAGTCGGCAGAAATCATCCGCATATTCAACAGCGCGTTCAACGAGCTGACCGGCAATACCCTGGACCTTTACCCCGCATCGCTGCGCCCGACCATCGAAGCACTGAACGAGCGCATCTACCCG‑GCGGTGAACAACGGTGTGTACCGCGCAGGCTTTGCTACCACGCAAGACGCT  >  NC_002947/4506619‑4506892

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