Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,142,898 G→A A31V (GCA→GTA)  hybC ← hydrogenase 2 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,142,8980GA97.4% 136.0 / ‑3.6 38A31V (GCA→GTA) hybChydrogenase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/26);  minor base T (0/1);  total (11/27)
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.87e-01

ACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAA  >  NC_000913/3142820‑3142970
                                                                              |                                                                        
acaATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACt                                                                        <  1:588429/81‑1 (MQ=255)
   atcatcCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTc                                                                     <  1:588430/81‑1 (MQ=255)
   atcatcCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTc                                                                     <  2:588431/81‑1 (MQ=255)
      atcCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGa                                                                   <  1:588433/80‑1 (MQ=255)
      atcCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGa                                                                   <  2:588432/80‑1 (MQ=255)
          aTGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAacga                                                              <  2:588434/81‑1 (MQ=255)
            gCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAacgacg                                                            >  1:588496/1‑81 (MQ=255)
              aTCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGcc                                                          <  1:588435/81‑1 (MQ=255)
              aTCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGc                                                           >  2:588487/1‑80 (MQ=255)
                 gcgcgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCAtt                                                       <  2:588437/81‑1 (MQ=255)
                 gcgcgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCAtt                                                       <  2:588436/81‑1 (MQ=255)
                    cgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTc                                                    <  2:588438/81‑1 (MQ=255)
                    cgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTc                                                    <  1:588440/81‑1 (MQ=255)
                    cgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTc                                                    <  1:588439/81‑1 (MQ=255)
                       aTCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGAt                                                 <  1:588441/81‑1 (MQ=255)
                           cgGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCg                                             <  2:588442/81‑1 (MQ=255)
                            ggTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATTCTTTCGAAACGACGCCATTTTCGATTTCGc                                            <  1:588443/81‑1 (MQ=255)
                             gTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGc                                            >  1:588507/1‑80 (MQ=255)
                                 ttCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATc                                        >  2:588489/1‑80 (MQ=255)
                                  tCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCg                                       <  1:588444/80‑1 (MQ=255)
                                  tCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGa                                      >  1:588493/1‑81 (MQ=255)
                                       aTCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCg                                  >  1:588494/1‑80 (MQ=255)
                                       aTCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCg                                  >  1:588513/1‑80 (MQ=255)
                                            ttCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGt                                 <  2:588447/76‑1 (MQ=255)
                                            ttCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGt                                 <  1:588445/76‑1 (MQ=255)
                                            ttCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTaaa                              <  2:588446/79‑1 (MQ=255)
                                            ttCCATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTa                                <  2:588448/77‑1 (MQ=255)
                                              ccATGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGcc                          >  2:588501/1‑81 (MQ=255)
                                                aTGCCGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGcccc                        <  1:588449/81‑1 (MQ=255)
                                                   ccGCGCCACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTc                      <  2:588450/80‑1 (MQ=255)
                                                        ccACATGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACg                <  2:588451/81‑1 (MQ=255)
                                                           caTGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGt             <  1:588452/81‑1 (MQ=255)
                                                           caTGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGt             <  1:588453/81‑1 (MQ=255)
                                                            aTGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGtt            <  1:588454/81‑1 (MQ=255)
                                                             tGGTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTa           <  2:588455/81‑1 (MQ=255)
                                                               gTACCGGAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTAcc         >  2:588511/1‑81 (MQ=255)
                                                                    ggAAGCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGAt     >  1:588504/1‑80 (MQ=255)
                                                                        gCCCATACTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCaa  >  1:588456/1‑79 (MQ=255)
                                                                              |                                                                        
ACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAA  >  NC_000913/3142820‑3142970

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