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% 85.3 / NA 21coding (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 . (12/9);  total (12/9)

TTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCAT  >  NC_000913/401410‑401673
                                                                                                                                      |                                                                                                                                 
ttATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggc                                                                                                                               <  2:119339/139‑2 (MQ=255)
 tATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAATAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              <  2:480617/139‑3 (MQ=255)
       aGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCGGGTAGAGGGGCGCTGTAc                                                                                                                       <  2:736761/139‑1 (MQ=255)
            gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGt                                                                                                                  >  2:225110/1‑139 (MQ=255)
             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTa                                                                                                                 >  2:235031/1‑139 (MQ=255)
                            agTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGc                                                                                                  >  2:457308/1‑139 (MQ=255)
                                   aGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTg                                                                                           >  1:757107/1‑139 (MQ=255)
                                               gcAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                                                               >  2:247004/1‑139 (MQ=255)
                                                    tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                                                          <  1:225110/139‑1 (MQ=255)
                                                              ggAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  1:247004/139‑1 (MQ=255)
                                                                     aTCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         <  1:457308/139‑1 (MQ=255)
                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                       <  1:353052/139‑1 (MQ=255)
                                                                            cTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              >  1:414565/1‑62 (MQ=255)
                                                                            cTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              <  2:414565/64‑3 (MQ=255)
                                                                               cAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTg                                               <  2:70152/139‑1 (MQ=255)
                                                                                  tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              >  2:127861/1‑56 (MQ=255)
                                                                                  tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              <  1:127861/58‑3 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa                                          >  2:701808/1‑139 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGGGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATcc                                       >  2:166296/1‑139 (MQ=255)
                                                                                        gCAATATGGCTCAAAATGACCATCCGCGGTTGATTGATCAGGTAGAGGTGCGCTGTACGAGGTAAAGCCCGACGCCAGTATTCCTTACGACGATACGGGGCTGCTGTGCGAATACGTAATGAAGTTATTGaattctgct                                      <  1:10136/139‑8 (MQ=255)
                                                                                            tatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTc                                  >  1:128861/1‑139 (MQ=255)
                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAg                                >  1:535337/1‑139 (MQ=255)
                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAg                                >  2:65566/1‑139 (MQ=255)
                                                                                                             aaCAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCtt                 <  1:10814/139‑1 (MQ=255)
                                                                                                                     ttgattgaTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg         >  1:710570/1‑139 (MQ=255)
                                                                                                                     ttgattgaTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg         <  2:530171/139‑1 (MQ=255)
                                                                                                                            atCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCAt  >  2:475739/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                 
TTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCAT  >  NC_000913/401410‑401673

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