Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,041,654	G→A	intergenic (+29/‑125)	dsbA → / → b3861	protein disulfide isomerase I/putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,041,654	0	G	A	100.0%	203.6 / NA	55	intergenic (+29/‑125)	dsbA/b3861	protein disulfide isomerase I/putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/53);  total (2/53)

ATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCA  >  NC_000913/4041564‑4041733                                                                                           |                                                                                 atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                                <  2:3125471/93‑1 (MQ=255) atatGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATa                                                                                <  1:445804/93‑1 (MQ=255)   atGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATAtt                                                                              <  2:3536702/93‑1 (MQ=255)   atGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATAtt                                                                              <  2:4029076/93‑1 (MQ=255)    tGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTc                                                                             <  2:2191189/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                            <  1:1421443/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                            <  2:4206014/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                            <  1:3473759/93‑1 (MQ=255)     ggATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCa                                                                            <  1:2653442/93‑1 (MQ=255)      gATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAc                                                                           <  1:1081671/93‑1 (MQ=255)       aTGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACg                                                                          <  2:2715366/93‑1 (MQ=255)       aTGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACg                                                                          <  1:2088241/93‑1 (MQ=255)        tGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACgg                                                                         <  2:385443/93‑1 (MQ=255)        tGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACgg                                                                         <  2:2128050/93‑1 (MQ=255)         gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACggg                                                                        <  1:13264/93‑1 (MQ=255)         gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACggg                                                                        <  2:1134986/93‑1 (MQ=255)          tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                       <  1:3580747/93‑1 (MQ=255)           ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑t                                                                     <  2:4499597/93‑1 (MQ=255)            tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tt                                                                    <  1:42162/93‑1 (MQ=255)            tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tt                                                                    <  2:3579131/93‑1 (MQ=255)             ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                                                   <  1:451177/93‑1 (MQ=255)             ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                                                   <  1:527036/93‑1 (MQ=255)              tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                                                  <  2:3519855/93‑1 (MQ=255)              tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                                                  <  2:3519856/93‑1 (MQ=255)               gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                                 <  2:412060/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                                <  2:2820366/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                                <  2:849789/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                                <  1:928726/93‑1 (MQ=255)                ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                                <  1:4229536/93‑1 (MQ=255)                 tCAGCAGTATGCTGATACAGTGAAATCTCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttt                                                               <  2:3793203/93‑2 (MQ=255)                 tCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttt                                                               <  1:498266/93‑2 (MQ=255)                 tCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttt                                                               <  2:3223758/93‑2 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                              <  1:879806/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                              <  2:1476132/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                              <  2:2093103/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                              <  1:1989427/93‑3 (MQ=255)                  cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                              <  1:3592655/93‑3 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                             <  2:1570734/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                             <  1:1895694/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                             <  2:461359/93‑4 (MQ=255)                   agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                             <  1:2439297/93‑4 (MQ=255)                    gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                                            <  2:921912/93‑5 (MQ=255)                    gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                                            <  2:450523/93‑5 (MQ=255)                    gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                                            <  2:3389643/93‑5 (MQ=255)                    gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                                            <  1:4225896/93‑5 (MQ=255)                      agTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                           <  2:1022147/93‑1 (MQ=255)                       gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                          <  1:1553296/93‑1 (MQ=255)                        tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                         <  2:1801506/93‑1 (MQ=255)                         aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTta                                                        <  1:2387308/93‑1 (MQ=255)                           gCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaat                                                      >  2:2374121/1‑93 (MQ=255)                           gCTGATACAGGGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaat                                                      <  1:3588848/93‑1 (MQ=255)                            cTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaata                                                     <  1:2225274/93‑1 (MQ=255)                             tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                    <  1:2140804/93‑1 (MQ=255)                                                                        aaTTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACattat         >  1:4026167/1‑93 (MQ=255)                                                                               tAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCa  <  1:1053309/93‑1 (MQ=255)                                                                                           |                                                                                 ATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCA  >  NC_000913/4041564‑4041733

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