Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	461,242	G→A	*785* (TAG→TAA)	lon →	DNA‑binding ATP‑dependent protease La

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	461,242	0	G	A	100.0%	145.5 / NA	45	*785* (TAG→TAA)	lon	DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/21);  total (24/21)

GAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGG  >  NC_000913/461123‑461373                                                                                                                        |                                                                                                                                    gagaTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                  <  2:280480/139‑1 (MQ=255)   gaTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACt                                                                                                                <  2:491134/139‑1 (MQ=255)    aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                  <  2:1395933/136‑1 (MQ=255)    aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                  >  1:1395933/1‑136 (MQ=255)    aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACta                                                                                                               >  1:1477381/1‑139 (MQ=255)    aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACta                                                                                                               >  1:811465/1‑139 (MQ=255)          gACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATaaaa                                                                                                         <  2:811465/139‑1 (MQ=255)             aaCGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACa                                                                                                      >  1:23040/1‑139 (MQ=255)             aaCGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACa                                                                                                      >  1:936634/1‑139 (MQ=255)                                  aTTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCggg                                                                                 <  2:23040/139‑1 (MQ=255)                                  aTTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCggg                                                                                 <  2:936634/139‑1 (MQ=255)                                    tCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                                               >  1:1214495/1‑123 (MQ=255)                                    tCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                                               <  2:1214495/123‑1 (MQ=255)                                    tCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCt                                                                               >  2:1096662/1‑139 (MQ=255)                                    tCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCCGGCAGGCTAATTCGGGCt                                                                               >  2:726476/1‑139 (MQ=255)                                         cTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCa                                                                          >  2:1080769/1‑139 (MQ=255)                                                gcgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCtttt                                                                   >  2:562574/1‑139 (MQ=255)                                                 cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttt                                                                  >  2:736252/1‑139 (MQ=255)                                                        aGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAACAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGtctc                                                           >  1:1193067/1‑139 (MQ=255)                                                                    ctctGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGAt                                               >  1:153279/1‑139 (MQ=255)                                                                     tctGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATg                                              <  2:985525/139‑1 (MQ=255)                                                                       tGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGc                                            >  2:675381/1‑139 (MQ=255)                                                                         gcgcTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                  >  2:392318/1‑67 (MQ=255)                                                                         gcgcTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                  <  1:392318/67‑1 (MQ=255)                                                                           gcTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                                        <  1:1413539/139‑1 (MQ=255)                                                                               cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                                        >  2:325319/1‑135 (MQ=255)                                                                               cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                                        <  1:325319/135‑1 (MQ=255)                                                                               cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCggg                                    >  2:210197/1‑139 (MQ=255)                                                                                 aaaTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGc                                                                            >  2:1079617/1‑97 (MQ=255)                                                                                 aaaTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGc                                                                            <  1:1079617/97‑1 (MQ=255)                                                                                      aaCCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGccc                             <  2:1477381/139‑1 (MQ=255)                                                                                        ccGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCtt                           <  1:1080769/139‑1 (MQ=255)                                                                                              ggTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAg                     <  1:1096662/139‑1 (MQ=255)                                                                                                tATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg                   <  1:736252/139‑1 (MQ=255)                                                                                                 aTGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCCGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAgggg                  <  1:726476/139‑1 (MQ=255)                                                                                                     aGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACTGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCg                                                          <  1:1043744/95‑1 (MQ=255)                                                                                                     aGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCg                                                          >  2:1043744/1‑95 (MQ=255)                                                                                                       gTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCt                                  <  1:792417/117‑1 (MQ=255)                                                                                                       gTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCt                                  >  2:792417/1‑117 (MQ=255)                                                                                                         tgtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTaa          >  1:1114944/1‑139 (MQ=255)                                                                                                          gtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg                   <  2:1277467/129‑1 (MQ=255)                                                                                                          gtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg                   >  1:1277467/1‑129 (MQ=255)                                                                                                               tGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAgggg      <  2:17067/137‑1 (MQ=255)                                                                                                               tGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAgggg      >  1:17067/1‑137 (MQ=255)                                                                                                                 cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtgg  >  2:1216986/1‑139 (MQ=255)                                                                                                                        |                                                                                                                                    GAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGG  >  NC_000913/461123‑461373

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: ‑