Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	401,440	T→G	S19A (TCG→GCG)	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,440	0	T	G	100.0%	211.5 / NA	60	S19A (TCG→GCG)	iraP	anti‑RssB factor, RpoS stabilzer during Pi starvation; anti‑adapter protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/39);  total (21/39)

ACATAATGTTAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTG  >  NC_000913/401342‑401533                                                                                                   |                                                                                              aCATAATGTTAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCGCATTGCTGAGTTGTTATTTAAGCTTGCCC‑AAAAGAAGAAGAGGCg                                                                                             <  2:1121639/100‑1 (MQ=255)      aTGTTAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGACGAAGAGGCGAAAg                                                                                         <  2:804969/100‑1 (MQ=255)       tGTTAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGa                                                                                        <  2:14035/100‑1 (MQ=255)        gTTAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGaa                                                                                       >  1:731288/1‑100 (MQ=255)         ttAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAAc                                                                                      >  2:172819/1‑100 (MQ=255)         ttAACTTCTCCATACTTAGGATAAGGAAATACAGAAATGAAAAATATCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGaan                                                                                      <  2:314866/100‑2 (MQ=255)           aaCTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACtg                                                                                    <  2:220309/100‑1 (MQ=255)           aaCTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACtg                                                                                    <  1:332999/100‑1 (MQ=255)           aaCTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACtg                                                                                    <  2:1214945/100‑1 (MQ=255)                   cATACTTTGGATAAGGAAATACAGACATGAATAATCTTATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCa                                                                            <  2:881127/100‑1 (MQ=255)                     tACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAgg                                                                          <  1:365249/100‑1 (MQ=255)                     tACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAg                                                                           >  2:575304/1‑99 (MQ=255)                     ccgatctgATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAgg                                                                          <  2:375620/93‑1 (MQ=255)                      aCTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGt                                                                         >  2:370493/1‑100 (MQ=255)                       cTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTa                                                                        <  1:19203/100‑1 (MQ=255)                       cTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTa                                                                        >  2:147882/1‑100 (MQ=255)                       cTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTa                                                                        <  2:1007604/100‑1 (MQ=255)                          tGGATAAGGAACGACACACATCAACAATCTCATTGCTGAGTTGTTATTTCAGCTTGCACCCCACGAAGACGGGGCCAAACAACTGTGTGCGCAGGTAg                                                                       <  2:760245/98‑1 (MQ=255)                          tGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGaa                                                                     <  2:658640/100‑1 (MQ=255)                          tGGATAAGGAAATACAGACATGAAAAATCTCAATGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGaa                                                                     >  1:1236549/1‑100 (MQ=255)                           ggATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAg                                                                    <  2:767234/100‑1 (MQ=255)                           ggATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAg                                                                    <  1:670315/100‑1 (MQ=255)                            gATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAggagatcgg                                                                   >  1:375623/1‑93 (MQ=255)                              tAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCtt                                                                 >  2:179690/1‑100 (MQ=255)                               aaGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCttt                                                                >  1:21143/1‑100 (MQ=255)                                aGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTg                                                               <  1:588362/100‑1 (MQ=255)                                 ggAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTgg                                                              <  2:376609/100‑1 (MQ=255)                                 ggAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTgg                                                              >  1:1126826/1‑100 (MQ=255)                                  gAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTgg                                                              >  2:789188/1‑99 (MQ=255)                                  gAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGga                                                             <  1:1116299/100‑1 (MQ=255)                                   aaaTACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGgag                                                            <  2:731283/100‑1 (MQ=255)                                    aaTACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGgaga                                                           <  1:536793/100‑1 (MQ=255)                                       aCAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTa                                                        <  2:21143/100‑1 (MQ=255)                                           aCATTAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGTGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGt                                                    <  2:608678/100‑1 (MQ=255)                                                 aaaaTCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGc                                              <  2:1126816/100‑1 (MQ=255)                                                 aaaaTCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGc                                              >  2:1122917/1‑100 (MQ=255)                                                      ctcATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGc                                         >  2:339497/1‑100 (MQ=255)                                                         aTTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTc                                      >  1:415640/1‑100 (MQ=255)                                                         aTTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTc                                      >  1:1149849/1‑100 (MQ=255)                                                            gCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCa                                   <  2:528068/100‑1 (MQ=255)                                                             cTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCaa                                  >  1:15286/1‑100 (MQ=255)                                                                aGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAatat                               <  2:154533/100‑1 (MQ=255)                                                                   tgttATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGgc                            >  2:752686/1‑100 (MQ=255)                                                                    gttATTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGgcg                           <  1:172821/100‑1 (MQ=255)                                                                       aTTTAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCaa                        >  2:526360/1‑100 (MQ=255)                                                                          tAAGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAAt                     <  1:179692/100‑1 (MQ=255)                                                                           aatcTTTCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAAGTAGAAGCTTTGGAGATTATCGTCACTGCATTGCTTCGCAATATGGCGCaaaa                      <  2:755494/95‑1 (MQ=255)                                                                           aaGCTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATg                    <  2:183772/100‑1 (MQ=255)                                                                              cTTGCCCAAAAAGAAGAAGAGGCGAAGGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGAcc                 <  1:373451/100‑1 (MQ=255)                                                                              cTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGAcc                 <  1:526364/100‑1 (MQ=255)                                                                              cTTGCCCAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGAcc                 <  1:370496/100‑1 (MQ=255)                                                                                   ccAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAg            >  1:546359/1‑100 (MQ=255)                                                                                    cAAAAAGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGc           <  1:109960/100‑1 (MQ=255)                                                                                      aaaaGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCgg         <  1:339500/100‑1 (MQ=255)                                                                                       aaaGAAGAAGAGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGt        <  2:97432/100‑1 (MQ=255)                                                                                         agaagaagaGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGttg      >  1:18487/1‑100 (MQ=255)                                                                                         agaagaagaGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGttg      <  2:619843/100‑1 (MQ=255)                                                                                             gaagaGGCGAAATAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGttgattg  >  1:855014/1‑100 (MQ=255)                                                                                             gaagaGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGttgattg  <  1:464224/100‑1 (MQ=255)                                                                                             gaagaGGCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGttgattg  <  1:1166101/100‑1 (MQ=255)                                                                                                   |                                                                                              ACATAATGTTAACTTCTCCATACTTTGGATAAGGAAATACAGACATGAAAAATCTCATTGCTGAGTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTG  >  NC_000913/401342‑401533

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