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

ATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGC  >  NC_000913/461126‑461374                                                                                                                     |                                                                                                                                     aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGc                                                                                                                    <  2:76331/135‑1 (MQ=255) aTTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGc                                                                                                                    >  1:76331/1‑135 (MQ=255)   tCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                   <  2:388801/134‑1 (MQ=255)   tCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCa                                                                                                                   >  1:388801/1‑134 (MQ=255)      tGACAACGTAATTGCCGAACTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATaaa                                                                                                           >  1:1019285/1‑139 (MQ=255)          aaCGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACa                                                                                                       >  1:1304203/1‑139 (MQ=255)           aCGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAg                                                                                                      >  2:1453070/1‑139 (MQ=255)            cGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAgg                                                                                                     >  1:1279393/1‑139 (MQ=255)                 ttGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                                                <  2:1304203/139‑1 (MQ=255)                 ttGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                                                <  1:802279/139‑1 (MQ=255)                       aTCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCt                                                                                          >  2:772038/1‑139 (MQ=255)                          tGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAAt                                                                                       <  1:1453070/139‑1 (MQ=255)                            gACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTc                                                                                     <  2:1019285/139‑1 (MQ=255)                               aTTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCggg                                                                                  <  2:1279393/139‑1 (MQ=255)                                      cTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCa                                                                           >  2:1808060/1‑139 (MQ=255)                                       tgtgAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                                                >  2:1487486/1‑117 (MQ=255)                                       tgtgAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTgg                                                                                                <  1:1487486/117‑1 (MQ=255)                                           aaGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCtt                                                                      >  1:289093/1‑139 (MQ=255)                                              cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATaaaaa                                                                                                         <  2:647548/101‑1 (MQ=255)                                              cgcATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATaaaaa                                                                                                         >  1:647548/1‑101 (MQ=255)                                                cATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAgg                                                                                            >  2:785475/1‑112 (MQ=255)                                                cATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAgg                                                                                            <  1:785475/112‑1 (MQ=255)                                                cATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCttttttt                                                                 >  1:901365/1‑139 (MQ=255)                                                  ttGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGt                                                               <  2:843365/139‑1 (MQ=255)                                                   tGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGtc                                                              >  1:1599321/1‑139 (MQ=255)                                                              tGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCCTTTTTTGTCTCGCTAAGTTa                                                   >  1:359462/1‑139 (MQ=255)                                                               gACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAg                                                  >  1:1355557/1‑139 (MQ=255)                                                                aCTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACtaat                                                                                                              >  1:818214/1‑78 (MQ=255)                                                                aCTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACtaat                                                                                                              <  2:818214/78‑1 (MQ=255)                                                                   ctGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATgg                                              <  2:289093/139‑1 (MQ=255)                                                                     ggCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACtaat                                                                                                              <  2:248950/73‑1 (MQ=255)                                                                     ggCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACtaat                                                                                                              >  1:248950/1‑73 (MQ=255)                                                                         cTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATc                                        <  2:1355557/139‑1 (MQ=255)                                                                            cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCgg                                                                                   <  1:1605976/93‑1 (MQ=255)                                                                            cAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCgg                                                                                   >  2:1605976/1‑93 (MQ=255)                                                                                 tGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGc                                <  2:901365/139‑1 (MQ=255)                                                                                  gAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGcc                               <  2:1599321/139‑1 (MQ=255)                                                                                  gAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCCTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGcc                               <  2:359462/139‑1 (MQ=255)                                                                                   aaCCGTCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGccc                              <  2:1651490/139‑1 (MQ=255)                                                                                        tCTGGTATGCAGGTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAACTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCttatt                         <  1:1808060/139‑1 (MQ=255)                                                                                                   ggTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGcc                                                                            >  1:716121/1‑77 (MQ=255)                                                                                                   ggTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGcc                                                                            <  2:716121/77‑1 (MQ=255)                                                                                                    gTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTg                                                                              >  1:76353/1‑74 (MQ=255)                                                                                                    gTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTg                                                                              <  2:76353/74‑1 (MQ=255)                                                                                                    gTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgtg                 <  2:1166168/135‑1 (MQ=255)                                                                                                    gTTGTGACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgtg                 >  1:1166168/1‑135 (MQ=255)                                                                                                      tgtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgt                  <  1:527946/132‑1 (MQ=255)                                                                                                      tgtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGgt                  >  2:527946/1‑132 (MQ=255)                                                                                                      tgtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTaa           >  2:261387/1‑139 (MQ=255)                                                                                                       gtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg                    <  2:1805179/129‑1 (MQ=255)                                                                                                       gtgACTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAggg                    >  1:1805179/1‑129 (MQ=255)                                                                                                          aCTGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAgggg       >  1:427730/1‑139 (MQ=255)                                                                                                            tGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                                         <  1:1145714/103‑1 (MQ=255)                                                                                                            tGCAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGAt                                         >  2:1145714/1‑103 (MQ=255)                                                                                                              cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGcc                               <  2:1798370/111‑1 (MQ=255)                                                                                                              cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGcc                               >  1:1798370/1‑111 (MQ=255)                                                                                                              cAAAATAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtgg   >  2:996916/1‑139 (MQ=255)                                                                                                               aaaaTAATGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGAtggc  <  1:996916/139‑1 (MQ=255)                                                                                                                     |                                                                                                                                     ATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGCAGGCTAATTCGGGCTTGCCAGCCTTTTTTTGTCTCGCTAAGTTAGATGGCGGATCGGGCTTGCCCTTATTAAGGGGTGTTGTAAGGGGATGGC  >  NC_000913/461126‑461374

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