Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,368,615	A→G	E179E (GAA→GAG)	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,615	0	A	G	100.0%	191.0 / NA	55	E179E (GAA→GAG)	pspA	phage shock protein A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/26);  total (29/26)

CGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATA  >  NC_000913/1368481‑1368749                                                                                                                                       |                                                                                                                                       cGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATggaggcgga                                                                                                                                   >  1:786456/1‑140 (MQ=255)             ttGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCca                                                                                                                       <  1:1041755/140‑1 (MQ=255)              tGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCcac                                                                                                                      >  2:155760/1‑140 (MQ=255)                aTGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAg                                                                                                                    >  1:578755/1‑140 (MQ=255)                    tACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTc                                                                                                                >  2:1129052/1‑140 (MQ=255)                    tACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTc                                                                                                                >  1:697690/1‑140 (MQ=255)                       gTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGt                                                                                                             >  1:870739/1‑140 (MQ=255)                          atcaGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGCATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTaaa                                                                                                          >  2:978426/1‑140 (MQ=255)                           tcaGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAAc                                                                                                         >  1:383665/1‑140 (MQ=255)                                   cAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAGACAAAAATCg                                                                                                 >  1:788271/1‑140 (MQ=255)                                    aaaCTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGc                                                                                                <  1:749096/140‑1 (MQ=255)                                           tcgCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGa                                                                                         >  1:190102/1‑140 (MQ=255)                                           tcgCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGa                                                                                         <  2:291821/140‑1 (MQ=255)                                                cgATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGt                                                                                    <  1:215905/140‑1 (MQ=255)                                                            tcAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGa                                                                        <  1:155760/140‑1 (MQ=255)                                                                  ggACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCg                                                                  <  2:190102/140‑1 (MQ=255)                                                                  ggACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCg                                                                  <  2:383665/140‑1 (MQ=255)                                                                   gACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCga                                                                 >  2:351619/1‑140 (MQ=255)                                                                        tGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGatgatg                                                            >  1:918093/1‑140 (MQ=255)                                                                        tGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGatgatg                                                            >  1:892983/1‑140 (MQ=255)                                                                         ggCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGc                                                           <  2:786456/140‑1 (MQ=255)                                                                         ggCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGCATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGc                                                           <  1:978426/140‑1 (MQ=255)                                                                          gCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCa                                                          <  1:1129052/140‑1 (MQ=255)                                                                          gCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCa                                                          >  1:89880/1‑140 (MQ=255)                                                                            aaaCTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAAt                                                        <  1:725534/140‑1 (MQ=255)                                                                             aaCTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATc                                                       >  1:1071778/1‑140 (MQ=255)                                                                                tGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTTGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCatc                                                    >  2:1063842/1‑138 (MQ=255)                                                                                        gCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACt                                            >  2:1117408/1‑140 (MQ=255)                                                                                          aaTGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTgg                                          <  2:892983/140‑1 (MQ=255)                                                                                          aaTGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTgg                                          >  1:570402/1‑140 (MQ=255)                                                                                            tGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGca                                        >  2:864255/1‑140 (MQ=255)                                                                                            tGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGca                                        <  2:788271/140‑1 (MQ=255)                                                                                            tGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGca                                        >  1:168073/1‑140 (MQ=255)                                                                                             ggCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGcac                                       <  1:1117408/140‑1 (MQ=255)                                                                                             ggCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGcac                                       <  1:646777/140‑1 (MQ=255)                                                                                                 cGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAAtt                                   >  2:950507/1‑140 (MQ=255)                                                                                                  gTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGa                                                  >  1:86445/1‑124 (MQ=255)                                                                                                  gTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGa                                                  <  2:86445/124‑1 (MQ=255)                                                                                                    ttCGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTaaa                                <  1:917466/140‑1 (MQ=255)                                                                                                      cGAATCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGATGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGGACAATTAAAAg                              <  1:762430/140‑1 (MQ=255)                                                                                                        aaTCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGcc                            >  2:747766/1‑140 (MQ=255)                                                                                                         aTCTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAGCTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCa                           >  2:680694/1‑140 (MQ=255)                                                                                                           cTTTCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaaa                         <  1:747766/140‑1 (MQ=255)                                                                                                            tttCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaaaa                        <  2:831895/140‑1 (MQ=255)                                                                                                            tttCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCaaaa                        <  2:89880/140‑1 (MQ=255)                                                                                                             ttCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAAt                       >  1:404768/1‑140 (MQ=255)                                                                                                              tCGAACGTCGTATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATg                      >  1:415129/1‑140 (MQ=255)                                                                                                                      cgtATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGa              <  2:697690/140‑1 (MQ=255)                                                                                                                        tATTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAca            <  2:918093/140‑1 (MQ=255)                                                                                                                         aTTGACCAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAcaa           <  2:168073/140‑1 (MQ=255)                                                                                                                               cAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAca            >  1:983422/1‑133 (MQ=255)                                                                                                                               cAGATGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAca            <  2:983422/133‑1 (MQ=255)                                                                                                                                  aTGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAc             >  2:787421/1‑129 (MQ=255)                                                                                                                                  aTGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGAc             <  1:787421/129‑1 (MQ=255)                                                                                                                                  aTGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAataata  >  2:468998/1‑140 (MQ=255)                                                                                                                                  aTGGAGGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAataata  >  2:1013614/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       CGCCAGCAGGCATTGATGTTACGTCATCAGGCGGCAAACTCGTCGCGCGATGTGCGTCGTCAGCTGGACAGTGGCAAACTGGATGAAGCAATGGCTCGTTTCGAATCTTTCGAACGTCGTATTGACCAGATGGAAGCGGAAGCAGAAAGCCACAGCTTCGGTAAACAAAAATCGCTGGACGATCAGTTTGCCGAACTGAAAGCCGATGATGCAATCAGCGAACAACTGGCACAATTAAAAGCCAAAATGAAGCAAGACAATCAATAATA  >  NC_000913/1368481‑1368749

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