Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,299,611 T→C 7.0% H217R (CAC→CGC)  napG ← ferredoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,299,6110TC 7.0% 171.0 / 3.8 57H217R (CAC→CGC) napGferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base C (2/2);  total (27/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.26e-01

CATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGT  >  NC_000913/2299465‑2299756
                                                                                                                                                  |                                                                                                                                                 
cATCCCCAGCACGAAGAACTGAAAAAGGCGACGTAAGACCAGCCAACGGTGACTGCGCCTCCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTg                                                                                                                                                  <  6:297835/148‑1 (MQ=255)
 aTCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGa                                                                                                                                                 <  1:315162/148‑1 (MQ=255)
   ccccAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGAcc                                                                                                                                               >  8:16418/1‑148 (MQ=255)
      cAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTaa                                                                                                                                            <  4:82585/148‑1 (MQ=255)
             aagaaCTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCt                                                                                                                                     >  8:274247/1‑148 (MQ=255)
             aagaaCTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACCTTAACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCt                                                                                                                                     >  2:73868/1‑148 (MQ=255)
              agaaCTGACAAAGGCGACGTAACACCATCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGAGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCtt                                                                                                                                    <  2:189729/148‑1 (MQ=255)
                     aCAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCa                                                                                                                              <  5:109965/147‑1 (MQ=255)
                         aGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGAc                                                                                                                         <  2:19464/148‑1 (MQ=255)
                          ggCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGCGACCTAACTCCCCTTTCGCCAGTGa                                                                                                                          >  4:29247/1‑146 (MQ=255)
                          ggCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCATCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGAGGTAATGGCGACCAAACTCCCCATTCGCCAGTGa                                                                                                                          >  6:279348/1‑146 (MQ=255)
                           gCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAg                                                                                                                       >  5:16879/1‑148 (MQ=255)
                               cGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCgg                                                                                                                    >  4:162646/1‑147 (MQ=255)
                                  aaCACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGcc                                                                                                                               <  4:77280/133‑1 (MQ=255)
                                  aaCACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGcc                                                                                                                               >  3:77280/1‑133 (MQ=255)
                                   acacCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAac                                                                                                                <  4:162396/147‑1 (MQ=255)
                                     acCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACAcc                                                                                                             <  6:180753/148‑1 (MQ=255)
                                        cgCCAACGGGGTCTTCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGGTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAAATCCCCTTTCGCCAGTGACAGCGGTAACACCTTg                                                                                                          <  6:69877/147‑1 (MQ=255)
                                           cAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCACGCCAGCCGAAGCGGTAATGGTGCACTAACTCCCATTTCGCCAGTGACAGCGGTAACAACTTGAtt                                                                                                       >  4:157985/1‑148 (MQ=255)
                                              cGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGcc                                                                                                    >  6:136729/1‑148 (MQ=255)
                                               ggTGCCTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCg                                                                                                   >  2:121874/1‑148 (MQ=255)
                                                   aCTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGgtt                                                                                                <  3:162646/147‑1 (MQ=255)
                                                      gcgcCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCCGCCGCAGCGGTAATGGTGCCCTAACTCCCCGTTCGCCAGTGACAGCGGTAACACCTTGATTGCCgggttac                                                                                             >  5:233239/1‑142 (MQ=255)
                                                          caccaGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGc                                                                                        <  6:16879/148‑1 (MQ=255)
                                                           accaGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCa                                                                                       >  1:316267/1‑148 (MQ=255)
                                                              aGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGgt                                                                                                 >  7:184453/1‑135 (MQ=255)
                                                              aGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGgt                                                                                                 <  8:184453/135‑1 (MQ=255)
                                                                ccTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTGCAGCCAGCCGAAGCGGTAATGGCGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCAt                                                                                  <  5:279348/148‑1 (MQ=255)
                                                                ccTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGCGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCAt                                                                                  <  3:29247/148‑1 (MQ=255)
                                                                       tctcCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACtttt                                                                            <  2:316267/147‑1 (MQ=255)
                                                                          ccAGCGCCTCGCGCCCGGCGTCACGTTTATGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCg                                                                         >  2:55611/1‑147 (MQ=255)
                                                                           cAGCGCCTCGCGCCCGTCGTCACGTTTACGATTTGCCATTGTTACCCTCCATACAGCCGAAGCGGTAATGGTTACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTc                                                                          <  1:73868/145‑1 (MQ=255)
                                                                           cAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCa                                                                       <  1:121874/148‑1 (MQ=255)
                                                                               gcCTCGCGCCCGGCGTCACGTTTATGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTa                                                                   <  1:55611/148‑1 (MQ=255)
                                                                                 cTCGCGCCCGGCGTCTCGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTAcc                                                                 <  2:9599/148‑1 (MQ=255)
                                                                                       cccGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCACTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAAcc                                                           >  1:321946/1‑148 (MQ=255)
                                                                                          ggCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGCCAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCgg                                                         >  4:267712/1‑147 (MQ=255)
                                                                                               cACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACg                                                                                     <  3:157512/114‑1 (MQ=255)
                                                                                               cACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACg                                                                                     >  4:157512/1‑114 (MQ=255)
                                                                                                            tGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGaa                                       <  7:274247/147‑1 (MQ=255)
                                                                                                            tGCCAGTGTTCCCCTCCAGCCAGCCGACGCGGTAATGGTGACCTACCTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATGGCCGGGTGTGCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGGGa                                        <  4:49685/146‑1 (MQ=255)
                                                                                                              ccATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCg                                    >  2:91009/1‑148 (MQ=255)
                                                                                                               cATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGt                                   >  1:187309/1‑148 (MQ=255)
                                                                                                               cATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGt                                   >  7:205387/1‑148 (MQ=255)
                                                                                                               cATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCAAACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGt                                   >  4:258984/1‑148 (MQ=255)
                                                                                                                 ttGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCGTACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCg                                 >  6:286117/1‑148 (MQ=255)
                                                                                                                 ttGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCg                                 >  6:59859/1‑148 (MQ=255)
                                                                                                                 ttGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCg                                 <  7:16418/148‑1 (MQ=255)
                                                                                                                  tGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGa                                        <  6:232288/140‑1 (MQ=255)
                                                                                                                  tGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGa                                        >  5:232288/1‑140 (MQ=255)
                                                                                                                           ccagccagCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCg                       >  4:54651/1‑148 (MQ=255)
                                                                                                                              gccagccGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGc                    <  3:54651/148‑1 (MQ=255)
                                                                                                                                     gAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGGGCTTa             <  5:59859/148‑1 (MQ=255)
                                                                                                                                       aGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTAcc           <  3:278673/148‑1 (MQ=255)
                                                                                                                                              aTGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTAcg     <  5:136729/147‑1 (MQ=255)
                                                                                                                                               tGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTAcgcg   <  5:42535/148‑1 (MQ=255)
                                                                                                                                                ggTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACTGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCgt  >  4:10761/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAGCCGAAGCGGTAATGGTGACCTAACTCCCCTTTCGCCAGTGACAGCGGTAACACCTTGATTGCCGGTTGTTCCAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGT  >  NC_000913/2299465‑2299756

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