Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,041,671	(T)6→7	intergenic (+46/‑108)	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,665	1	.	T	100.0%	106.8 / NA	29	intergenic (+40/‑114)	dsbA/b3861	protein disulfide isomerase I/putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (3/26);  total (3/26)

TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCA  >  NC_000913/4041573‑4041730                                                                                              |                                                                  tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                    <  2:3909230/93‑1 (MQ=255)  ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGC‑t                                                                  <  2:4723719/93‑1 (MQ=255)  ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑t                                                                  <  2:6651758/93‑1 (MQ=255)   tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tt                                                                 <  1:4764421/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                                                <  2:1026264/93‑1 (MQ=255)    ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttt                                                                <  2:3988259/93‑1 (MQ=255)     tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                                               <  1:2074507/93‑1 (MQ=255)     tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttt                                                               <  1:7430343/93‑1 (MQ=255)      gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  2:2731165/93‑1 (MQ=255)      gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  2:4090169/93‑1 (MQ=255)      gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  2:6396754/93‑1 (MQ=255)      gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  1:3322499/93‑1 (MQ=255)      gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttt                                                              <  1:2547441/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  1:384467/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  2:3688686/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  1:4423446/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  2:7087095/93‑1 (MQ=255)       ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttt                                                             <  1:2542141/93‑1 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:599164/93‑3 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:5450894/93‑3 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:5107459/93‑3 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:7220882/93‑3 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:4690199/93‑3 (MQ=255)         cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑ttttttta                                                           <  2:7206717/93‑3 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAATAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  1:3395823/93‑4 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:7500769/93‑4 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:4846640/93‑4 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:6272823/93‑4 (MQ=255)          agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttat                                                          <  2:6834227/93‑4 (MQ=255)           gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC‑tttttttatt                                                         <  2:1286216/93‑5 (MQ=255)             agTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                        <  1:3283474/93‑1 (MQ=255)             agTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                        <  2:3678096/93‑1 (MQ=255)              gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:1086984/93‑1 (MQ=255)              gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:3269379/93‑1 (MQ=255)              gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:608877/93‑1 (MQ=255)              gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                       <  1:972897/93‑1 (MQ=255)               tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                      <  2:1307840/93‑1 (MQ=255)               tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                      <  1:4558408/93‑1 (MQ=255)               tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                      <  2:1632553/93‑1 (MQ=255)                aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTta                                                     <  2:4862367/93‑1 (MQ=255)                 tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  1:6527821/93‑1 (MQ=255)                 tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  1:4461691/93‑1 (MQ=255)                 tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  2:1820286/93‑1 (MQ=255)                 tGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaa                                                    <  2:2574690/93‑1 (MQ=255)                  gCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaat                                                   <  2:1797766/93‑1 (MQ=255)                   cTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaata                                                  <  2:2191865/93‑1 (MQ=255)                    tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  2:3834225/93‑1 (MQ=255)                    tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  2:3065511/93‑1 (MQ=255)                    tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  1:4241208/93‑1 (MQ=255)                    tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaataa                                                 <  1:3538297/93‑1 (MQ=255)                     gATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaa                                                <  2:5122377/93‑1 (MQ=255)                         cAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatata                                            <  1:3879526/93‑1 (MQ=255)                            tGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatataaat                                         >  1:2332700/1‑93 (MQ=255)                             gAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAtaaatataaata                                        <  2:3374021/93‑1 (MQ=255)                                 tatCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACAtt                                    <  1:3988715/93‑1 (MQ=255)                                  atCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTc                                   <  1:7029741/93‑1 (MQ=255)                                             aaaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCAt                        >  1:560180/1‑93 (MQ=255)                                                             tAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACatt        <  1:6773556/93‑1 (MQ=255)                                                                   aCATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCa  >  1:6590349/1‑93 (MQ=255)                                                                                              |                                                                  TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGC‑TTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCA  >  NC_000913/4041573‑4041730

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