Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,053,037	A→G	V109A (GTC→GCC)	yccM ←	putative 4Fe‑4S membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,053,037	0	A	G	100.0%	165.7 / NA	44	V109A (GTC→GCC)	yccM	putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/15);  total (29/15)

ACGCCGCGCAGAGGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACT  >  NC_000913/1052893‑1053182                                                                                                                                                 |                                                                                                                                                  aCGCCGCGCAGAGGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATa                                                                                                                                              >  1:48621/1‑150 (MQ=255)  cgccgcGCAGAGGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATaa                                                                                                                                             >  2:186391/1‑150 (MQ=255)   gccgcGCAGAGGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATaaa                                                                                                                                            >  2:232253/1‑150 (MQ=255)           gagGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCa                                                                                                                                    >  1:315535/1‑150 (MQ=255)            agGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAAt                                                                                                                                  <  2:543837/151‑1 (MQ=255)                   tCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGc                                                                                                                            >  1:678991/1‑150 (MQ=255)                      agccagcGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATcg                                                                                                                        >  2:193523/1‑151 (MQ=255)                         cagcGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCgcgg                                                                                                                     >  2:289114/1‑151 (MQ=255)                            cGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCgcggc                                                                                                                    >  1:561205/1‑149 (MQ=255)                            cGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCgcggc                                                                                                                    <  2:561188/149‑1 (MQ=255)                                gAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTg                                                                                                              >  2:28125/1‑151 (MQ=255)                                 aaGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTgg                                                                                                             >  1:512387/1‑151 (MQ=255)                                  aGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTTCCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTgga                                                                                                            <  2:442964/151‑1 (MQ=255)                                  aGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTgg                                                                                                             >  1:643871/1‑150 (MQ=255)                                        caTTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGgacggacg                                                                                                      >  1:155948/1‑151 (MQ=255)                                              cGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACt                                                                                                >  2:60349/1‑151 (MQ=255)                                              cGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACt                                                                                                <  2:126702/151‑1 (MQ=255)                                              cGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACt                                                                                                <  2:158465/151‑1 (MQ=255)                                                  cAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGAcc                                                                                            >  1:502440/1‑151 (MQ=255)                                                   aaaCAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCt                                                                                           >  2:449012/1‑151 (MQ=255)                                                     aCAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGt                                                                                         >  2:71374/1‑151 (MQ=255)                                                       aGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGtt                                                                                        >  1:318250/1‑150 (MQ=255)                                                       aGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTa                                                                                       <  1:582041/151‑1 (MQ=255)                                                           tGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACaa                                                                                    <  1:526052/150‑1 (MQ=255)                                                               aCCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCt                                                                                >  1:169207/1‑150 (MQ=255)                                                                 cGAGACCGCCGATTAATTCAGAAAGCGTACCAACCGGTCATAACCATGAGCAGAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGGCTGGACGGACGGGAAGTGATCTGTTACatagc                                                                                 >  2:109979/1‑143 (MQ=255)                                                                  gagaTCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTatat                                                                            <  2:334910/151‑1 (MQ=255)                                                                      tCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTc                                                                         >  1:527315/1‑150 (MQ=255)                                                                          cGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGAtt                                                                    >  1:670426/1‑151 (MQ=255)                                                                                   cAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGcc                                                            <  1:565693/150‑1 (MQ=255)                                                                                        aGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCa                                                       >  1:29455/1‑150 (MQ=255)                                                                                            tACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCgg                                                  >  2:350697/1‑151 (MQ=255)                                                                                                 aCCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTagcc                                             >  1:260901/1‑151 (MQ=255)                                                                                                   cGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTagcc                                             >  1:139770/1‑149 (MQ=255)                                                                                                   cGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTagcc                                             >  1:544049/1‑149 (MQ=255)                                                                                                    gggCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTagccagc                                          >  1:33254/1‑151 (MQ=255)                                                                                                      gCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCt                                        <  1:614480/151‑1 (MQ=255)                                                                                                           aCCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGAt                                   <  1:71376/151‑1 (MQ=255)                                                                                                                 aGCAAAAGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTcc                             <  2:75793/151‑1 (MQ=255)                                                                                                                     aaaGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGt                         <  1:275017/151‑1 (MQ=255)                                                                                                                       aGGCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCt                       <  1:507679/151‑1 (MQ=255)                                                                                                                        ggCCTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTg                      >  1:481101/1‑151 (MQ=255)                                                                                                                          ccTTTTTGAGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGc                    <  1:350032/151‑1 (MQ=255)                                                                                                                                  aGTAATAGACTGATGGCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATg            >  2:415879/1‑151 (MQ=255)                                                                                                                                            tgatggCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACt  >  1:416441/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ACGCCGCGCAGAGGAATATCCAGCCAGCGGGGAAGGACACATTGCCGACCAAACAGTTTGTTACCGAGATCGCCGATTAATTCAGAAAGCGTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACT  >  NC_000913/1052893‑1053182

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