Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,486	2 bp→GC	coding (316‑317/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,486	0	C	G	100.0%	165.3 / NA	64	R106A (CGG→GCG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (24/40);  total (24/40)
*	NC_000913	4,001,487	0	G	C	100.0%	163.0 / NA	65	R106A (CGG→GCG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (24/41);  total (24/41)

GTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGA  >  NC_000913/4001422‑4001556                                                                 ||                                                                      ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                             >  1:1786733/2‑76 (MQ=255)  tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                            <  1:1754562/76‑1 (MQ=255)  gcGCCGATGGAAATCATCCCCAGCACAAGGTTCACCAGCGGGTTAATAAAGTAGCCAAGGCTCGCTTCCAGCatat                                                            <  2:423225/75‑1 (MQ=255)   cGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg                                                           <  2:2735737/76‑1 (MQ=255)    gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                          >  1:4595171/1‑76 (MQ=255)    gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                          >  1:350960/1‑76 (MQ=255)    gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                          <  1:286615/76‑1 (MQ=255)      cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTg                                                        >  1:915403/1‑76 (MQ=255)       gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa                                                       >  2:4649541/1‑76 (MQ=255)          gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTg                                                    <  2:3987662/76‑1 (MQ=255)            aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt                                                  <  2:1591555/76‑1 (MQ=255)             aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                 >  1:3705781/1‑76 (MQ=255)             aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                 >  1:446076/1‑76 (MQ=255)              atcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCa                                                <  2:3791457/76‑1 (MQ=255)              atcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCa                                                <  1:3615215/76‑1 (MQ=255)                cacccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                              <  2:640739/73‑1 (MQ=255)                 atcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                 <  1:682445/72‑1 (MQ=255)                 atcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGCGATTGTTc                                                 >  2:682445/1‑72 (MQ=255)                    cccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                          <  1:1793793/76‑1 (MQ=255)                    cccAGCACAATGTGCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                          <  1:938492/76‑1 (MQ=255)                           caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                   <  2:823897/76‑1 (MQ=255)                           caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                   >  1:3822935/1‑76 (MQ=255)                            aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAAc                                  <  1:2003740/76‑1 (MQ=255)                              tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                >  2:391931/1‑76 (MQ=255)                              tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                <  1:2174347/76‑1 (MQ=255)                               gTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                               <  1:2097950/76‑1 (MQ=255)                               gTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                               <  2:4208611/76‑1 (MQ=255)                                  cACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGc                            <  1:4784532/76‑1 (MQ=255)                                  cACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGc                            >  2:1124745/1‑76 (MQ=255)                                   aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCCGATAAACAGAAGcc                           >  2:54285/1‑76 (MQ=255)                                   aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                           >  1:1282767/1‑76 (MQ=255)                                    ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCa                          >  1:133536/1‑76 (MQ=255)                                     cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                         <  2:3237708/76‑1 (MQ=255)                                      ggggggTTAATAAAGTAACAAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGt                        <  2:546567/73‑1 (MQ=255)                                      aGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGt                        <  2:1680390/76‑1 (MQ=255)                                       gCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                       <  2:1986873/76‑1 (MQ=255)                                        cGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTg                      >  1:3538868/1‑76 (MQ=255)                                         gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                     >  2:4553585/1‑76 (MQ=255)                                          gagtaaaaaAGGAACCAGGGGCCGCTTCCAGCATATGGTGATTTTTCACCGCCCAGATAAACAGTAGCCAGTTGcc                    <  2:184900/73‑1 (MQ=255)                                           gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                   >  2:4037634/1‑76 (MQ=255)                                            ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                  >  1:1003210/1‑76 (MQ=255)                                            ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                  <  2:3235636/76‑1 (MQ=255)                                            ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGCTGccac                  >  2:4383647/1‑76 (MQ=255)                                             taataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                 <  2:3796655/76‑1 (MQ=255)                                             taataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                 <  1:3251297/76‑1 (MQ=255)                                             taataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                 <  2:4003237/76‑1 (MQ=255)                                               ataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCaa               >  2:4384099/1‑76 (MQ=255)                                                aaaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCATATAAACAGTAGCCAGTTGCCACCAAt              <  1:4822582/75‑1 (MQ=255)                                                 aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc             >  2:1614537/1‑76 (MQ=255)                                                  aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACATTAGCCAGTTGCCACCAATCa            <  2:4641158/76‑1 (MQ=255)                                                  aaGTAACCAAGGCTCGCTTCCAGCATATGGTGAGTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa            <  1:4751650/76‑1 (MQ=255)                                                  aaGGAACCAAGGCTGGCATCCAGGTTATGGTGGTTGTTCACCGCCCATATAAACAGTAGCCAGTTGCCACCAATCa            <  2:2001806/76‑1 (MQ=255)                                                     tAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa         >  1:758136/1‑76 (MQ=255)                                                      aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc        <  1:4291606/76‑1 (MQ=255)                                                      aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc        <  1:1510769/76‑1 (MQ=255)                                                       aCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg       <  1:1588052/76‑1 (MQ=255)                                                       aCCAAGGCTCGCTTCCAGCATAGGGTGATTGTTCGCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg       <  2:4543339/76‑1 (MQ=255)                                                        ccAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGCTGCCACCAATCAGCACgg      >  2:4524816/1‑76 (MQ=255)                                                         cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc     <  1:1557092/76‑1 (MQ=255)                                                          aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa    <  1:757413/76‑1 (MQ=255)                                                          aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa    <  1:712255/76‑1 (MQ=255)                                                          aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa    >  1:4105357/1‑76 (MQ=255)                                                          aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa    <  1:1775999/76‑1 (MQ=255)                                                            ggCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga  <  2:1235808/76‑1 (MQ=255)                                                            ggCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga  <  1:127820/76‑1 (MQ=255)                                                                 ||                                                                      GTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGA  >  NC_000913/4001422‑4001556

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑