Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	988,952	C→T	G12S (GGC→AGC)	asnS ←	asparaginyl tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	988,952	0	C	T	100.0%	223.4 / NA	70	G12S (GGC→AGC)	asnS	asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/37);  total (33/37)

ATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTGGCATCTGCAATCAG  >  NC_000913/988825‑989074                                                                                                                                |                                                                                                                           atgaCAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGt                                                                                                                 >  1:99203/1‑139 (MQ=255) atgaCAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGt                                                                                                                 >  1:134206/1‑139 (MQ=255)   gaCAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCg                                                                                                               >  2:15721/1‑139 (MQ=255)      aGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCt                                                                                                            >  1:1195362/1‑139 (MQ=255)               aGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGcaca                                                                                                   <  2:176616/139‑1 (MQ=255)               aGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGcaca                                                                                                   <  2:637111/139‑1 (MQ=255)                ggATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACaa                                                                                                  <  1:622627/139‑1 (MQ=255)                 gATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAAc                                                                                                 <  1:933910/139‑1 (MQ=255)                 gATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAAc                                                                                                 <  2:1195362/139‑1 (MQ=255)                   tCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGc                                                                                               <  2:99203/139‑1 (MQ=255)                   tCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGc                                                                                               <  2:134206/139‑1 (MQ=255)                      aaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTgg                                                                                                                        >  2:307206/1‑111 (MQ=255)                      aaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTgg                                                                                                                        <  1:307206/111‑1 (MQ=255)                      aaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGcaca                                                                                                   <  2:1266214/132‑1 (MQ=255)                      aaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGcaca                                                                                                   >  1:1266214/1‑132 (MQ=255)                           ggAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCataata                                                                                       >  1:1060249/1‑139 (MQ=255)                           ggAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCataata                                                                                       >  2:125970/1‑139 (MQ=255)                            gAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAAc                                                                                                 >  2:1230192/1‑128 (MQ=255)                            gAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAAc                                                                                                 <  1:1230192/128‑1 (MQ=255)                            gAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAatat                                                                                      >  1:442427/1‑139 (MQ=255)                                 gTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctct                                                                                 <  2:691367/139‑1 (MQ=255)                                 gTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctct                                                                                 <  2:442427/139‑1 (MQ=255)                                 gTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctct                                                                                 >  2:788438/1‑139 (MQ=255)                                  tCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctctc                                                                                <  1:366602/139‑1 (MQ=255)                                    aTAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTg                                                                              <  1:788438/139‑1 (MQ=255)                                          ggCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtc                                                                                    >  2:82404/1‑127 (MQ=255)                                          ggCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtc                                                                                    <  1:82404/127‑1 (MQ=255)                                            cGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACTGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGt                                                                      >  2:907179/1‑139 (MQ=255)                                            cGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGt                                                                      >  2:1121130/1‑139 (MQ=255)                                                    gagaTGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCTCAACGCTCATAATATTCTCTCTGTTAATAGTCGGaaaaa                                                              <  1:15721/139‑1 (MQ=255)                                                           cAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCt                                                                                                            >  2:59423/1‑86 (MQ=255)                                                           cAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCt                                                                                                            <  1:59423/86‑1 (MQ=255)                                                           cAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGc                                                                                               <  2:11135/99‑1 (MQ=255)                                                           cAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGc                                                                                               >  1:11135/1‑99 (MQ=255)                                                              cTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACtt                                                    >  1:498598/1‑139 (MQ=255)                                                                 ttGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTc                                                 <  1:1121130/139‑1 (MQ=255)                                                                        tcGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGc                                                                                               <  1:757969/86‑1 (MQ=255)                                                                        tcGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGc                                                                                               >  2:757969/1‑86 (MQ=255)                                                                        tcGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctc                                                                                  <  2:500692/99‑1 (MQ=255)                                                                        tcGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctc                                                                                  >  1:500692/1‑99 (MQ=255)                                                                        tcGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTc                                                 >  1:1470575/1‑132 (MQ=255)                                                                        tcGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTc                                                 <  2:1470575/132‑1 (MQ=255)                                                                         cggcggGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTc                                                 <  2:169859/131‑1 (MQ=255)                                                                         cggcggGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTc                                                 >  1:169859/1‑131 (MQ=255)                                                                                  cgtacCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATggggg                                 <  2:1072461/138‑1 (MQ=255)                                                                                  cgtacCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATggggg                                 >  1:1072461/1‑138 (MQ=255)                                                                                         aTCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGcaca                                                                                                   >  2:570347/1‑65 (MQ=255)                                                                                         aTCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGcaca                                                                                                   <  1:570347/65‑1 (MQ=255)                                                                                          tCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCataa                                                                                         <  2:555413/74‑1 (MQ=255)                                                                                          tCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCataa                                                                                         >  1:555413/1‑74 (MQ=255)                                                                                          tCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATTTTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTAt                        >  2:695373/1‑139 (MQ=255)                                                                                          tCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATaaa                                                         <  2:259332/106‑1 (MQ=255)                                                                                          tCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATaaa                                                         >  1:259332/1‑106 (MQ=255)                                                                                             gcgcACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGtt                     >  1:237948/1‑139 (MQ=255)                                                                                              cgcACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctc                                                                                  >  1:589019/1‑77 (MQ=255)                                                                                              cgcACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtctc                                                                                  <  2:589019/77‑1 (MQ=255)                                                                                               gcACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTAc                   <  2:237948/139‑1 (MQ=255)                                                                                                cACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACg                                                                                                >  1:79853/1‑61 (MQ=255)                                                                                                cACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACg                                                                                                <  2:79853/61‑1 (MQ=255)                                                                                                cACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGaa                                                                 >  1:1289806/1‑92 (MQ=255)                                                                                                cACGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGaa                                                                 <  2:1289806/92‑1 (MQ=255)                                                                                                 aCGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCataat                                                                                        >  2:1365210/1‑68 (MQ=255)                                                                                                 aCGGTGACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCataat                                                                                        <  1:1365210/68‑1 (MQ=255)                                                                                                      gACTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTGGCAt            <  1:1146279/139‑1 (MQ=255)                                                                                                       aCTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTgg               <  1:511472/135‑1 (MQ=255)                                                                                                       aCTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTgg               >  2:511472/1‑135 (MQ=255)                                                                                                       aCTTCGCTGTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTGGCATc           <  2:498598/139‑1 (MQ=255)                                                                                                               gTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtc                                                                                    <  1:1415795/58‑1 (MQ=255)                                                                                                               gTCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATtc                                                                                    >  2:1415795/1‑58 (MQ=255)                                                                                                                tCAACGGCTACACGGCTCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTGGCATCTGCAATCAg  <  1:907179/139‑1 (MQ=255)                                                                                                                                |                                                                                                                           ATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACGCTCATAATATTCTCTCTGTTAATAGTCGGAAAAAATAAACACTTGTCCACCCGTAAATGGGGGGATACCTATGTTACCTGGCATCTGCAATCAG  >  NC_000913/988825‑989074

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