Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,368,320	A→G	K81R (AAA→AGA)	pspA →	phage shock protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,368,320	0	A	G	100.0%	170.3 / NA	48	K81R (AAA→AGA)	pspA	phage shock protein A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/29);  total (19/29)

ATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAAAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAA  >  NC_000913/1368184‑1368456                                                                                                                                         |                                                                                                                                         aTGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTgagaga                                                                                                                                       <  1:524169/140‑1 (MQ=255)  tGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGagagag                                                                                                                                      >  1:1166168/1‑140 (MQ=255)   ggAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTgagagaga                                                                                                                                     >  2:970661/1‑140 (MQ=255)   ggAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTgagagaga                                                                                                                                     >  2:541097/1‑140 (MQ=255)   ggAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTgagagaga                                                                                                                                     >  1:251289/1‑140 (MQ=255)      aGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGagagagag                                                                                                                                  <  2:835719/140‑1 (MQ=255)      aGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGagagagag                                                                                                                                  <  1:197236/140‑1 (MQ=255)      aGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGagagagag                                                                                                                                  <  1:541097/140‑1 (MQ=255)           cacTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGAt                                                                                                                             <  2:843949/140‑1 (MQ=255)             cTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCt                                                                                                                           <  2:1087945/140‑1 (MQ=255)               ggTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTgg                                                                                                                         >  1:1075941/1‑140 (MQ=255)               ggTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTgg                                                                                                                         >  1:226087/1‑140 (MQ=255)               ggTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTgg                                                                                                                         >  1:184074/1‑140 (MQ=255)                gTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGc                                                                                                                        <  2:605554/140‑1 (MQ=255)                   gAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACg                                                                                                                     >  1:443264/1‑140 (MQ=255)                          gTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCg                                                                                                              >  1:1083924/1‑140 (MQ=255)                          gTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCg                                                                                                              >  2:780513/1‑140 (MQ=255)                          gTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCg                                                                                                              >  2:192167/1‑140 (MQ=255)                                  tCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGa                                                                                                      <  2:443264/140‑1 (MQ=255)                                  tCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGa                                                                                                      <  2:251289/140‑1 (MQ=255)                                      cgcgcgcgTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGaaaaa                                                                                                  >  1:819156/1‑140 (MQ=255)                                                ggCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGa                                                                                        <  2:226087/140‑1 (MQ=255)                                                ggCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGa                                                                                        <  2:577949/140‑1 (MQ=255)                                                ggCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGACAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGa                                                                                        >  2:265791/1‑140 (MQ=255)                                                     aaaagaaaCAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGAt                                                                                   <  2:1075941/140‑1 (MQ=255)                                                          aaaCAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGGCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGAt                                                                              <  2:717295/140‑1 (MQ=255)                                                          aaaCAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGAt                                                                              <  1:1089823/140‑1 (MQ=255)                                                           aaCAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGAtt                                                                             >  2:444842/1‑140 (MQ=255)                                                            aCAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTa                                                                            <  1:879710/140‑1 (MQ=255)                                                               gCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGt                                                                         <  2:184074/140‑1 (MQ=255)                                                                   aCTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCt                                                                     <  1:780513/140‑1 (MQ=255)                                                                          gTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACAt                                                              <  2:1166168/140‑1 (MQ=255)                                                                                aaCAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTg                                                        <  2:819156/140‑1 (MQ=255)                                                                                 aCAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGa                                                       >  1:492293/1‑140 (MQ=255)                                                                                 aCAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGa                                                       <  2:492293/140‑1 (MQ=255)                                                                                   aaGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACg                                                     <  1:1116038/140‑1 (MQ=255)                                                                                     gCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCt                                                   >  1:883904/1‑140 (MQ=255)                                                                                     gCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCt                                                   >  1:582120/1‑140 (MQ=255)                                                                                     gCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCt                                                   >  2:273354/1‑140 (MQ=255)                                                                                       gTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCtgg                                                 <  1:858325/140‑1 (MQ=255)                                                                                        tCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGa                                                       >  1:658958/1‑133 (MQ=255)                                                                                        tCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGa                                                       <  2:658958/133‑1 (MQ=255)                                                                                                tGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATa                                        <  1:970661/140‑1 (MQ=255)                                                                                                      tGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTgg                                  <  2:831352/140‑1 (MQ=255)                                                                                                              aGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCa                            <  2:52333/138‑1 (MQ=255)                                                                                                              aGGAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCa                            >  1:52333/1‑138 (MQ=255)                                                                                                                gAAAAAGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATgaa                        >  1:1065841/1‑140 (MQ=255)                                                                                                                 aaaaaGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGAt                                         <  1:200081/122‑1 (MQ=255)                                                                                                                 aaaaaGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGAt                                         >  2:200081/1‑122 (MQ=255)                                                                                                                 aaaaaGCCGAACTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATgaag                       <  1:273354/140‑1 (MQ=255)                                                                                                                          aaCTGGCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAGCAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGAtt              >  2:281857/1‑140 (MQ=255)                                                                                                                              ggCGCTGCTGAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCg          <  1:444842/140‑1 (MQ=255)                                                                                                                                      tgAGAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGaa  <  2:582120/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         ATGGAAGATACACTGGTTGAAGTACGTTCTACTTCGGCGCGCGCGTTGGCAGAAAAGAAACAGCTGACTCGCCGTATTGAACAAGCGTCGGCGCGTGAGGTTGAATGGCAGGAAAAAGCCGAACTGGCGCTGCTGAAAGAGAGAGAGGATCTGGCACGTGCAGCGTTAATTGAAAAACAGAAACTGACCGATCTGATTAAGTCCCTGGAACATGAAGTGACGCTGGTGGACGATACGCTGGCACGCATGAAGAAAGAGATTGGCGAGCTGGAA  >  NC_000913/1368184‑1368456

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