Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,142,568 T→C H141R (CAC→CGC)  hybC ← hydrogenase 2 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,142,5680TC97.3% 122.9 / NA 37H141R (CAC→CGC) hybChydrogenase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/19);  minor base A (0/1);  total (17/20)
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.98e-01

CTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATGGTCATGCGTGGTGTGC  >  NC_000913/3142443‑3142696
                                                                                                                             |                                                                                                                                
cTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACg                                                                                                                    <  2:347956/140‑1 (MQ=255)
 ttCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGc                                                                                                                   >  1:568114/1‑140 (MQ=255)
   ccGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCt                                                                                                                 >  1:646603/1‑140 (MQ=255)
      gcggcAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTc                                                                                                              >  1:461940/1‑140 (MQ=255)
                cATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCgg                                                                                                    <  1:613980/140‑1 (MQ=255)
                     ccGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCt                                                                                               <  1:651520/140‑1 (MQ=255)
                              cccAGTAGCCATTAGCGAAAATCCCAAACTGACCGCTAGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTggg                                                                                      <  1:264396/140‑1 (MQ=255)
                                   tAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGc                                                                                 >  2:155525/1‑140 (MQ=255)
                                      ccATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCATCTtg                                                                              <  1:126414/140‑1 (MQ=255)
                                      ccATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTtg                                                                              <  1:475390/140‑1 (MQ=255)
                                           aGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTtgcagtg                                                                         <  1:356453/140‑1 (MQ=255)
                                             cGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTtgcagtgca                                                                       <  2:528726/140‑1 (MQ=255)
                                             cGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTtgcagtgca                                                                       >  1:6082/1‑140 (MQ=255)
                                                  aTACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGt                                                                  <  1:471199/140‑1 (MQ=255)
                                                          cTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCa                                                          >  2:375495/1‑140 (MQ=255)
                                                            gACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCAcc                                                        <  2:6082/140‑1 (MQ=255)
                                                             aCCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCATGAAGTGCAGAAGTGATGTCCAccc                                                       >  1:228025/1‑140 (MQ=255)
                                                               cGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAg                                                     >  1:374709/1‑140 (MQ=255)
                                                                          aGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGa                                          <  2:461940/140‑1 (MQ=255)
                                                                                    tCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAg                                <  2:646603/140‑1 (MQ=255)
                                                                                     cTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGa                               <  2:66813/140‑1 (MQ=255)
                                                                                           cTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGaa                         <  2:500847/140‑1 (MQ=255)
                                                                                              aaCTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACaa                      >  2:552645/1‑140 (MQ=255)
                                                                                               aCTTTGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAat                     <  2:553119/140‑1 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCTGTCCAGCGCCGGAAGCTGATACAAATGAACAATATgg                 >  1:215429/1‑140 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGaa                         <  1:101608/132‑1 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGaa                         >  2:101608/1‑132 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATgg                 >  2:413426/1‑140 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATgg                 >  2:436474/1‑140 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATgg                 <  1:409280/140‑1 (MQ=255)
                                                                                                   tGGTGAACTCTTCCGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATgg                 >  1:401213/1‑140 (MQ=255)
                                                                                                    ggTGAACTCTTCCGGCCTGTTCTGGAGCCAGGTCGAAGCGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATGGt                <  1:299594/140‑1 (MQ=255)
                                                                                                               ccGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGTAAGCTGAAAGAGATAAACAATATGGTCATGCgtggtg     >  1:242029/1‑140 (MQ=255)
                                                                                                               ccGGACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGCAAGCTGATAGACATGAACAATATGGTCATGCgtggtg     >  1:266756/1‑140 (MQ=255)
                                                                                                                  gACTGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATGGTCATGCGTGGTGTgc  <  2:374709/140‑1 (MQ=255)
                                                                                                                     tGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATGGt                <  2:432339/123‑1 (MQ=255)
                                                                                                                     tGTTCAGGCGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATGGt                >  1:432339/1‑123 (MQ=255)
                                                                                                                             |                                                                                                                                
CTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGCAGAAGTGATGTCCACCCAGTCCAGCGCCGAAAGCTGATAGAAATGAACAATATGGTCATGCGTGGTGTGC  >  NC_000913/3142443‑3142696

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: