Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,035,337 T→A 6.8% L91Q (CTG→CAG)  hyaC → hydrogenase 1, b‑type cytochrome subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,035,3370TA 6.8% 172.1 / 3.7 59L91Q (CTG→CAG) hyaChydrogenase 1, b‑type cytochrome subunit
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CGCTACCTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGAT  >  NC_000913/1035190‑1035481
                                                                                                                                                   |                                                                                                                                                
cGCTACCTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTg                                                                                                                                                 >  3:278058/1‑149 (MQ=255)
 gCTACCTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGt                                                                                                                                                >  5:294046/1‑149 (MQ=255)
  cTACCTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGtt                                                                                                                                               <  2:335693/149‑1 (MQ=255)
     ccTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTCCATCAGGTTACTTCACTTCAGCGCGGGGCTGGTTTTTACCGTGGTTTTGCTGAGGCGGATCTACTGGGCTGTTGTTGGCAATCGATACTCCCGCGAGCTGttt                                                                                                                                              <  4:299122/147‑1 (MQ=255)
      cTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATc                                                                                                                                           <  6:2722/149‑1 (MQ=255)
          cGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGc                                                                                                                                       <  6:279676/149‑1 (MQ=255)
               gcggcgAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTa                                                                                                                                  >  2:165900/1‑149 (MQ=255)
                   cgaggcgaCGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGc                                                                                                                              <  1:165900/149‑1 (MQ=255)
                         gaCGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTaaaa                                                                                                                        <  7:54511/149‑1 (MQ=255)
                                   ttCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCtggtgg                                                                                                                >  4:92732/1‑147 (MQ=255)
                                       atatGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGc                                                                                                          >  8:255915/1‑149 (MQ=255)
                                        tatGTGCTACTTCAGGTTATTTCACTTCTGCGCCGGGATGTTTTTTACCGTGGTTTTGCTGATGAGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCTCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTTGTGGCAGGGc                                                                                                          <  4:278058/148‑1 (MQ=255)
                                             gCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTgg                                                                                                    <  6:294046/149‑1 (MQ=255)
                                                 cATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGt                                                                                                                            <  1:367096/121‑1 (MQ=255)
                                                 cATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGt                                                                                                                            >  2:367096/1‑121 (MQ=255)
                                                   acAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGATTTGCTGATGCTGATCTACTGCGCTTTTGTTCACAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGTCAGGTCGTGTGGTATGa                                                                                               >  5:2745/2‑148 (MQ=255)
                                                      ggTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATc                                                                                           >  4:291122/1‑149 (MQ=255)
                                                         tAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGc                                                                                        >  3:253396/1‑149 (MQ=255)
                                                            ttaaCTTCAGCGCTTGTTTGGTTTTTACAGTTGTTTTTCTGATGCTGATCTACTGGGCTTTTGTTGGCAATCGATACTCGCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGGAGGGCGTGTGGTATGAAATCCGCTg                                                                                      <  6:2745/145‑1 (MQ=255)
                                                            ttCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGc                                                                                                                              >  8:76924/1‑108 (MQ=255)
                                                            ttCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGc                                                                                                                              <  7:76924/108‑1 (MQ=255)
                                                             tCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGc                                                                                        <  2:242627/145‑1 (MQ=255)
                                                             tCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGc                                                                                        >  1:242627/1‑145 (MQ=255)
                                                               aCTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGa                                                                                               <  6:214716/136‑1 (MQ=255)
                                                               aCTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGa                                                                                               >  5:214716/1‑136 (MQ=255)
                                                                             aTGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCAGTTTATCGTGCCGGTATGGCGTAAAAGCtg                                                                                                                    >  3:110014/1‑101 (MQ=255)
                                                                             aTGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCAGTTTATCGTGCCGGTATGGCGTAAAAGCtg                                                                                                                    >  3:105123/1‑101 (MQ=255)
                                                                             aTGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCAGTTTATCGTGCCGGTATGGCGTAAAAGCtg                                                                                                                    <  4:105123/101‑1 (MQ=255)
                                                                             aTGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCAGTTTATCGTGCCGGTATGGCGTAAAAGCtg                                                                                                                    <  4:110014/101‑1 (MQ=255)
                                                                                      aCCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTAt                                                                                                 <  6:135181/111‑1 (MQ=255)
                                                                                      aCCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTAt                                                                                                 >  5:135181/1‑111 (MQ=255)
                                                                                      aCCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAg                                                           >  7:203359/1‑149 (MQ=255)
                                                                                        cGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATc                                                                                           <  1:117007/115‑1 (MQ=255)
                                                                                        cGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATc                                                                                           >  2:117007/1‑115 (MQ=255)
                                                                                             ttttGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGTTATGAAATCAGCTGGTATCTGTTTCTGTCAAAACGTCCGAGTGCCGa                                                     >  3:13660/1‑148 (MQ=255)
                                                                                                     aTGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAgg                                                                                                            <  1:307116/85‑1 (MQ=255)
                                                                                                     aTGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAgg                                                                                                            >  2:307116/1‑85 (MQ=255)
                                                                                                         ggATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTTCCGATATAGGCCATAAt                                        >  5:327036/1‑149 (MQ=255)
                                                                                                          gATCTACTGGGCTTTTGTTGGCAATCGATAATCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTTCCGATATAGGCCATAATc                                       <  6:327036/149‑1 (MQ=255)
                                                                                                              tACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCAt                                   >  3:276747/1‑149 (MQ=255)
                                                                                                                 tGGGCTTTTGATGGCAATCGATACTCCCGCGAGCTGTATATCGTGCCGGAATGGCCTAAAAGCTGGTGGCAGGGCCAGTGGTATGAAATCAGCTGGTATCTCTTTCTGGCAAAACGTCCGAGTGACGAAATAGGCCATAATCCCATCGc                                >  2:11219/1‑149 (MQ=255)
                                                                                                                   ggCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGcc                               >  1:379293/1‑148 (MQ=255)
                                                                                                                      ttttGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAgg                           >  7:22238/1‑149 (MQ=255)
                                                                                                                        ttgttgGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAgg                           <  4:253396/147‑1 (MQ=255)
                                                                                                                          gttgGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAgg                           <  4:233501/145‑1 (MQ=255)
                                                                                                                          gttgGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAgg                           >  3:233501/1‑145 (MQ=255)
                                                                                                                               cAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGt                  <  8:203359/149‑1 (MQ=255)
                                                                                                                                aaTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTc                                                                           <  6:109187/91‑1 (MQ=255)
                                                                                                                                aaTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTc                                                                           >  5:109187/1‑91 (MQ=255)
                                                                                                                                aaTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGtt                 >  5:218432/1‑149 (MQ=255)
                                                                                                                                 aTCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTc                <  8:22238/149‑1 (MQ=255)
                                                                                                                                    gATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGc             >  8:364191/1‑149 (MQ=255)
                                                                                                                                        cTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTAt          <  7:255915/148‑1 (MQ=255)
                                                                                                                                           ccGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTcc      >  1:223717/1‑149 (MQ=255)
                                                                                                                                           ccGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTcc      >  2:201298/1‑149 (MQ=255)
                                                                                                                                           ccGAGAGCTGTTTATAGTTCAGGTATGTCGTAAATGCTGGTGGCAGGGCGTGTGGTATGAAATCCGATGGTATCTGTTTCTTGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCATTAGGCGATGTTCGGCTATTTcc      <  4:13660/149‑1 (MQ=255)
                                                                                                                                            cgcgAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCTGGCGGCGATGTTCGGCTATTTCCt     <  1:11219/149‑1 (MQ=255)
                                                                                                                                            cgcgAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGGGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCt     <  3:92732/149‑1 (MQ=255)
                                                                                                                                               gAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGAt  <  3:291122/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGCTACCTTCCGTCAGCGGCGAGGCGACGTATCTGTTCTATATGGGCTACATCAGGTTAATTCACTTCAGCGCCGGGATGGTTTTTACCGTGGTTTTGCTGATGCGGATCTACTGGGCTTTTGTTGGCAATCGATACTCCCGCGAGCTGTTTATCGTGCCGGTATGGCGTAAAAGCTGGTGGCAGGGCGTGTGGTATGAAATCCGCTGGTATCTGTTTCTGGCAAAACGTCCGAGTGCCGATATAGGCCATAATCCCATCGCCCAGGCGGCGATGTTCGGCTATTTCCTGAT  >  NC_000913/1035190‑1035481

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