Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,528,926	A→G	R417R (CGT→CGC)	ligA ←	DNA ligase, NAD(+)‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,528,926	0	A	G	100.0%	179.0 / NA	53	R417R (CGT→CGC)	ligA	DNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/30);  total (23/30)

TTTTTCAACCAGCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCG  >  NC_000913/2528782‑2529069                                                                                                                                                 |                                                                                                                                                tttttCAACCAGCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTc                                                                                                                                            >  1:675646/1‑150 (MQ=255)     tCAACCAGCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACg                                                                                                                                        <  2:370500/150‑1 (MQ=255)          cAGCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGTAAACATAGTGTTTCATCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCCGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAg                                                                                                                                    >  1:522227/1‑149 (MQ=255)           aGCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAAc                                                                                                                                 >  2:771788/1‑151 (MQ=255)            gCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAAcc                                                                                                                                >  1:721268/1‑151 (MQ=255)             cTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAAcc                                                                                                                                <  2:384559/150‑1 (MQ=255)               gATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACca                                                                                                                               <  2:844020/149‑1 (MQ=255)                aTCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACAt                                                                                                                            <  1:745643/151‑1 (MQ=255)                aTCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACAt                                                                                                                            <  2:721311/151‑1 (MQ=255)                   gctgatTTTGTCGCCCATTCCGTCAACATCCATCGCACAGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATAc                                                                                                                          >  1:448784/3‑150 (MQ=255)                   gatgatTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATAcc                                                                                                                         >  1:787455/1‑151 (MQ=255)                    atgatTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCg                                                                                                                        <  2:624236/151‑1 (MQ=255)                    atgatTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCg                                                                                                                        <  1:74531/151‑1 (MQ=255)                      gatTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGa                                                                                                                      <  2:288778/151‑1 (MQ=255)                           tGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATg                                                                                                                 >  2:212318/1‑151 (MQ=255)                               gCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTc                                                                                                             >  2:206514/1‑151 (MQ=255)                                cccATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCggc                                                                                                             >  2:635700/1‑148 (MQ=255)                                 ccATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCgg                                                                                                           >  1:200594/1‑151 (MQ=255)                                     tCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAAt                                                                                                       >  1:598477/1‑151 (MQ=255)                                      ccGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATa                                                                                                      <  2:630895/151‑1 (MQ=255)                                       cGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATAc                                                                                                     <  1:374016/151‑1 (MQ=255)                                        gTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATAc                                                                                                     >  1:753232/1‑150 (MQ=255)                                           aaCATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAAcc                                                                                                 >  1:739024/1‑151 (MQ=255)                                            aCATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCt                                                                                                <  1:568069/151‑1 (MQ=255)                                              atccatcGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCa                                                                                              <  2:50241/151‑1 (MQ=255)                                                ccatcGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCAc                                                                                             <  2:465169/150‑1 (MQ=255)                                                    cGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGt                                                                                         <  1:802690/150‑1 (MQ=255)                                                    cGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGt                                                                                         <  1:63693/150‑1 (MQ=255)                                                     gCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTAt                                                                                       >  2:303447/1‑151 (MQ=255)                                                       aCGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATc                                                                                      <  2:812622/150‑1 (MQ=255)                                                       aCGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCt                                                                                     <  1:102366/151‑1 (MQ=255)                                                         ggcggGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTc                                                                                   >  1:749163/1‑151 (MQ=255)                                                             ggAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCggg                                                                               >  1:413475/1‑151 (MQ=255)                                                               aaaCAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCg                                                                             >  2:22389/1‑151 (MQ=255)                                                                     aGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGa                                                                       <  2:836314/151‑1 (MQ=255)                                                                     aGTGTTTCAGCGACTCTGTACGCTGCGCGCCGCAAATCAGGCCACCGGTTCAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATCCCGGACAATGCGGCGGGAATACAACCTCACGGGTATGTTCCGGGCGTTCa                                                                         <  1:245409/149‑1 (MQ=255)                                                                      gtgtTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGaa                                                                      <  1:150783/151‑1 (MQ=255)                                                                        gtTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGaaa                                                                     <  2:84765/150‑1 (MQ=255)                                                                          ttCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGc                                                                   <  2:10515/150‑1 (MQ=255)                                                                            cAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCac                                                                 <  2:407534/150‑1 (MQ=255)                                                                                        aCGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAaccacc                                                    >  2:174235/1‑151 (MQ=255)                                                                                             gcgcACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCGCCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGTGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCg                                                <  1:657363/150‑1 (MQ=255)                                                                                                cACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGAt                                            >  1:79022/1‑151 (MQ=255)                                                                                                    gCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACg                                        >  1:146912/1‑151 (MQ=255)                                                                                                    gCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACg                                        >  1:375694/1‑151 (MQ=255)                                                                                                           aGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGcc                                   >  2:476133/1‑149 (MQ=255)                                                                                                                                aCCGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTAt            >  2:730653/1‑151 (MQ=255)                                                                                                                                 ccGCTTCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATc           <  1:171955/151‑1 (MQ=255)                                                                                                                                 ccGCTTCACCGTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATc           <  1:745050/151‑1 (MQ=255)                                                                                                                                      tCACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGcc        <  1:843081/149‑1 (MQ=255)                                                                                                                                       cACCTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAAt     <  2:195689/151‑1 (MQ=255)                                                                                                                                         ccTTCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATg    <  2:756918/150‑1 (MQ=255)                                                                                                                                            tCCACGCGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCg  <  2:119034/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                TTTTTCAACCAGCTGATCGATGATTTTGTCGCCCATTCCGTCAACATCCATCGCACGGCGGGAAACAAAGTGTTTCAGCGACTCTTTACGCTGCGCACCGCAAATCAGGCCACCGGTACAGCGGGCAACCGCTTCACCTTCCACACGCTCAACGTCAGAACCACATACCGGACAATGCGTCGGGAATACAACCTCACGGGTATCTTCCGGGCGTTCAGAAAGCACGACGTTAACCACCTGCGGGATCACGTCGCCAGCGCGGCGAATCACCACTTTATCGCCAATGCG  >  NC_000913/2528782‑2529069

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