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%	206.4 / NA	63	*785* (TAG→TAA)	lon	DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/31);  total (32/31)

TCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGG  >  NC_000913/461116‑461357                                                                                                                               |                                                                                                                    tCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCa                                                                                                         <  1:750182/139‑1 (MQ=255)      aaGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATg                                                                                                    <  2:1360120/139‑1 (MQ=255)           aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACta                                                                                               >  2:292551/1‑139 (MQ=255)           aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACta                                                                                               >  1:1678610/1‑139 (MQ=255)                      cGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAgg                                                                                    >  1:342400/1‑139 (MQ=255)                           ttGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                               <  1:292551/139‑1 (MQ=255)                             gCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCa                                                                             >  1:234378/1‑139 (MQ=255)                                 aTCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATCATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCt                                                                         >  1:134011/1‑139 (MQ=255)                                 aTCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCt                                                                         >  1:1193190/1‑139 (MQ=255)                                   cTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTaa                                                                       >  2:429741/1‑139 (MQ=255)                                    tGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAAt                                                                      <  2:1193190/139‑1 (MQ=255)                                           tCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCt                                                               >  2:1529311/1‑139 (MQ=255)                                             aTCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTg                                                             >  2:1207653/1‑139 (MQ=255)                                             aTCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTg                                                             >  2:903217/1‑139 (MQ=255)                                              tCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGc                                                            <  2:1678610/139‑1 (MQ=255)                                                cTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTg                                                             >  2:665983/1‑136 (MQ=255)                                                cTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTg                                                             <  1:665983/136‑1 (MQ=255)                                                   tgAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGAc                                                                                                                 >  1:1130302/1‑81 (MQ=255)                                                   tgAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGcc                                                       >  1:61831/1‑139 (MQ=255)                                                   tgAAGCGCATTGAGGAAGCTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGAc                                                                                                                 <  2:1130302/81‑1 (MQ=255)                                                     aaGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                  <  2:1379526/94‑1 (MQ=255)                                                     aaGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                  >  1:1379526/1‑94 (MQ=255)                                                       gcgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCtttt                                                   >  2:1620078/1‑139 (MQ=255)                                                        cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCaaaa                                                                                                      >  1:1490032/1‑87 (MQ=255)                                                        cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCaaaa                                                                                                      <  2:1490032/87‑1 (MQ=255)                                                        cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCaaa                                                                                                       >  1:480133/1‑86 (MQ=255)                                                        cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCaaa                                                                                                       <  2:480133/86‑1 (MQ=255)                                                        cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttt                                                  >  2:1413325/1‑139 (MQ=255)                                                        cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCGAATTCGGGCTTGCCAGCCttttt                                                  >  2:1626855/1‑139 (MQ=255)                                                         gcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCtttttt                                                 >  2:830899/1‑139 (MQ=255)                                                         gcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCtttttt                                                 >  1:21538/1‑139 (MQ=255)                                                          cATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttt                                                >  1:712721/1‑139 (MQ=255)                                                            ttGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGt                                              <  1:1620078/139‑1 (MQ=255)                                                            ttGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGt                                              <  2:61831/139‑1 (MQ=255)                                                            ttGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGt                                              <  2:234378/139‑1 (MQ=255)                                                             tGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGtc                                             >  2:369097/1‑139 (MQ=255)                                                                       cTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGtt                                   <  2:68786/139‑1 (MQ=255)                                                                         gACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAg                                 <  1:1529311/139‑1 (MQ=255)                                                                           ctctGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGAt                               <  1:1413325/139‑1 (MQ=255)                                                                            tctGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATg                              <  1:830899/139‑1 (MQ=255)                                                                               ggCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACt                                                                                                >  1:1048306/1‑70 (MQ=255)                                                                               ggCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACt                                                                                                <  2:1048306/70‑1 (MQ=255)                                                                               ggCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTCTTTGTCTCGCTAAGTTAGATGGCg                           >  2:855426/1‑139 (MQ=255)                                                                                  gcTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                        <  1:429741/139‑1 (MQ=255)                                                                                      cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGTCCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCggg                    >  2:1226656/1‑139 (MQ=255)                                                                                      cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAAt                                                                      >  2:917289/1‑89 (MQ=255)                                                                                      cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAAt                                                                      <  1:917289/89‑1 (MQ=255)                                                                                        aaaTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCt                  <  2:344472/139‑1 (MQ=255)                                                                                          aTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTg                <  2:21538/139‑1 (MQ=255)                                                                                          aTGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTg                <  2:134011/139‑1 (MQ=255)                                                                                              aCCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                  <  2:316784/53‑1 (MQ=255)                                                                                              aCCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                  >  1:316784/1‑53 (MQ=255)                                                                                                      gTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAgg    <  1:369097/139‑1 (MQ=255)                                                                                                       tATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTCTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg   <  1:855426/139‑1 (MQ=255)                                                                                                        aTGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAgggg  <  1:903217/139‑1 (MQ=255)                                                                                                                tgtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCg                                          >  2:1541746/1‑91 (MQ=255)                                                                                                                tgtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCg                                          <  1:1541746/91‑1 (MQ=255)                                                                                                                   gACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAACCAGGGCTGGCAGGCTAATTCGGGCTGGCCAGcc                                                       <  2:888795/75‑1 (MQ=255)                                                                                                                   gACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGcc                                                       >  1:888795/1‑75 (MQ=255)                                                                                                                        cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGc                                                            >  2:703186/1‑65 (MQ=255)                                                                                                                        cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGc                                                            <  1:703186/65‑1 (MQ=255)                                                                                                                        cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCt                                        <  2:1271674/85‑1 (MQ=255)                                                                                                                        cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCt                                        >  1:1271674/1‑85 (MQ=255)                                                                                                                               |                                                                                                                    TCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGG  >  NC_000913/461116‑461357

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