Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,413,542	A→T	Y67N (TAT→AAT)	ydaC ←	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,542	0	A	T	100.0%	184.3 / NA	56	Y67N (TAT→AAT)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (25/31);  total (25/31)

AAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTT  >  NC_000913/1413399‑1413684                                                                                                                                                |                                                                                                                                               aaGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATCATCAGCTATGAAGTTTTAATTTGAATTCAATg                                                                                                                                           <  1:975026/149‑1 (MQ=255)   gCATCAGGCCGTATTCAACGATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCa                                                                                                                                         >  2:147098/1‑149 (MQ=255)           ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                 <  1:828763/149‑1 (MQ=255)                     ccATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGt                                                                                                                       <  1:438534/149‑1 (MQ=255)                        tCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTTACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtc                                                                                                                    <  2:169001/149‑1 (MQ=255)                           gcgctTCATTAAAGAGTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCg                                                                                                                 <  1:1132193/146‑1 (MQ=255)                             tctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCa                                                                                                               >  2:778106/1‑149 (MQ=255)                              ctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCaa                                                                                                              <  2:1015634/149‑1 (MQ=255)                              ctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCaa                                                                                                              <  1:632876/149‑1 (MQ=255)                              ctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCaa                                                                                                              <  1:632855/149‑1 (MQ=255)                               ttCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAt                                                                                                             >  1:151144/1‑149 (MQ=255)                                   ttAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaa                                                                                                         >  2:1152900/1‑149 (MQ=255)                                   ttAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaa                                                                                                         >  1:512245/1‑149 (MQ=255)                                   ttAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaa                                                                                                         >  1:512236/1‑149 (MQ=255)                                   ttAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaa                                                                                                         >  1:400129/1‑149 (MQ=255)                                        gATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTa                                                                                                    <  1:546575/149‑1 (MQ=255)                                              ttGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGCTTCTCGCAATTAAAATTTATCAGtt                                                                                              <  1:646604/149‑1 (MQ=255)                                               ggggCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATTTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGttt                                                                                             <  1:115750/146‑1 (MQ=255)                                                   ccATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATCATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGACGTTTCTCGCAATTAAAATTTATCCGTTTTACt                                                                                         >  2:573107/1‑149 (MQ=255)                                                   ccATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACt                                                                                         <  2:1142003/149‑1 (MQ=255)                                                   ccATCGAAGCGCCTCGGGCAGGCGGCAAGGTTAGGCTGAACTTGACGCCCGATGTGGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGTCTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACt                                                                                         <  2:860962/149‑1 (MQ=255)                                                     aTCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGACATTGACGCCTGATGTTGATTCTTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGCGGACGGAAGTTTCTCGCAATACAAATTTATCAGGTTTagtat                                                                                       >  1:46760/1‑145 (MQ=255)                                                       cGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCt                                                                                     <  2:61113/149‑1 (MQ=255)                                                       cGCATCACCTCTGGCAGGCGCCAACGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCt                                                                                     <  1:1236241/149‑1 (MQ=255)                                                            cACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTGCTTctct                                                                                <  2:827213/149‑1 (MQ=255)                                                              ccTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTGTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctct                                                                              >  1:857711/1‑149 (MQ=255)                                                               ctctGGCAGGCGCCAATGTTAGACTGAAATTTGCGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                             <  1:247951/149‑1 (MQ=255)                                                               ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGGTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                             >  1:1240562/1‑149 (MQ=255)                                                                tctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTgg                                                                            >  2:509308/1‑149 (MQ=255)                                                                tcgggCAGGTGACAATGTTAGACGGAAATTGACGTCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTgg                                                                            <  2:372302/146‑1 (MQ=255)                                                                   ggCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGaaa                                                                         <  2:1213220/149‑1 (MQ=255)                                                                    gCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTGTCTCGCGATTAAAATTTATCAGTTGTACTTTCTGCTCTCTGGCAAc                                                                        >  1:639361/1‑149 (MQ=255)                                                                      aGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGc                                                                      >  2:630847/1‑149 (MQ=255)                                                                               tGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGCAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGcttctt                                                             <  1:58134/149‑1 (MQ=255)                                                                                 ttAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCtttt                                                           >  2:897080/1‑149 (MQ=255)                                                                                 ttAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAGTAAAATTTATCAATATTACTTTCTACTCTCTGGAACCGCCTGCTTAtttt                                                           >  1:168384/1‑149 (MQ=255)                                                                                    gACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTa                                                        <  1:616466/149‑1 (MQ=255)                                                                                       tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCt                                                     <  1:597569/149‑1 (MQ=255)                                                                                               aCGCCCGATGTTAATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGCCTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAAGCCCTGCTTCTTTTTTTCCTGAGAACAt                                             >  1:1006882/1‑149 (MQ=255)                                                                                                    cGATGTTGATTATTAATAATCAGTTATTAAGTTTTAATTTTTATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAtttttt                                        <  1:117512/149‑1 (MQ=255)                                                                                                    cGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGATCAtttttt                                        >  1:152736/1‑149 (MQ=255)                                                                                                     taTGTTGATTATTAAGAATCAGCTATGAAGTTTTAATGTGACTTTAATTTAATTCTCTTGGATTGAATTTGCTCGCAATTAAAATTTATCAGTTTTACTTTCTTCTCTCTGGAAACGCCTGTTTCTTTTGTACCTTAGAGCATTTTTTc                                       <  2:448458/148‑1 (MQ=255)                                                                                                          tGATTATTAATAATCAGCTATGAAGTTTCAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCAtt                                  >  2:632369/1‑149 (MQ=255)                                                                                                            attattAACAATCAGCTATTAAGTTTTCATTTGAATTCCATGCAATTCTCGAGGACTGAAGGTTCTCCCAATTAAAATTTACCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCt                                <  1:85463/149‑1 (MQ=255)                                                                                                              tattaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGa                              >  1:387840/1‑149 (MQ=255)                                                                                                                    taatCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGt                        <  1:934346/149‑1 (MQ=255)                                                                                                                    taatCAGCTATGAAGTTTTAATGTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAACCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGt                        <  2:322813/149‑1 (MQ=255)                                                                                                                      atCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCGTTTTTTCGCATTCTGATTGCGttc                      >  2:441991/1‑148 (MQ=255)                                                                                                                                 cgTTTTAATTTGTATTCAATGTAATTCTCGAGGACTGTAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAtttt           <  2:52489/148‑1 (MQ=255)                                                                                                                                  gTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGATTCTTTTATACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTg          <  2:782535/149‑1 (MQ=255)                                                                                                                                   ttttATTTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGa         <  2:367379/149‑1 (MQ=255)                                                                                                                                   ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAACTTTCTCAGTTTTACTTTCTGCTCTCTGGAAACGCATGCTTCTCTTTTAACTGAGAGCGTTTTTTCGCATTCTGAGTTCGTTAGTTTAGAttttca         >  2:81569/1‑147 (MQ=255)                                                                                                                                   ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGa         >  1:1198823/1‑149 (MQ=255)                                                                                                                                     ttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAat       <  1:509308/149‑1 (MQ=255)                                                                                                                                        aTTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATc    >  2:254100/1‑149 (MQ=255)                                                                                                                                          ttGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCtt  >  1:1150050/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                               AAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTT  >  NC_000913/1413399‑1413684

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