Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,359,379	T→C	V61A (GTA→GCA)	gltD →	glutamate synthase, 4Fe‑4S protein, small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,359,379	0	T	C	100.0%	121.8 / NA	35	V61A (GTA→GCA)	gltD	glutamate synthase, 4Fe‑4S protein, small subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/16);  total (19/16)

GCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGTACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACGAGTCTGCCCGCAAGACCGTCTGTGCGAAGGTTCCTGCACT  >  NC_000913/3359233‑3359524                                                                                                                                                   |                                                                                                                                                  gcgcGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGcaca                                                                                                                                                <  1:95871/150‑1 (MQ=255)        gATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTAc                                                                                                                                         <  2:350219/150‑1 (MQ=255)         aTCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACAt                                                                                                                                       <  2:134662/151‑1 (MQ=255)          tCCGCCAAAGAACCCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACa                                                                                                                                        <  2:617586/149‑1 (MQ=255)             gCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATccc                                                                                                                                    <  2:153389/150‑1 (MQ=255)             gCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCg                                                                                                                                   >  1:619565/1‑151 (MQ=255)                         gCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGa                                                                                                                         >  2:499890/1‑149 (MQ=255)                           tGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCt                                                                                                                     >  2:485762/1‑151 (MQ=255)                            gAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTc                                                                                                                    >  2:720235/1‑151 (MQ=255)                             aaGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTc                                                                                                                    <  1:323882/150‑1 (MQ=255)                                aTCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGcc                                                                                                                 >  2:385249/1‑150 (MQ=255)                                  ccGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAAc                                                                                                              <  2:730121/151‑1 (MQ=255)                                                 tgttgAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTAtttt                                                                                                <  1:386520/150‑1 (MQ=255)                                                       aaTTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAgcg                                                                                         >  2:188184/1‑151 (MQ=255)                                                                 ccGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCg                                                                                                      >  2:149865/1‑128 (MQ=255)                                                                 ccGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTa                                                                                                    <  1:149832/130‑1 (MQ=255)                                                                   gTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCg                                                                             >  2:181687/1‑151 (MQ=255)                                                                   gTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCg                                                                             <  2:197191/151‑1 (MQ=255)                                                                   gTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCg                                                                             >  2:458642/1‑151 (MQ=255)                                                                       tCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCAcc                                                                         >  1:184349/1‑151 (MQ=255)                                                                        ccGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCa                                                                        <  2:6173/151‑1 (MQ=255)                                                                         cGAATGCCAGGCCAAAGCGCATGCTGACCGCTGCCTTTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCGCCAg                                                                       <  2:586614/151‑1 (MQ=255)                                                                          gAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGCAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGCACTGCCGcccc                                                                         >  1:100628/1‑148 (MQ=255)                                                                          gAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGa                                                                      <  1:207895/151‑1 (MQ=255)                                                                               ccaggccaAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCaa                                                                  >  1:745816/1‑150 (MQ=255)                                                                                      aaaGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGc                                                          >  1:233705/1‑151 (MQ=255)                                                                                      aaaGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGc                                                          >  2:506641/1‑151 (MQ=255)                                                                                         gcgcAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCgg                                                       <  1:485710/151‑1 (MQ=255)                                                                                               gCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGttt                                                 >  1:467474/1‑151 (MQ=255)                                                                                                cTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTg                                                >  1:209294/1‑151 (MQ=255)                                                                                                     cGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGAc                                           <  2:72924/151‑1 (MQ=255)                                                                                                      gCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACg                                          <  1:777099/151‑1 (MQ=255)                                                                                                       cTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACg                                          >  1:674731/1‑150 (MQ=255)                                                                                                            tGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACGAGTCTg                                    >  1:90629/1‑151 (MQ=255)                                                                                                                                       ggAAATGCCCGGCACACAACTACATCCCGAACTGGCTGAAGCGCGCCAACGAGTGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACGAGTCTGCCCGCAAGACCGTCTGTGCGAAGGtt          <  1:366302/150‑1 (MQ=255)                                                                                                                                              cccGGCACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACGAGTCTGCCCGCAAGACCGTCTGTGCGAAGGTTCCTGCAct  >  1:75357/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAGGCCAAAGCGCAGGCTGACCGCTGCCTGTCGTGCGGCAACCCATACTGCGAGTGGAAATGCCCGGTACACAACTACATCCCGAACTGGCTGAAGCTCGCCAACGAGGGGCGTATTTTTGAAGCGGCGGAACTGTCGCACCAGACCAACACCCTGCCGGAAGTTTGCGGACGAGTCTGCCCGCAAGACCGTCTGTGCGAAGGTTCCTGCACT  >  NC_000913/3359233‑3359524

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