Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,922,067	T→C	K125R (AAA→AGA)	atpI ←	ATP synthase accessory factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,922,067	0	T	C	100.0%	200.1 / NA	56	K125R (AAA→AGA)	atpI	ATP synthase accessory factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/26);  total (30/26)

ATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATCGCCAGAACTTTGAA  >  NC_000913/3921930‑3922200                                                                                                                                          |                                                                                                                                      attgTCCAGAAGGTGGCTGGGGGGGTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGGTGCCTTTTACCCTctg                                                                                                                                     >  2:935223/1‑140 (MQ=255)   tgTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTctgtt                                                                                                                                   >  1:824455/1‑140 (MQ=255)    gTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTctgttg                                                                                                                                  <  2:373679/140‑1 (MQ=255)    gTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCCTTTACCCTctg                                                                                                                                     >  2:171955/1‑137 (MQ=255)      ccAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTgttgtt                                                                                                                                >  2:617993/1‑140 (MQ=255)                gCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGcc                                                                                                                      >  1:716824/1‑140 (MQ=255)                    gggggTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGC‑GCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGc                                                                                                                 <  2:177358/140‑1 (MQ=255)                      gggTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGcc                                                                                                                <  1:826223/140‑1 (MQ=255)                      gggTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGcc                                                                                                                <  1:688061/140‑1 (MQ=255)                       ggTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCa                                                                                                               <  2:806691/140‑1 (MQ=255)                               ggATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGa                                                                                                       <  2:621566/140‑1 (MQ=255)                                     accaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAaccacc                                                                                                 >  1:469999/1‑140 (MQ=255)                                     accaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAac                                                                                                     <  2:749904/136‑1 (MQ=255)                                     accaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAac                                                                                                     >  1:749904/1‑136 (MQ=255)                                      ccaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAaccacc                                                                                                 >  1:143665/1‑139 (MQ=255)                                      ccaGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAaccacc                                                                                                 <  2:143665/139‑1 (MQ=255)                                       caGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAg                                                                                               <  1:895032/140‑1 (MQ=255)                                         gCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCa                                                                                             >  2:684430/1‑140 (MQ=255)                                          cGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCAc                                                                                            <  1:617993/140‑1 (MQ=255)                                             gaATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCaa                                                                                         <  2:469999/140‑1 (MQ=255)                                               aTGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCt                                                                                                         <  2:431683/122‑1 (MQ=255)                                               aTGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCt                                                                                                         >  1:431683/1‑122 (MQ=255)                                                 gTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAcc                                                                                     <  1:186036/140‑1 (MQ=255)                                                    cGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAc                                                                                      <  2:566850/136‑1 (MQ=255)                                                    cGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAc                                                                                      >  1:566850/1‑136 (MQ=255)                                                        ggTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAc                                                                                      <  1:463656/132‑1 (MQ=255)                                                        ggTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAAc                                                                                      >  2:463656/1‑132 (MQ=255)                                                         gTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAAc                                                                             >  1:570922/1‑140 (MQ=255)                                                            cAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGAt                                                                          >  1:833621/1‑140 (MQ=255)                                                                tGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGc                                                                      <  1:684430/140‑1 (MQ=255)                                                                tGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGc                                                                      <  2:716824/140‑1 (MQ=255)                                                                    ggTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCa                                                                  >  2:22236/1‑140 (MQ=255)                                                                                ggtgTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTCTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCt                                                      >  2:968912/1‑140 (MQ=255)                                                                                 gtgtCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCtt                                                     >  2:249775/1‑140 (MQ=255)                                                                                   gtCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTa                                                   >  2:156097/1‑140 (MQ=255)                                                                                   gtCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTa                                                   >  1:270156/1‑140 (MQ=255)                                                                                    tCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATAcc                                                          >  1:870798/1‑132 (MQ=255)                                                                                    tCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATAcc                                                          <  2:870798/132‑1 (MQ=255)                                                                                           tAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGcc                                           >  1:238920/1‑140 (MQ=255)                                                                                              tCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAAc                                        >  2:103380/1‑140 (MQ=255)                                                                                              tCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAAc                                        >  1:559910/1‑140 (MQ=255)                                                                                                cTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGc                                      >  2:598716/1‑140 (MQ=255)                                                                                                 tGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGcc                                     <  1:598716/140‑1 (MQ=255)                                                                                                 tGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGcc                                     <  1:22236/140‑1 (MQ=255)                                                                                                  gcggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGcca                                    <  2:559910/140‑1 (MQ=255)                                                                                                   cggcgTCATATTTTCTGAAGCCATGATGCCTTTTCCCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccac                                   <  2:358672/140‑1 (MQ=255)                                                                                                   cggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccac                                   >  1:597057/1‑140 (MQ=255)                                                                                                   cggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccac                                   >  2:424490/1‑140 (MQ=255)                                                                                                   cggcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGGATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccac                                   >  1:48195/1‑140 (MQ=255)                                                                                                     gcgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccacca                                 <  2:570922/140‑1 (MQ=255)                                                                                                      cgTCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccaccac                                <  2:833621/140‑1 (MQ=255)                                                                                                        tCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccaccacca                              <  1:849073/140‑1 (MQ=255)                                                                                                        tCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccaccacca                              >  2:760162/1‑140 (MQ=255)                                                                                                        tCATATTTTCTGAAGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGccaccacca                              <  1:760162/140‑1 (MQ=255)                                                                                                                    aaGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAaca                        <  1:858844/134‑1 (MQ=255)                                                                                                                    aaGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAaca                        >  2:858844/1‑134 (MQ=255)                                                                                                                     aGCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACAt                 >  2:459286/1‑140 (MQ=255)                                                                                                                      gCCATGATGCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATc                >  1:394913/1‑140 (MQ=255)                                                                                                                          tgatgCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATCGCCa            >  2:7830/1‑140 (MQ=255)                                                                                                                          tgatgCCTTTTACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATCGCCa            <  2:394913/140‑1 (MQ=255)                                                                                                                                    tACCCTCTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATCGCCAGAACTTTGaa  >  1:255569/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                      ATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAAGGTTATTCAGGTGGTGTCCTATGTAATCCTGCGGCGTCATATTTTCTGAAGCCATGATGCCTTTTACCCTTTGTTGTTAATTACAGCCGGTGCCAGTATCTGAACCACCAGCACCAAAACCCACGTAACGATCAGCGGCAAGAATACCGCCTTTAAAACCGCCAACGCCACCACCAGTAACACCAACATCGCCAGAACTTTGAA  >  NC_000913/3921930‑3922200

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