Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	45,545	A→T	9.6%	K28M (AAG→ATG)	fixX →	putative ferredoxin FixX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	45,545	0	A	T	9.6%	125.5 / 3.9	53	K28M (AAG→ATG)	fixX	putative ferredoxin FixX
Reads supporting (aligned to +/- strand):  ref base A (28/18);  new base T (2/3);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CGGCAAGAAAGTGGGCTTCATCAATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGG  >  NC_000913/45403‑45681                                                                                                                                               |                                                                                                                                         cGGCGAGCAGGTGGGCTTCATCACTCTAATCTAGGAGGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGaa                                                                                                                                          <  1:315988/143‑1 (MQ=255)  ggCAAGAAAGTGGGCTTCATCAATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGaat                                                                                                                                         >  1:445488/1‑142 (MQ=255)     aaGAAAGTGGGCTTCATCAATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGcc                                                                                                                                      >  1:444040/1‑143 (MQ=255)        aaaGTGGGCTTCATCAATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGAt                                                                                                                                   <  1:763899/143‑1 (MQ=255)                 ttCATCAATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATaaa                                                                                                                          <  1:713112/143‑1 (MQ=255)                  tcatcaATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGGTaacc                                                                                                                         >  2:462755/1‑141 (MQ=255)                     tcaAGCTAAGCAAGGAGGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTGAATGTGGACGTCAAACTGGGCGGCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAgg                                                                                                                      <  2:366841/143‑1 (MQ=255)                             aTCAATGATGGCATGTAAGGAGTGACCGTTTTATTACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAg                                                                                                              <  2:740971/143‑1 (MQ=255)                               cAAGGAGGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAgct                                                                                                            <  1:313325/143‑1 (MQ=255)                                  ggATGGCATGAAAGGAGTGAACGTTTTATGACTTCTCCCGTAAATGTGGACGTCAAACTGAGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACAATGTTGTGAAGGACGATACTAATGAACAGGCGCTGAAggagca                                                                                                         >  2:451527/1‑138 (MQ=255)                                         aTGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCCATGTGGACGTCAAACTGGGCGTCAATAAATTCTATGTCGATGTTGTGTTTAATTGCAGAAATGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGATGGTGaaa                                                                                                  >  1:573728/1‑143 (MQ=255)                                          tGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAg                                                                                                 >  1:382230/1‑143 (MQ=255)                                          tGAAAGGAGTGACCGTTTTATGACTCCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATACATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAg                                                                                                 <  2:789313/143‑1 (MQ=255)                                             aaGGAGTGACCGTTTTATGACTTCTCCCGTCAAAGTAGAAGCCAAACTGGGCCTCAGTAAATTCAATGTCGATGAAGAGAATCCGCACCTTGTAGTCAATGCTGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGCAAGCGt                                                                                              >  1:107953/1‑143 (MQ=255)                                             aaGGAGTGACAGGTGTATGGCTTCTCCCGTCAATGTGGAAGTCAAACTGGGCGTCAAGAAATGCAATGTCGATGAAGAGCATCAGCACATTGTTGTTAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGt                                                                                              <  1:183424/143‑1 (MQ=255)                                                   tGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGACGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCg                                                                                        >  1:409475/1‑143 (MQ=255)                                                   tGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACGGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACAGTGCGGTGAAGGCCGATGCTGATAAACAGGCGCGGGAGCGGCTGGTGAAAGCGTGCCCCg                                                                                        >  2:317367/1‑143 (MQ=255)                                                     aCCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCa                                                                                      >  1:86467/1‑143 (MQ=255)                                                     aCCGTTTTAGGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCAGCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCa                                                                                      <  2:119350/143‑1 (MQ=255)                                                          tttATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCt                                                                                 <  1:452020/143‑1 (MQ=255)                                                              cgACTTCCCCCGTAAATGTGGACGTCAAACGGGGCGGCAATAAATTGAAGGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTAc                                                                             <  1:196704/142‑1 (MQ=255)                                                                aCTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAACTTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGTCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCACGTCTGTAcca                                                                           >  2:30709/1‑141 (MQ=255)                                                                aCTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCGGGTGTGTACaa                                                                           >  1:516570/1‑143 (MQ=255)                                                                  ttCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGCGGCGGATAAACAGGCGCTGGAGCTGCTGGTGACAGCGTGCCCAGCAGGGCTGTACaaga                                                                         >  2:40287/1‑143 (MQ=255)                                                                     tcCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGATTCTACAAGACGc                                                                      >  2:43748/1‑143 (MQ=255)                                                                        cGTAAAAGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACACGCAGCggg                                                                   >  2:135999/1‑143 (MQ=255)                                                                         gTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGa                                                                  >  1:643101/1‑143 (MQ=255)                                                                         gTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTGGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGCAAGCGTGCCCCGCAGGACTGTGCACGAAGCAGGa                                                                  >  2:372993/1‑143 (MQ=255)                                                                           ccaTGTGGACGTCAACCCGGTCGCCAATGAATTCAATGTCGATGAAGCGCATCCGCACATTGTTGTGAGGGTCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATg                                                                <  2:608772/141‑1 (MQ=255)                                                                            aaTGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAggcgga                                                               >  2:134099/1‑139 (MQ=255)                                                                               gtgGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCCGTACAAGAAGCAGGATGACgg                                                            >  2:638040/1‑143 (MQ=255)                                                                                 cgACGTCAACCTGGGCGTCAATGCACTCACTGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATCAACCGGCGCTGGAGTTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCa                                                          <  1:112468/142‑1 (MQ=255)                                                                                  gACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGATGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAg                                                         >  1:496777/1‑143 (MQ=255)                                                                                  gACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGATGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAg                                                         >  1:115688/1‑143 (MQ=255)                                                                                  gACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGCAAGCGTGCCCCGCAGGTCTGAACAAGAGGCACGATGACGGCAg                                                         >  1:272043/1‑143 (MQ=255)                                                                                    cGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAGGCAGGATGACGGCAgtg                                                       >  1:64214/1‑143 (MQ=255)                                                                                       cAAACTGGTCGTAAATAAAGTAAATGTTGATGAAGAGCATCCTCATATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTgc                                                    <  2:32135/143‑1 (MQ=255)                                                                                        aaaCTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTgcg                                                   >  2:663521/1‑143 (MQ=255)                                                                                        aaaCTGGGCGTCAATAAATTCAATGTCGATAAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGACAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGAGGACGGCGGTGTgcg                                                   >  1:578872/1‑143 (MQ=255)                                                                                                     aTAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGcc                                      >  1:689475/1‑143 (MQ=255)                                                                                                      tAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCg                                     <  1:531081/143‑1 (MQ=255)                                                                                                               tGTCGTTGAAGAGCATCCGCACATTGTTGTGATGGCCGATGCTGATGAACAGGCTCTTGTGCTGCTGGTGAAAGCTTTTCCCGCTGGTCTGTACAAGAACCATGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTgg                            <  1:686238/143‑1 (MQ=255)                                                                                                                gacgctGAAGAGCATCCGCACATTGTTGTGACGTCCTATGCTGATAAACAGGCGCTGGAGCCCCTGGTGAAAGCGTGCCCCTCAGTTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGa                           <  1:246427/138‑1 (MQ=255)                                                                                                                gacGATGAAGAGCATCCGCACATGGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGGGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGa                           <  2:108584/141‑1 (MQ=255)                                                                                                                 tCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTgggg                          >  2:26848/1‑141 (MQ=255)                                                                                                                    aTGTAGAGCATTCGTTTATTGTTTTGAAGGTCGATGCTGAGAAACAGGGGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGAAGGATGACGGAAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGc                       <  1:369349/143‑1 (MQ=255)                                                                                                                    aTGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCCTCGATTACGCCGGATGTGTGGAGTGc                       >  1:751235/1‑143 (MQ=255)                                                                                                                    aTGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGGCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGc                       >  1:659026/1‑143 (MQ=255)                                                                                                                     tGAAGAGCATCGGCAGATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCg                      <  2:247713/143‑1 (MQ=255)                                                                                                                         gagCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCAc                  >  2:408262/1‑143 (MQ=255)                                                                                                                              tCCGCACATTCTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGGGCTGCTGGTGAAGGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTc             <  1:315889/143‑1 (MQ=255)                                                                                                                               ccGCACATTGTTGTGATGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCg            <  2:496777/143‑1 (MQ=255)                                                                                                                               ccGCACATTGTTGTGATGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCg            <  2:115688/143‑1 (MQ=255)                                                                                                                                cGCACATTGCTGTAAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGACAGCGTGCCCCGCAGGCCCGTACAAGAAGCAGGATGACGGCGGTGTGCGCCTCGATTGCGACGGATGCGTGGAGTGGGGCACCGGTCGc           >  2:434971/1‑143 (MQ=255)                                                                                                                                         ttgtGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTctgggg  <  1:638040/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                         CGGCAAGAAAGTGGGCTTCATCAATCTAATCAAGGATGGCATGAAAGGAGTGACCGTTTTATGACTTCTCCCGTCAATGTGGACGTCAAACTGGGCGTCAATAAATTCAATGTCGATGAAGAGCATCCGCACATTGTTGTGAAGGCCGATGCTGATAAACAGGCGCTGGAGCTGCTGGTGAAAGCGTGCCCCGCAGGTCTGTACAAGAAGCAGGATGACGGCAGTGTGCGCTTCGATTACGCCGGATGTCTGGAGTGCGGCACCTGTCGCATTCTGGGG  >  NC_000913/45403‑45681

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑