Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 300,698 A→C 7.4% A79A (GCT→GCG)  paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,6980AC 7.4% 188.9 / 7.1 54A79A (GCT→GCG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base C (2/2);  total (24/30)
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.44e-01

GCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTC  >  NC_000913/300550‑300805
                                                                                                                                                     |                                                                                                           
gCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACa                                                                                                             <  4:892355/149‑1 (MQ=255)
 cGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAg                                                                                                            <  4:1225178/149‑1 (MQ=255)
        gccgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTAc                                                                                                     <  2:25485/149‑1 (MQ=255)
                      ttCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCttttt                                                                                       >  3:1343637/1‑149 (MQ=255)
                       tCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCG‑TACTGGCTGTCATTACCGCCAGTaacgcc                                                                                                 >  1:1269053/1‑132 (MQ=38)
                       tCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCG‑TACTGGCTGTCATTACCGCCAGTaacgcc                                                                                                 <  2:1269053/137‑6 (MQ=38)
                          aaGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgcgc                                                                                   <  3:253429/149‑1 (MQ=255)
                           aGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGcgcg                                                                                  >  3:352627/1‑149 (MQ=255)
                             gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCgcggc                                                                                <  3:777702/149‑1 (MQ=255)
                              tctcGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCgcggcg                                                                               >  1:244832/1‑149 (MQ=255)
                                    ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTa                                                                         <  4:352627/149‑1 (MQ=255)
                                     cACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGGCG‑TACTGGCTGTCATTACCGCCAGTaac                                                                                                    <  3:1007322/121‑3 (MQ=37)
                                     cACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGGCG‑TACTGGCTGTCATTACCGCCAGTaac                                                                                                    >  4:1007322/1‑119 (MQ=37)
                                     cACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTAt                                                                        <  4:1343637/149‑1 (MQ=255)
                                      aCTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGt                                                                          >  4:25622/1‑146 (MQ=255)
                                      aCTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGt                                                                          <  3:25622/146‑1 (MQ=255)
                                      aCTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATc                                                                       >  4:137048/1‑149 (MQ=255)
                                          ccAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAgg                                                                   <  1:766128/149‑1 (MQ=255)
                                             gcgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCttttt                                                                                       >  4:489563/1‑126 (MQ=255)
                                             gcgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCttttt                                                                                       <  3:489563/126‑1 (MQ=255)
                                                cAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTg                                                             >  3:1125742/1‑149 (MQ=255)
                                                       cTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtttt                                                                                        <  1:627266/115‑1 (MQ=255)
                                                       cTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtttt                                                                                        >  2:627266/1‑115 (MQ=255)
                                                       cTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTc                                                      >  4:139045/1‑149 (MQ=255)
                                                         gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCt                                                    <  3:139045/149‑1 (MQ=255)
                                                           tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCttt                                                  <  2:244832/149‑1 (MQ=255)
                                                            gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtttt                                                                                        >  1:1069895/1‑110 (MQ=255)
                                                            gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtttt                                                                                        <  2:1069895/110‑1 (MQ=255)
                                                            gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCttt                                                  <  4:746635/148‑1 (MQ=255)
                                                            gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCttt                                                  >  3:746635/1‑148 (MQ=255)
                                                            gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                 >  3:879339/1‑149 (MQ=255)
                                                                gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCaa                                             <  3:137048/149‑1 (MQ=255)
                                                                 tGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAAt                                            <  4:879339/149‑1 (MQ=255)
                                                                    tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                 <  3:537255/141‑1 (MQ=255)
                                                                    tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                 >  4:537255/1‑141 (MQ=255)
                                                                    tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGc                                         >  4:454613/1‑149 (MQ=255)
                                                                    tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGc                                         >  1:545312/1‑149 (MQ=255)
                                                                              ggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAg                                                                                                        <  1:633832/76‑1 (MQ=255)
                                                                              ggCCGCCTAACAGCCTGGCGGTGTTTTGGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAg                                                                                                        >  2:633832/1‑76 (MQ=255)
                                                                               gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgacga                              <  3:454613/149‑1 (MQ=255)
                                                                                   ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGAtata                          <  4:746540/149‑1 (MQ=255)
                                                                                            cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggc                 <  2:545312/149‑1 (MQ=255)
                                                                                                          gTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGctt   <  4:1125742/149‑1 (MQ=255)
                                                                                                               cTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGctt   <  2:33352/144‑1 (MQ=255)
                                                                                                               cTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGctt   >  1:33352/1‑144 (MQ=255)
                                                                                                                        gTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGata                            <  1:923936/110‑1 (MQ=255)
                                                                                                                        gTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGata                            >  2:923936/1‑110 (MQ=255)
                                                                                                                           cccc‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttc  >  1:615912/1‑133 (MQ=255)
                                                                                                                           cccc‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttc  <  2:615912/133‑1 (MQ=255)
                                                                                                                            ccc‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg       <  1:894144/127‑1 (MQ=255)
                                                                                                                            ccc‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg       >  2:894144/1‑127 (MQ=255)
                                                                                                                                        tGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAggg                                                                  <  1:644969/57‑1 (MQ=255)
                                                                                                                                        tGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAggg                                                                  >  2:644969/1‑57 (MQ=255)
                                                                                                                                        tGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                 >  4:1409943/1‑74 (MQ=255)
                                                                                                                                        tGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                 <  3:1409943/74‑1 (MQ=255)
                                                                                                                                                     |                                                                                                           
GCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCC‑GGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTC  >  NC_000913/300550‑300805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑