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% 232.3 / NA 57coding (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 . (31/26);  total (31/26)

GCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAA  >  NC_000913/401422‑401676
                                                                                                                          |                                                                                                                                    
gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggc                                                                                                                                  <  2:202328/127‑2 (MQ=255)
gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggc                                                                                                                                  >  1:202328/1‑126 (MQ=255)
gCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGt                                                                                                                     >  2:1354374/1‑139 (MQ=255)
 cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACg                                                                                                                         >  2:1370511/1‑134 (MQ=255)
 cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACg                                                                                                                         <  1:1370511/134‑1 (MQ=255)
          aagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                           <  2:1645499/139‑1 (MQ=255)
                 gTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCa                                                                                                    <  2:1258409/139‑1 (MQ=255)
                  tCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCAt                                                                                                   >  1:261417/1‑139 (MQ=255)
                           cTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          <  2:519402/139‑1 (MQ=255)
                            tgtgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          <  1:1164162/138‑1 (MQ=255)
                            tgtgtgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          >  2:1164162/1‑138 (MQ=255)
                                tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                     <  1:1354374/139‑1 (MQ=255)
                                tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                     <  1:1267761/139‑1 (MQ=255)
                                        tAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctg                                                                             >  1:744694/1‑139 (MQ=255)
                                            aGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCtgctgc                                                                         >  1:773933/1‑139 (MQ=255)
                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                               <  1:15630/112‑1 (MQ=255)
                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                               >  2:15630/1‑112 (MQ=255)
                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                    >  1:406306/1‑139 (MQ=255)
                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                    >  1:515210/1‑139 (MQ=255)
                                                 tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCACGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                    >  2:1097798/1‑139 (MQ=255)
                                                  ggAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                   <  2:261417/139‑1 (MQ=255)
                                                           cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          <  2:1290463/107‑1 (MQ=255)
                                                           cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          >  1:1290463/1‑107 (MQ=255)
                                                             tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                        >  1:387079/1‑139 (MQ=255)
                                                                cTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAgggg                                                                                                                                   <  1:36681/62‑1 (MQ=255)
                                                                cTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAgggg                                                                                                                                   >  2:36681/1‑62 (MQ=255)
                                                                cTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTa                                                     >  2:509061/1‑139 (MQ=255)
                                                                 tGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGTTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCa                                                                                                       <  1:1029214/88‑1 (MQ=255)
                                                                 tGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCa                                                                                                       >  2:1029214/1‑88 (MQ=255)
                                                                     aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                             <  1:629984/94‑1 (MQ=255)
                                                                     aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                             >  2:629984/1‑94 (MQ=255)
                                                                           cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATcc                                          >  1:180130/1‑139 (MQ=255)
                                                                                tatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctgct                                                                           >  1:740696/1‑101 (MQ=255)
                                                                                tatGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctgct                                                                           <  2:740696/101‑1 (MQ=255)
                                                                                  tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGGTGCCAGCATTCCTGACGACGATACGGAGCTGCTGTGCGATTACGTAAAGAAGTTATTGaa                                                >  2:1156425/1‑126 (MQ=255)
                                                                                  tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGTGCGATTACGTAAAGAAGTTATTGaa                                                <  1:1156425/126‑1 (MQ=255)
                                                                                   ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACg                                                                                                                         >  1:986945/1‑52 (MQ=255)
                                                                                   ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACg                                                                                                                         <  2:986945/52‑1 (MQ=255)
                                                                                   ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGt                                                                                                                     <  2:996526/56‑1 (MQ=255)
                                                                                   ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGt                                                                                                                     >  1:996526/1‑56 (MQ=255)
                                                                                    gcgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGCAGTTATTGAAGctgtctcttataca                                 >  1:1091099/1‑126 (MQ=255)
                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                           >  2:1273278/1‑64 (MQ=255)
                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGAt                                                                                                           <  1:1273278/64‑1 (MQ=255)
                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                                                      >  2:159910/1‑101 (MQ=255)
                                                                                     cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                                                      <  1:159910/101‑1 (MQ=255)
                                                                                              aCCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt                       <  2:180130/139‑1 (MQ=255)
                                                                                              aCCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt                       <  2:935949/139‑1 (MQ=255)
                                                                                                    aGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          <  2:1078141/66‑1 (MQ=255)
                                                                                                    aGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacga                                                                                          >  1:1078141/1‑66 (MQ=255)
                                                                                                      cGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                          <  2:933431/96‑1 (MQ=255)
                                                                                                      cGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaaga                                                          >  1:933431/1‑96 (MQ=255)
                                                                                                        gTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGt                                      >  2:412806/1‑114 (MQ=255)
                                                                                                        gTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGt                                      <  1:412806/114‑1 (MQ=255)
                                                                                                        gTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACa             >  2:323887/1‑139 (MQ=255)
                                                                                                         ttgattgaTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg            <  1:509061/139‑1 (MQ=255)
                                                                                                         ttgattgaTCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAg            >  2:1576658/1‑139 (MQ=255)
                                                                                                                   aGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGc           <  2:548440/130‑1 (MQ=255)
                                                                                                                   aGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGc           >  1:548440/1‑130 (MQ=255)
                                                                                                                   aGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACATAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATaaa  >  2:1618777/1‑139 (MQ=255)
                                                                                                                     gTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                               <  1:217547/44‑1 (MQ=255)
                                                                                                                     gTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCt                                                                                               >  2:217547/1‑44 (MQ=255)
                                                                                                                          |                                                                                                                                    
GCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCATAAA  >  NC_000913/401422‑401676

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