Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	459,834	C→A	A316E (GCG→GAG)	lon →	Lon protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	459,834	0	C	A	100.0%	82.6 / NA	33	A316E (GCG→GAG)	lon	Lon protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/16);  total (17/16)

GTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCC  >  NC_000913/459709‑459957                                                                                                                              |                                                                                                                            gTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCCCATCAAGGTCCAAAAAGACCTGCGTCAGGAGCa                                                                                                                          >  2:110718/1‑129 (MQ=255)       gAAGCTGAAAATTATGTCTCCGATCTCGACAGAAGCGACCGTAGAACAAGGTTATAACGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCCAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAAt                                                                                                                    >  1:188608/1‑129 (MQ=255)            tGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCAGTGGAATGCGCGTAGCAAGGTAAAAAAAGACCTGCGTCAGGAGCAGGGAAACCtag                                                                                                               >  2:103202/1‑127 (MQ=255)            tGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGATTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCCTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTg                                                                                                               >  1:492550/1‑129 (MQ=255)             gtatatgatgTCTCCGATGTCGGCAGAACCTGCCGTAGTGCGCTGTTAGATGGACTGTATGGTACAGGTGCCGTGGGATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGa                                                                                                              <  2:84574/125‑1 (MQ=255)              aaaaTGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGTCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGAt                                                                                                             <  1:400025/129‑1 (MQ=255)                 atgatgTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGAATGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATAcc                                                                                                          >  1:403539/1‑129 (MQ=255)                  tgatgTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGATACAGGTGCCGTGGAATGCGCGTAGCAAGGTAAAAAAAGACCAGCTTCAGGAGCAGGAAATCCTTGATACCg                                                                                                         >  1:417265/1‑129 (MQ=255)                        ctcCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGTGCGTATTAATGTTAATATAGATTTTCGTTAGGAGCAGGAAAACCTTGAGACCGACCAtt                                                                                                   <  1:27209/129‑1 (MQ=255)                         tcCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTa                                                                                                  >  2:563558/1‑129 (MQ=255)                                  gggaGAAGCGAGCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAAGCCTTGATACCGACCATTATGGTCTGGa                                                                                         <  2:66864/126‑1 (MQ=255)                                  ggCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGTGTTAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGa                                                                                         <  1:652644/129‑1 (MQ=255)                                       aaGCGACCGTAGTACGTGGTTATATCGACTGAATGGTACACGTGCCGTGCAAAGCACGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCACTATGGTCCGGAgaacg                                                                                    >  2:524438/1‑125 (MQ=255)                                            acggTACTGCGTGTGTTTTTCGACTGGATGGTACAGGTGCCGTGGAATTCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGaaa                                                                               <  2:269047/126‑1 (MQ=255)                                            aCCGCAGTTCCTTGTAATATAGACTGGATGCTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGCAAATCCTTGATACCTACCATTATGGTCTGGAGCGCGTGaaa                                                                               <  2:269165/129‑1 (MQ=255)                                                          ttATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATTCCGTTCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGa                                                                 <  1:603358/129‑1 (MQ=255)                                                           tataTCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAATACCTGCGTCAGGAGCAGGAAATTTTTTATACTGACTATTATGGTCTGGAGCGCGTGAAAGATCGAAGCCTTGAg                                                                <  1:110718/129‑1 (MQ=255)                                                             tatCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTa                                                              <  2:325701/129‑1 (MQ=255)                                                                  aCTGGATGGTACAGGTGCTGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGATCTTGATTTCTTTGATTCTGTTTATTATGGTTTGGAGGGCGTGAAAGATCGAATCCGGGAGTATCTTg                                                         <  1:225326/129‑1 (MQ=255)                                                                   cTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGc                                                        >  2:9976/1‑129 (MQ=255)                                                                           tACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACTGACCATTATGGTTTTGAGCGCGTGAAAGATCGAATCCTTGAGTATCTGGCGGTTCaaa                                                <  1:501302/129‑1 (MQ=255)                                                                                 ggCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACTTATTTTTATTGTTTTTTGTGCTTGAAAAATCGAATCCATGAGTATCTGGCGGTTCAAAGCCgtg                                          <  1:223489/128‑1 (MQ=255)                                                                                         aaTGCGCGTAGCAAGGTCACAAAAGACCTGCATCAAAAAAAAAAAAACATAGATACCGAACATTATGGTCTGGAGCACGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGGcaacaac                                  >  1:158634/1‑128 (MQ=255)                                                                                             cgcgTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCTTTGAGTATCTTGCGGTTAAAAGCCGTGTCAACAAAATCa                              <  1:294170/129‑1 (MQ=255)                                                                                                 tAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAgga                          >  2:334901/1‑128 (MQ=255)                                                                                                 tAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCAACAAGCTCCGGAGCGAGCACGAAACAGGATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAggg                          >  2:442008/1‑129 (MQ=255)                                                                                                  aGCAAGGTCACAAAAGACCTGCGTCAGGAGCAGGAAATCCTTAATACAGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGa                         >  1:617413/1‑129 (MQ=255)                                                                                                     aaGGACAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGACTGGAACGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCCAACCCGTGTCAGCAAAATCAAGGGAcgg                      >  2:613904/1‑127 (MQ=255)                                                                                                          cagaaaaGACCTGCGGCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCct                 <  2:417265/126‑1 (MQ=255)                                                                                                                     tGCGTCAGGAGCAGGAAATCATTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCAGTGTCAACAAAATCAAGGGCCCGATACTGTGCCGGGTAg      >  2:460847/1‑129 (MQ=255)                                                                                                                      gCGTCAGGAGCAGGAACTCCTTGATACCGACCATTATGAACTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAgg     >  1:35050/1‑129 (MQ=255)                                                                                                                        gTCAGGAGGAGGAAATCGTGGATACCGACCATTATGGTCGGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTGGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGgc   <  2:492550/129‑1 (MQ=255)                                                                                                                         tCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGAGCGGCTGCCTGGTAGGgcc  >  2:347854/1‑129 (MQ=255)                                                                                                                              |                                                                                                                            GTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCC  >  NC_000913/459709‑459957

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑