Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,761,556	A→C	L79* (TTA→TGA)	sufD ←	Fe‑S cluster scaffold complex subunit SufD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,761,556	0	A	C	100.0%	304.5 / NA	75	L79* (TTA→TGA)	sufD	Fe‑S cluster scaffold complex subunit SufD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (37/38);  total (37/38)

GCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCT  >  NC_000913/1761509‑1761602                                                |                                               gCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCAc                                              <  1:443029/50‑1 (MQ=255)  cACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACg                                             <  1:5813933/50‑1 (MQ=255)   aCGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGc                                            >  1:2949937/1‑50 (MQ=255)     gTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCt                                           >  1:4819728/1‑49 (MQ=255)     gTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTa                                          <  1:4965450/50‑1 (MQ=255)     gTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTa                                          <  1:179087/50‑1 (MQ=255)     gTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTa                                          <  1:2129509/50‑1 (MQ=255)       aaCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                        >  1:1552009/1‑50 (MQ=255)       aaCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                        >  1:4242488/1‑50 (MQ=255)        aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAgg                                       <  1:2862759/50‑1 (MQ=255)        aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAgg                                       >  1:4059319/1‑50 (MQ=255)        aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                        <  1:518913/49‑1 (MQ=255)        aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                        <  1:4231962/49‑1 (MQ=255)        aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                        <  1:1822817/49‑1 (MQ=255)         cGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGc                                      <  1:1584186/50‑1 (MQ=255)          gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCa                                     >  1:4531869/1‑50 (MQ=255)          gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCa                                     >  1:538875/1‑50 (MQ=255)           cccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCAt                                    <  1:649636/50‑1 (MQ=255)           cccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCAt                                    >  1:3444508/1‑50 (MQ=255)            ccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCAt                                    >  1:3523185/1‑49 (MQ=255)            ccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCAt                                    >  1:5367555/1‑49 (MQ=255)            ccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCAt                                    >  1:3028065/1‑49 (MQ=255)             cATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCa                                  <  1:693239/50‑1 (MQ=255)             cATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCa                                  >  1:5690593/1‑50 (MQ=255)              aTCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCAc                                 >  1:1618963/1‑50 (MQ=255)              aTCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCAc                                 >  1:6475394/1‑50 (MQ=255)               tCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACg                                <  1:2510794/50‑1 (MQ=255)               tCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACg                                <  1:3816914/50‑1 (MQ=255)                 gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                              >  1:2469647/1‑50 (MQ=255)                 gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                              <  1:654196/50‑1 (MQ=255)                 gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                              <  1:3762769/50‑1 (MQ=255)                  aCAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctg                             >  1:6296724/1‑50 (MQ=255)                   cAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgc                            >  1:528600/1‑50 (MQ=255)                    aaaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgct                           <  1:390086/50‑1 (MQ=255)                    aaaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgct                           >  1:3725279/1‑50 (MQ=255)                    aaaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgct                           <  1:4523021/50‑1 (MQ=255)                     aaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGctgctg                          >  1:2047445/1‑50 (MQ=255)                     aaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGctgctg                          <  1:436903/50‑1 (MQ=255)                      acacCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCtgt                         >  1:5005860/1‑50 (MQ=255)                       cacCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCtgtg                        >  1:676582/1‑50 (MQ=255)                        acCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTgg                       <  1:617353/50‑1 (MQ=255)                         ccAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTggg                      <  1:830881/50‑1 (MQ=255)                          cAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTgggg                     <  1:5873002/50‑1 (MQ=255)                          cAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTggg                      >  1:3943232/1‑49 (MQ=255)                           aGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGa                    >  1:2217905/1‑50 (MQ=255)                            gCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGat                   <  1:1615888/50‑1 (MQ=255)                             ccGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGata                  >  1:6232584/1‑50 (MQ=255)                              cGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGatat                 >  1:4931636/1‑50 (MQ=255)                               gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                >  1:4664303/1‑50 (MQ=255)                               gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                <  1:4418923/50‑1 (MQ=255)                               gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                >  1:2726914/1‑50 (MQ=255)                                cACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtct               >  1:5747682/1‑50 (MQ=255)                                 aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc              >  1:464109/1‑50 (MQ=255)                                 aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc              >  1:4595110/1‑50 (MQ=255)                                 aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc              <  1:59514/50‑1 (MQ=255)                                  cGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctct             <  1:2552112/50‑1 (MQ=255)                                  cGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctct             >  1:1797455/1‑50 (MQ=255)                                   ggAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctctc            >  1:1012428/1‑50 (MQ=255)                                    gAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTcc           <  1:564491/50‑1 (MQ=255)                                     aGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTccc          >  1:5016838/1‑50 (MQ=255)                                     aGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTccc          >  1:802861/1‑50 (MQ=255)                                     aGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTccc          <  1:2911975/50‑1 (MQ=255)                                      gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg         <  1:5666975/50‑1 (MQ=255)                                      gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg         <  1:3337306/50‑1 (MQ=255)                                         tAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCaa      >  1:1101026/1‑50 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCaa      <  1:4150252/48‑1 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAAt     <  1:5901298/49‑1 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAAt     <  1:4179947/49‑1 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAAt     <  1:161904/49‑1 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAAt     <  1:920211/49‑1 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATg    <  1:4048201/50‑1 (MQ=255)                                           aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATg    >  1:3911617/1‑50 (MQ=255)                                             gTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCt  <  1:6476302/50‑1 (MQ=255)                                             gTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCt  >  1:3214074/1‑50 (MQ=255)                                             gTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCt  <  1:2932789/50‑1 (MQ=255)                                                |                                               GCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCT  >  NC_000913/1761509‑1761602

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