Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	461,007	T→G	I707S (ATC→AGC)	lon →	Lon protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	461,007	0	T	G	100.0%	172.1 / NA	55	I707S (ATC→AGC)	lon	Lon protease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (32/23);  total (32/23)

CATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCC  >  NC_000913/460876‑461146                                                                                                                                    |                                                                                                                                            cATCCACGTCCACGTCCCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTg                                                                                                                                  <  1:540833/143‑1 (MQ=255) cATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTg                                                                                                                                  <  1:938140/143‑1 (MQ=255) cATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTg                                                                                                                                  <  1:155811/143‑1 (MQ=255)  atacacgtccacgtACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGTCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGAGCACTCTGCGTgg                                                                                                                                 >  1:164696/4‑143 (MQ=255)  aTCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTgg                                                                                                                                 >  1:787415/1‑143 (MQ=255)     cacgtccacgtACCGGGAGGTGCGACGCCGAAGGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTCCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCa                                                                                                                              <  2:929792/143‑1 (MQ=255)     cacgtccacgtACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGGGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCa                                                                                                                              >  2:500524/1‑143 (MQ=255)      acgtccacgtACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAg                                                                                                                             >  1:295236/1‑143 (MQ=255)         tccacgtACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGCCAGGTa                                                                                                                          >  2:953040/1‑143 (MQ=255)                gcgggAGGGTGTGTCGCCGAAAGTTGGTGCGGGGGCCGGTATTGGTGTGTGCACCGCGCTGGTTTATTGTCTCGCCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGa                                                                                                                   <  2:65298/140‑1 (MQ=255)                            ttacccGAAAGCTGGTCCGAGTTCCGGTATGGCTATGTGCACCGCGCTGCTTTCTTGCCTGAACGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTga                                                                                                       <  1:499903/139‑1 (MQ=255)                             aCGCCGATCGATGGTCCGCGTGGCGGTGCTGCGATGTGCACCGCGCTGGTTGCTTGCCTGACCGGTAAGTCGGTTGGCGCCGATGTGGCACTGACCGGGGAGAGCTCCCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaa                                                                                                      <  2:717046/143‑1 (MQ=255)                             aCGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaa                                                                                                      <  1:708476/143‑1 (MQ=255)                                 cgcaaGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGCTCGGTGGTCTgaaagac                                                                                                  >  2:763111/4‑142 (MQ=255)                                 cGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCCGGTACTGCCGATCGGTGGTTTgaaagaa                                                                                                  >  1:673597/1‑143 (MQ=255)                                 cGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGGTTgaaagac                                                                                                  >  2:668848/1‑142 (MQ=255)                                 cGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGCCAGGTAGCGCCGATCGGTGGTTGgaaagac                                                                                                  >  2:996067/1‑142 (MQ=255)                                   aaaGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTTGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGaaaa                                                                                                <  1:837705/143‑1 (MQ=255)                                       aTGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGGTTGAAAGAAAAAcgc                                                                                            >  1:284469/1‑141 (MQ=255)                                         ggTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGGTTGAAAGAAAAACTCCt                                                                                          >  1:751527/1‑143 (MQ=255)                                          gTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATATGGCAATGACCGGTGAGAGCACTCTGCGTGGACAGGTACTGGCGATCGGTGGTCTCAAAGAAAAACTCCTg                                                                                         >  2:483011/1‑143 (MQ=255)                                           tCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGCGCACTCTGCGCGGTCAGGTACTGCCGATCGGTGGCCTGAAAGAAAAACTCCTgg                                                                                        >  2:241569/1‑143 (MQ=255)                                           tCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTgg                                                                                        >  2:885925/1‑143 (MQ=255)                                              gggcgCCGCTATTGCTATGTGCACCGCGCTGGTCTCGTGCCTGACCGGTAACCCGGTTCGTGCCGAGGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCAGTGGTTTGAAAGAAAAACTCCTGGCAg                                                                                     <  2:937564/139‑1 (MQ=255)                                              gAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaa                                                                                                     <  2:609008/127‑1 (MQ=255)                                              gAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaa                                                                                                     >  1:609008/1‑127 (MQ=255)                                               aGTGCCGGTATTGCTATGGGCACCGCGCGGGTTGCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGGGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAgc                                                                                    <  1:445397/143‑1 (MQ=255)                                                   ccGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCAt                                                                                >  2:800783/1‑143 (MQ=255)                                                      gTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcgc                                                                             >  2:720327/1‑143 (MQ=255)                                                         ttgcggcgtgCACCGCGCTGGTCTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCgcggc                                                                          <  2:321567/136‑1 (MQ=255)                                                                 tgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTaa                                                                  >  1:780964/1‑143 (MQ=255)                                                                 tgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAAATCCTGGCAGCGCAGCGCGGCggggtcaa                                                                  >  1:900368/1‑138 (MQ=255)                                                                    aCCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAAc                                                               >  1:488454/1‑143 (MQ=255)                                                                        cgcTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACCCTACGAATTTAGTTTCTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTg                                                           <  1:374138/143‑1 (MQ=255)                                                                             gTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAt                                                      >  2:997185/1‑143 (MQ=255)                                                                              tttCTGGCGTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAtt                                                     <  2:491174/143‑1 (MQ=255)                                                                               ttCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCCATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCATGGCAGCGCATCGCGGCGGGCTTAAAACAGAGCTAATTc                                                    >  2:739352/1‑143 (MQ=255)                                                                               ttCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGCTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTc                                                    >  1:90905/1‑143 (MQ=255)                                                                                  ttGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATGGCGGCGGGCTTCACACGGTCCTAATTCCGt                                                 >  2:760635/1‑143 (MQ=255)                                                                                          ccGGGAACCCGGGTCGGGCCGATGTGGCGATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAAt                                         <  2:990391/143‑1 (MQ=255)                                                                                           cGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGGTTAAAACAGTGCTAATTCCGTTCGacaaga                                        >  2:557938/1‑141 (MQ=255)                                                                                                aCCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAAcgc                                   <  1:901068/143‑1 (MQ=255)                                                                                                 cccGGTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAAcgcg                                  >  2:827932/1‑143 (MQ=255)                                                                                                    ggTTCGTGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCGGTACGAAAATAAACGCGATc                               >  2:389618/1‑143 (MQ=255)                                                                                                          tGCCGATGTGGCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCAGCGCGGCGGGAATAAAACAGTACTAATTCCGTTCGACAATAAACGCGATCAGGAag                         >  2:127960/1‑143 (MQ=255)                                                                                                            ccGACGCGGCAATGACCGGTGCCAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACCCCTTGCAGCGCATCGCGGCGTGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAagag                       <  1:585050/143‑1 (MQ=255)                                                                                                                agtgGTAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGTCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTc                   <  1:935785/142‑1 (MQ=255)                                                                                                                  gggCAATGACAGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCt                 <  1:641368/142‑1 (MQ=255)                                                                                                                  gcgCACTGGCCGGGGAGAGGACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCt                 <  2:459108/141‑1 (MQ=255)                                                                                                                    gCAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGa               >  1:770040/1‑143 (MQ=255)                                                                                                                     cAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTg                                                                                                                       <  1:528641/38‑1 (MQ=255)                                                                                                                     cAATGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTg                                                                                                                       >  2:528641/1‑38 (MQ=255)                                                                                                                      aaTGACCGGTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGCAAATAAACACGATCTGGAAGAGATTCCTGACa             >  2:805571/1‑143 (MQ=255)                                                                                                                             ggTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAAt      <  1:545333/143‑1 (MQ=255)                                                                                                                              gTGAGAGCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAGTAAACGCGATCTGGAAGAGATTCCTGACAACGTAAtt     >  1:506832/1‑143 (MQ=255)                                                                                                                                 agagCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTAAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAAGGTacctgac  >  2:132146/1‑137 (MQ=255)                                                                                                                                    |                                                                                                                                            CATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCC  >  NC_000913/460876‑461146

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