Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	335,321	A→G	A14A (GCA→GCG)	yahD →	ankyrin repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	335,321	0	A	G	100.0%	158.8 / NA	43	A14A (GCA→GCG)	yahD	ankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/23);  total (20/23)

GGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCACAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCAGGCATTAATTGATGCCGGAG  >  NC_000913/335180‑335466                                                                                                                                              |                                                                                                                                                  ggCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGt                                                                                                                                          <  1:381950/151‑1 (MQ=255)  gCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTg                                                                                                                                         <  2:452325/151‑1 (MQ=255)     agaTACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAAGGTGat                                                                                                                                       >  1:12078/1‑150 (MQ=255)        tACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAAATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATc                                                                                                                                    <  2:140598/150‑1 (MQ=255)          cGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGAt                                                                                                                                 >  1:276938/1‑151 (MQ=255)              gAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATaa                                                                                                                               >  1:665967/1‑149 (MQ=255)                  tAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTaa                                                                                                                          >  2:329205/1‑150 (MQ=255)                   aaaTTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTaaa                                                                                                                         >  1:541271/1‑150 (MQ=255)                   aaaTTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTaaa                                                                                                                         <  2:176832/150‑1 (MQ=255)                     aTTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAAc                                                                                                                      >  1:135708/1‑151 (MQ=255)                      tttAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAAcc                                                                                                                     >  1:549663/1‑151 (MQ=255)                       ttAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCt                                                                                                                    >  2:171912/1‑151 (MQ=255)                       ttAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCt                                                                                                                    >  1:250909/1‑151 (MQ=255)                          aGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGcc                                                                                                                 >  2:543711/1‑151 (MQ=255)                             gggTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTg                                                                                                              <  2:343906/151‑1 (MQ=255)                                ttCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTgcg                                                                                                            <  1:373856/150‑1 (MQ=255)                                       aTCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGgtgt                                                                                                    <  2:618091/151‑1 (MQ=255)                                       aTCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGgtgt                                                                                                    <  1:512043/151‑1 (MQ=255)                                         cAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTc                                                                                                   <  2:342721/150‑1 (MQ=255)                                         cAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCg                                                                                                  >  1:63889/1‑151 (MQ=255)                                             ttATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGata                                                                                               <  2:18474/150‑1 (MQ=255)                                              tATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATAtt                                                                                             >  2:232403/1‑151 (MQ=255)                                                 cGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAAt                                                                                          >  1:324773/1‑151 (MQ=255)                                                 cGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAAt                                                                                          >  2:422183/1‑151 (MQ=255)                                                                   aTCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAgggg                                                                        >  1:163967/1‑151 (MQ=255)                                                                               cacaAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAAtt                                                            >  2:476875/1‑151 (MQ=255)                                                                                     aaCCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTg                                                      <  1:372359/151‑1 (MQ=255)                                                                                             gAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGttt                                              >  1:471742/1‑151 (MQ=255)                                                                                              aaTGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTta                                             <  1:263919/151‑1 (MQ=255)                                                                                                 gtgAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTtata                                           <  2:240821/150‑1 (MQ=255)                                                                                                   gAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCa                                        <  1:628659/151‑1 (MQ=255)                                                                                                     aTGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGc                                      <  1:441927/151‑1 (MQ=255)                                                                                                       gAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCa                                     >  2:5152/1‑150 (MQ=255)                                                                                                             aaaaaaTCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCt                              >  2:60164/1‑151 (MQ=255)                                                                                                              aaaaaTCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCtt                             <  2:549668/151‑1 (MQ=255)                                                                                                               aaaaTCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTg                            <  1:68619/151‑1 (MQ=255)                                                                                                               aaaaTCTACCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTAAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTg                            >  2:448701/1‑151 (MQ=255)                                                                                                                      aCCTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCa                      <  2:221642/150‑1 (MQ=255)                                                                                                                       ccTGCCGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCGATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGtt                        <  1:517040/147‑1 (MQ=255)                                                                                                                            cGACTATTTATTGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCAGGCATTa               <  2:539260/151‑1 (MQ=255)                                                                                                                                 atttattGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCAGGCATTAAtt            <  1:616114/149‑1 (MQ=255)                                                                                                                                 atttattGGCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCAGGCATTAATTGa          <  1:610111/151‑1 (MQ=255)                                                                                                                                         gCTGCGCAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCAGGCATTAATTGATGCCGGAg  <  2:135712/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  GGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATTTATTGGCTGCACAACAGGGTGATATCGATAAAGTAAAAACCTGCCTTGCGCTGGGTGTCGATATTAATACCTGCGATCGTCAGGGGAAAACGGCAATTACGCTGGCAAGTTTATATCAGCAATATGCTTGCGTTCAGGCATTAATTGATGCCGGAG  >  NC_000913/335180‑335466

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