Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,653,652	G→A	S88S (AGC→AGT)	rspB ←	putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,653,652	0	G	A	100.0%	212.7 / NA	66	S88S (AGC→AGT)	rspB	putative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/30);  total (36/30)

TTTCCACGCATTTTTTGCCGGAACCACGGCATATTCACTGAAACCACCGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAGCTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACGATAA  >  NC_000913/1653526‑1653787                                                                                                                               |                                                                                                                                        tttCCACGCATTTTTTGCCGGAACCACGGCATATTCACTGAAACCACCGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGa                                                                                                                             >  2:176148/1‑139 (MQ=255)         cATTTTTTGCCGGAACCACGGCATATTCACTGAAACCACCGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGc                                                                                                                     <  2:342018/139‑1 (MQ=255)                               atatTCACTGAAACCACCGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGgcg                                                                                               <  1:176148/139‑1 (MQ=255)                                      cTGAAACCACCGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCa                                                                                        <  2:910536/139‑1 (MQ=255)                                           accaccGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCt                                                                                   >  1:571360/1‑139 (MQ=255)                                              accGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCa                                                                                <  1:232139/139‑1 (MQ=255)                                                 gTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCAccc                                                                             <  2:709810/139‑1 (MQ=255)                                                                tAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCa                                                                                        >  2:254711/1‑113 (MQ=255)                                                                tAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCa                                                                                        <  1:254711/113‑1 (MQ=255)                                                                      aGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAg                                                        <  2:571360/139‑1 (MQ=255)                                                                       gCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAaca                                                                                                                 <  1:205788/81‑1 (MQ=255)                                                                       gCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAaca                                                                                                                 >  2:205788/1‑81 (MQ=255)                                                                       gCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGa                                                       >  1:186452/1‑139 (MQ=255)                                                                         cAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTcc                                                                                         <  2:619015/103‑1 (MQ=255)                                                                         cAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTcc                                                                                         >  1:619015/1‑103 (MQ=255)                                                                         cAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAAt                                                     >  2:709127/1‑139 (MQ=255)                                                                         cAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAAt                                                     >  2:258946/1‑139 (MQ=255)                                                                         cAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCACACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCAGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAAt                                                     >  2:366556/1‑139 (MQ=255)                                                                           gtgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGc                                                                                                                     >  1:789424/1‑73 (MQ=255)                                                                           gtgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGc                                                                                                                     <  2:789424/73‑1 (MQ=255)                                                                            tgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCAcc                                                                                                         >  1:249969/1‑84 (MQ=255)                                                                            tgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCAcc                                                                                                         <  2:249969/84‑1 (MQ=255)                                                                            tgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCa                                                  >  1:534652/1‑139 (MQ=255)                                                                            tgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCa                                                  >  2:324425/1‑139 (MQ=255)                                                                            tgtCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCa                                                  >  2:685971/1‑139 (MQ=255)                                                                                gTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGAc                                              >  2:750787/1‑139 (MQ=255)                                                                                 tACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGAcc                                             >  1:1435/1‑139 (MQ=255)                                                                                    aaaCGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAAt                                          <  2:899468/139‑1 (MQ=255)                                                                                      taagagacagtACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGa                                        <  1:549354/129‑1 (MQ=255)                                                                                                tACCTATAGAGCACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGata                              <  1:324425/139‑1 (MQ=255)                                                                                                           cACGGATAGCATTGTCCACAACTGACCACCGGATCGACAGGAACGCGGTCAGCGACTCTGGCGCTTGTCACCCCTTTTCCCACTTCATCAATTACGCCAGAGAATTCATGACCAATGACGCGCGGatatgtgtctctta                   >  2:549354/1‑129 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCAGCGGATAGACAGCAACACGTTGACCGACTCTGGGGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAg                   >  2:600634/1‑139 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCAc                                                                                                          >  1:475216/1‑52 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCAc                                                                                                          <  2:475216/52‑1 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCAcc                                                                              <  1:140805/80‑1 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCAcc                                                                              <  1:441849/80‑1 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCAcc                                                                              >  2:441849/1‑80 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCAcc                                                                              >  2:140805/1‑80 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAg                   >  1:731154/1‑139 (MQ=255)                                                                                                           cACGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAg                   >  2:661618/1‑139 (MQ=255)                                                                                                            aCGGATAGCTATCCCCACAACTGACCACAGGATCGACAGCAAGACGTTTGCCGACTCTGGCGCTTTCTACGCCTTGAACCACTGCATCAATGACGAGAAAGAATTCATGACCAATTACGAGCGGATATTTCGCAAAAgg                  >  2:258820/1‑139 (MQ=255)                                                                                                            aCGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAgg                  >  2:197092/1‑139 (MQ=255)                                                                                                             cGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGAct                                                                                                     <  1:848735/55‑1 (MQ=255)                                                                                                             cGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGAct                                                                                                     >  2:848735/1‑55 (MQ=255)                                                                                                             cGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGgc                                                                                                <  1:856398/60‑1 (MQ=255)                                                                                                             cGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGgc                                                                                                >  2:856398/1‑60 (MQ=255)                                                                                                             cGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGAcg                                      <  1:772211/118‑1 (MQ=255)                                                                                                             cGGATAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGAcg                                      >  2:772211/1‑118 (MQ=255)                                                                                                                aTAGCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTa              >  2:354239/1‑139 (MQ=255)                                                                                                                   gCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGAcg                                      >  2:674919/1‑112 (MQ=255)                                                                                                                   gCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGAcg                                      <  1:674919/112‑1 (MQ=255)                                                                                                                   gCAATGCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGc           >  1:469608/1‑139 (MQ=255)                                                                                                                        gCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACg                                                                                      >  1:70771/1‑59 (MQ=255)                                                                                                                        gCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACg                                                                                      <  2:70771/59‑1 (MQ=255)                                                                                                                        gCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACg      <  1:354239/139‑1 (MQ=255)                                                                                                                        gCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACg      <  2:1435/139‑1 (MQ=255)                                                                                                                        gCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACg      <  1:197092/139‑1 (MQ=255)                                                                                                                        gCCCACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACg      <  1:750787/139‑1 (MQ=255)                                                                                                                         cccACAACTGACCACCGGATCGACAGCAACACGTTCAcc                                                                                                         <  1:904596/39‑1 (MQ=255)                                                                                                                         cccACAACTGACCACCGGATCGACAGCAACACGTTCAcc                                                                                                         >  2:904596/1‑39 (MQ=255)                                                                                                                         cccACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCAt                                                 >  2:822799/1‑95 (MQ=255)                                                                                                                         cccACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCAt                                                 <  1:822799/95‑1 (MQ=255)                                                                                                                         cccACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACGa     >  2:21758/1‑139 (MQ=255)                                                                                                                          ccACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCa                                                                                        <  2:317474/55‑1 (MQ=255)                                                                                                                          ccACAACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCa                                                                                        >  1:317474/1‑55 (MQ=255)                                                                                                                            acaACTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACGATaa  <  1:661618/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                        TTTCCACGCATTTTTTGCCGGAACCACGGCATATTCACTGAAACCACCGTCAGCGTGCACACCTAATACAGCCAGTGTCGTACAAACGTTCGGTTTACCTATAGAGCACGGATAGCAATGCCCACAGCTGACCACCGGATCGACAGCAACACGTTCACCGACTCTGGCGCTTTCCACGCCTTCACCCACTGCATCAATGACGCCAAAGAATTCATGACCAATGACGCGCGGATATTTCGCAAAAGGATTATGCCCACGATAA  >  NC_000913/1653526‑1653787

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