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_000913401,5440G.100.0% 171.3 / NA 41coding (159/261 nt)iraPanti‑RssB factor, RpoS stabilzer during Pi starvation, anti‑adapter protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/21);  total (20/21)

GTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAAGT  >  NC_000913/401406‑401678
                                                                                                                                          |                                                                                                                                      
gttgttATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTagag                                                                                                                                        <  1:344929/139‑1 (MQ=255)
     tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                   <  2:431894/139‑3 (MQ=255)
     tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                   <  1:355817/139‑3 (MQ=255)
       tttAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcgct                                                                                                                                 <  1:546365/139‑5 (MQ=255)
               tGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTCCGCAATATGGCCCAAAATGACCAACAGTGGTTGATTGATCAGGTAGAGGAGCGCTGTACGAgg                                                                                                                        >  1:254815/1‑139 (MQ=255)
                gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGt                                                                                                                       >  2:583356/1‑139 (MQ=255)
                 cgcAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTa                                                                                                                      >  1:158446/3‑139 (MQ=255)
                          aagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                             <  2:780792/139‑1 (MQ=255)
                           agaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                   >  1:350302/1‑115 (MQ=255)
                           agaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                   <  2:350302/117‑3 (MQ=255)
                            gaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGc                                                                                                           <  2:13726/139‑1 (MQ=255)
                                              tgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGAcgacga                                                                                         >  1:56936/1‑139 (MQ=255)
                                                tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                       >  2:489559/1‑139 (MQ=255)
                                                   gcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                                                                    >  2:226415/1‑139 (MQ=255)
                                                     aGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAg                                                                                  >  2:133116/1‑139 (MQ=255)
                                                        tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                                                               >  1:350607/1‑139 (MQ=255)
                                                        tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                                                               <  2:37966/139‑1 (MQ=255)
                                                                ttGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGa                                                                       <  1:133116/139‑1 (MQ=255)
                                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCa                                                                                                         >  2:192209/1‑104 (MQ=255)
                                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCa                                                                                                         <  1:192209/104‑1 (MQ=255)
                                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                               <  1:93267/114‑1 (MQ=255)
                                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                               >  2:93267/1‑114 (MQ=255)
                                                                  ggAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                     <  2:745817/139‑1 (MQ=255)
                                                                      aTTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGt                                                                 >  2:457096/1‑139 (MQ=255)
                                                                             tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                          >  1:320613/1‑139 (MQ=255)
                                                                                 tGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTAt                                                      >  2:647446/1‑139 (MQ=255)
                                                                                     aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                                  <  2:90198/139‑1 (MQ=255)
                                                                                     aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                                  <  1:94523/139‑1 (MQ=255)
                                                                                      tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAg                                                 <  2:320613/139‑1 (MQ=255)
                                                                                          tCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATc                                             <  2:231789/139‑1 (MQ=255)
                                                                                           cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                               >  1:680319/1‑120 (MQ=255)
                                                                                           cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                               <  2:680319/120‑1 (MQ=255)
                                                                                              aaTATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCg                                         <  2:56936/139‑1 (MQ=255)
                                                                                                    gcgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTa                                   <  1:312736/139‑1 (MQ=255)
                                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTa                                                       <  2:206188/118‑1 (MQ=255)
                                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTa                                                       >  1:206188/1‑118 (MQ=255)
                                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTa                                   <  1:756718/138‑1 (MQ=255)
                                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTa                                   >  2:756718/1‑138 (MQ=255)
                                                                                                      gcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGcc                                                                                                          <  2:506312/66‑1 (MQ=255)
                                                                                                      gcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGcc                                                                                                          >  1:506312/1‑66 (MQ=255)
                                                                                                          aaTGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCTAGCATTCCTGACGACGATACGGAGCTGCTgcg                                                                          <  2:729551/94‑1 (MQ=255)
                                                                                                          aaTGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTACTGACGACGATACGGAGCTGCTgcg                                                                          >  1:729551/1‑94 (MQ=255)
                                                                                                             gACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                   <  1:308146/35‑3 (MQ=255)
                                                                                                             gACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                                   >  2:308146/1‑33 (MQ=255)
                                                                                                                 aaCAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCtt                      <  1:583356/139‑1 (MQ=255)
                                                                                                                         ttgattgaTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg              <  1:489559/139‑1 (MQ=255)
                                                                                                                         ttgattgaTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg              <  2:339717/139‑1 (MQ=255)
                                                                                                                                   aGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATaaa    >  2:81896/1‑139 (MQ=255)
                                                                                                                                     gTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAAGt  >  2:793889/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GTTGTTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAAGT  >  NC_000913/401406‑401678

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