Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	401,548	(G)5→4	coding (163/261 nt)	iraP →	anti‑RssB factor, RpoS stabilzer during Pi starvation, anti‑adapter protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	401,544	0	G	.	100.0%	231.1 / NA	55	coding (159/261 nt)	iraP	anti‑RssB factor, RpoS stabilzer during Pi starvation, anti‑adapter protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (28/27);  total (28/27)

TATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAA  >  NC_000913/401411‑401676                                                                                                                                      |                                                                                                                                     tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                 <  1:640754/139‑3 (MQ=255)       aGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                          >  1:5379/1‑139 (MQ=255)            gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGt                                                                                                                     >  1:819543/1‑139 (MQ=255)            gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGt                                                                                                                     >  1:754055/1‑139 (MQ=255)                  aaaGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGcc                                                                                                               >  2:56234/1‑139 (MQ=255)                       agaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                           <  2:349795/138‑1 (MQ=255)                       agaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                           >  1:349795/1‑138 (MQ=255)                            agTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGc                                                                                                     >  1:485043/1‑139 (MQ=255)                              tCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                 >  1:142315/1‑108 (MQ=255)                              tCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                 <  2:142315/110‑3 (MQ=255)                               cGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCAtt                                                                                                  <  2:388258/139‑1 (MQ=255)                                     aaCTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgac                                                                                            <  1:60424/139‑1 (MQ=255)                                       cTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAAGGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          <  1:905567/139‑1 (MQ=255)                                        tgtgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacgac                                                                                         >  1:113208/1‑139 (MQ=255)                                         gtgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGacgacg                                                                                        >  1:529499/1‑139 (MQ=255)                                            tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                     >  2:130038/1‑139 (MQ=255)                                            tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                     >  1:781595/1‑139 (MQ=255)                                              cgcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACg                                                                                   >  2:742622/1‑139 (MQ=255)                                              cgcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACg                                                                                   <  1:423041/139‑1 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggc                                                                                                                                  >  2:509246/1‑86 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggc                                                                                                                                  <  1:509246/87‑2 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAg                                                                                >  1:494080/1‑136 (MQ=255)                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAg                                                                                <  2:494080/136‑1 (MQ=255)                                                         gCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTgcg                                                                        <  2:819543/139‑1 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGccc                                                                                                              >  2:866988/1‑97 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGccc                                                                                                              <  1:866988/97‑1 (MQ=255)                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                    >  2:883113/1‑139 (MQ=255)                                                              ggAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                   <  1:897249/139‑1 (MQ=255)                                                                     aTCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                            <  1:742622/139‑1 (MQ=255)                                                                     aTCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                            <  2:529499/139‑1 (MQ=255)                                                                      tCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                 <  2:192757/70‑3 (MQ=255)                                                                      tCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                 >  1:192757/1‑68 (MQ=255)                                                                      tCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaag                                                           <  1:56234/139‑1 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGcc                                                                                                               >  1:484159/1‑86 (MQ=255)                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGcc                                                                                                               <  2:484159/86‑1 (MQ=255)                                                                         tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                        >  1:826963/1‑139 (MQ=255)                                                                                  tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                          >  2:157292/1‑64 (MQ=255)                                                                                  tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                          <  1:157292/64‑1 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                          >  1:226640/1‑62 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTAc                                                                                                                          <  2:226640/62‑1 (MQ=255)                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa                                             >  1:556850/1‑139 (MQ=255)                                                                                      tCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATc                                           <  2:781595/139‑1 (MQ=255)                                                                                           atatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGt                                      <  2:113208/139‑1 (MQ=255)                                                                                                gcgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTa                                 <  1:58906/139‑1 (MQ=255)                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                     >  1:273606/1‑86 (MQ=255)                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                     <  2:273606/86‑1 (MQ=255)                                                                                                      aaTGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGt                           >  2:214243/1‑139 (MQ=255)                                                                                                        tGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTa                         >  1:607790/1‑139 (MQ=255)                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt                       <  2:30411/139‑1 (MQ=255)                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt                       <  1:418183/139‑1 (MQ=255)                                                                                                           ccAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGt                                                                                                                     <  2:968776/44‑1 (MQ=255)                                                                                                           ccAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGt                                                                                                                     >  1:968776/1‑44 (MQ=255)                                                                                                               cAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCtttt                  >  2:623217/1‑139 (MQ=255)                                                                                                               cAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCtt                    >  2:116159/1‑137 (MQ=255)                                                                                                               cAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCtt                    <  1:116159/137‑1 (MQ=255)                                                                                                                    gTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACa             >  2:48661/1‑139 (MQ=255)                                                                                                                     ttgattgaTCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg            <  2:754055/139‑1 (MQ=255)                                                                                                                      tgattgatCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGc           <  1:315207/139‑1 (MQ=255)                                                                                                                       gattgatCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCt          <  1:678617/139‑1 (MQ=255)                                                                                                                        attgatCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTg         <  2:5379/139‑1 (MQ=255)                                                                                                                              cAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATaa   >  1:922756/1‑139 (MQ=255)                                                                                                                               aGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATaaa  >  1:830034/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     TATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAA  >  NC_000913/401411‑401676

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