Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,433,765	G→A	A26T (GCA→ACA)	def →	peptide deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,433,765	0	G	A	100.0%	215.8 / NA	66	A26T (GCA→ACA)	def	peptide deformylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/35);  total (31/35)

AGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGC  >  NC_000913/3433645‑3433900                                                                                                                         |                                                                                                                                        agaagaAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTAtcg                                                                                                                       <  1:889016/139‑1 (MQ=255)        ataatCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTa                                                                                                                          <  1:563502/129‑1 (MQ=255)        ataatCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTa                                                                                                                          >  2:563502/1‑129 (MQ=255)        ataatCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCgatga                                                                                                                <  2:87964/139‑1 (MQ=255)        ataatCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCgatga                                                                                                                >  1:87964/1‑139 (MQ=255)              tttCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGtt                                                                                                          <  1:880331/139‑1 (MQ=255)                            acacATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACg                                                                                            <  2:572704/139‑1 (MQ=255)                                aTCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCaga                                                                                        >  1:1194077/1‑139 (MQ=255)                                aTCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCaga                                                                                        >  2:56165/1‑139 (MQ=255)                                  ctctGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCagaag                                                                                      >  1:103775/1‑139 (MQ=255)                                  ctctGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCagaag                                                                                      >  1:481152/1‑139 (MQ=255)                                   tctGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCagaaga                                                                                     <  2:1194077/139‑1 (MQ=255)                                   tctGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCagaaga                                                                                     <  2:338567/139‑1 (MQ=255)                                          aTTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTAtt                                                                              >  2:1503970/1‑139 (MQ=255)                                          aTTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTAtt                                                                              >  2:193162/1‑139 (MQ=255)                                               tGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGAt                                                                                                 >  1:261593/1‑115 (MQ=255)                                               tGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGAt                                                                                                 <  2:261593/115‑1 (MQ=255)                                               tGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCt                                                                         >  1:497323/1‑139 (MQ=255)                                                    gTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCa                                                                                          <  2:756407/117‑1 (MQ=255)                                                    gTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCa                                                                                          >  1:756407/1‑117 (MQ=255)                                                        tGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAAc                                                                <  2:103775/139‑1 (MQ=255)                                                             gtgtTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAgg                                                           <  1:274717/139‑1 (MQ=255)                                                                  aCATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGat                                                      <  2:1358042/139‑1 (MQ=255)                                                                        tCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcgg                                                                    >  1:1412869/1‑119 (MQ=255)                                                                        tCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcgg                                                                    <  2:1412869/119‑1 (MQ=255)                                                                        tCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCAt                                                <  2:497323/139‑1 (MQ=255)                                                                         ccGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATc                                               <  1:56165/139‑1 (MQ=255)                                                                          cGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTAt                                                                                                                         >  2:320580/1‑64 (MQ=255)                                                                          cGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTAt                                                                                                                         <  1:320580/64‑1 (MQ=255)                                                                          cGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcg                                                                     >  1:722274/1‑116 (MQ=255)                                                                          cGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcg                                                                     <  2:722274/116‑1 (MQ=255)                                                                             aCGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAAtt                                                                                                                                 <  1:1074659/53‑1 (MQ=255)                                                                             aCGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAAtt                                                                                                                                 >  2:1074659/1‑53 (MQ=255)                                                                             aCGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCagaaga                                                                                     <  1:299784/97‑1 (MQ=255)                                                                             aCGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCagaaga                                                                                     >  2:299784/1‑97 (MQ=255)                                                                              cGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGt                                          <  1:1503970/139‑1 (MQ=255)                                                                                aGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAAcc                                                               >  1:1246383/1‑116 (MQ=255)                                                                                aGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAAcc                                                               <  2:1246383/116‑1 (MQ=255)                                                                                 gCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATc                                       <  1:343234/139‑1 (MQ=255)                                                                                 gCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATc                                       >  2:513087/1‑139 (MQ=255)                                                                                  cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTa                                                                                                                          <  2:1024630/55‑1 (MQ=255)                                                                                  cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTa                                                                                                                          >  1:1024630/1‑55 (MQ=255)                                                                                  cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGc                                                                                           >  2:530030/1‑86 (MQ=255)                                                                                  cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGc                                                                                           <  1:530030/86‑1 (MQ=255)                                                                                      ttCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTg                                                        >  2:1137164/1‑117 (MQ=255)                                                                                      ttCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTg                                                        <  1:1137164/117‑1 (MQ=255)                                                                                       tCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAAt                                                                                                                                  <  2:658813/42‑1 (MQ=255)                                                                                       tCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATACAGAAAt                                                                                                                                  >  1:658813/1‑42 (MQ=255)                                                                                                     aaCCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATg                          <  1:1090318/132‑1 (MQ=255)                                                                                                     aaCCGGTAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATg                          >  2:1090318/1‑132 (MQ=255)                                                                                                           tAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGa                            <  1:1425449/124‑1 (MQ=255)                                                                                                           tAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGa                            >  2:1425449/1‑124 (MQ=255)                                                                                                           tAGAAGAAGTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCg             >  2:1082230/1‑139 (MQ=255)                                                                                                                 aagTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGa       >  1:907221/1‑139 (MQ=255)                                                                                                                 aagTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGa       <  1:513087/139‑1 (MQ=255)                                                                                                                 aagTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGa       <  1:193162/139‑1 (MQ=255)                                                                                                                  agTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTa                                                                                <  2:835195/65‑1 (MQ=255)                                                                                                                  agTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTa                                                                                >  1:835195/1‑65 (MQ=255)                                                                                                                  agTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcgg                                                                    >  1:283278/1‑77 (MQ=255)                                                                                                                  agTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcgg                                                                    <  2:283278/77‑1 (MQ=255)                                                                                                                  agTGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGaa      >  2:1179719/1‑139 (MQ=255)                                                                                                                    tGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCa                                      >  1:1434961/1‑105 (MQ=255)                                                                                                                    tGAATACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCa                                      <  2:1434961/105‑1 (MQ=255)                                                                                                                      aaTACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGata                                                     >  2:180936/1‑88 (MQ=255)                                                                                                                      aaTACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGata                                                     <  1:180936/88‑1 (MQ=255)                                                                                                                      aaTACAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGAAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTGCGGAAAACCGTGACGAACGGc  <  2:907221/139‑1 (MQ=255)                                                                                                                         |                                                                                                                                        AGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGC  >  NC_000913/3433645‑3433900

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