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.98.2% 221.5 / ‑2.6 55coding (159/261 nt)iraPanti‑RssB factor, RpoS stabilzer during Pi starvation, anti‑adapter protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (27/27);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TTATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCAT  >  NC_000913/401410‑401673
                                                                                                                                      |                                                                                                                                 
ttATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTagag                                                                                                                                   >  2:814883/1‑135 (MQ=255)
ttATTTAAGCTTGCCCAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTagag                                                                                                                                   <  1:814883/135‑1 (MQ=255)
             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              >  2:762291/1‑125 (MQ=255)
             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAggggcg                                                                                                                              <  1:762291/127‑3 (MQ=255)
             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTa                                                                                                                 >  1:290458/1‑139 (MQ=255)
                  aaaGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGcc                                                                                                            <  1:441645/139‑1 (MQ=255)
                     gaagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGTAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGa                                                                                                         <  2:93586/139‑1 (MQ=255)
                     gaagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGa                                                                                                         <  2:290458/139‑1 (MQ=255)
                                            tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  <  2:567329/139‑1 (MQ=255)
                                                      gAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctgct                                                                        >  1:1496726/1‑139 (MQ=255)
                                                           tttGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                                                   <  2:1496726/139‑1 (MQ=255)
                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGTCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                 >  1:944231/1‑139 (MQ=255)
                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                 >  2:1169497/1‑139 (MQ=255)
                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                 >  2:1009224/1‑139 (MQ=255)
                                                              ggAGATTATCGTCACTGCAATGCTTCGCAATATGTCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  2:944231/139‑1 (MQ=255)
                                                                    tATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                          >  1:141625/1‑139 (MQ=255)
                                                                    tATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                          >  2:818895/1‑139 (MQ=255)
                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                          <  1:258313/104‑1 (MQ=255)
                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                          >  2:258313/1‑104 (MQ=255)
                                                                         tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                     <  1:818895/139‑1 (MQ=255)
                                                                          cACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGt                                                    >  1:235985/1‑139 (MQ=255)
                                                                               cAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTg                                               <  1:1169497/139‑1 (MQ=255)
                                                                                aaTGCTTCGCAATATGGTGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTCAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         >  1:956392/1‑128 (MQ=255)
                                                                                aaTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         <  2:956392/128‑1 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacg                                                                                        <  1:136149/96‑1 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacg                                                                                        >  2:136149/1‑96 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             <  2:235985/139‑1 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             <  1:1009224/139‑1 (MQ=255)
                                                                                  tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAg                                            >  2:1190432/1‑139 (MQ=255)
                                                                                   gCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGc                                           >  1:975903/1‑139 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                                                               <  1:1529588/102‑1 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                                                               >  2:1529588/1‑102 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa                                          <  2:141625/139‑1 (MQ=255)
                                                                                      aagagacagggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGTGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATc                                        <  2:398122/130‑1 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                     >  1:1237576/1‑125 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                     <  2:1237576/125‑1 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATcc                                       >  2:243483/1‑139 (MQ=255)
                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  <  2:786124/89‑1 (MQ=255)
                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  >  1:786124/1‑89 (MQ=255)
                                                                                               ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGt                               <  1:243483/139‑1 (MQ=255)
                                                                                               ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGt                               >  1:1353380/1‑139 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                >  1:57250/1‑104 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  2:57250/104‑1 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             >  1:510725/1‑123 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             <  2:510725/123‑1 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGAGGGGG‑GCTGTACGAGGTAAAGCACGATGACAGAATTCCTGAAGAGGATACGGAGCTGATGGGCGAATAGGTAAAGAAGATATTGAAGCATCCTAGTCAGTaactggctctt                    >  1:398122/1‑130 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         <  1:993203/102‑1 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         >  2:993203/1‑102 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt                    <  2:761915/139‑1 (MQ=255)
                                                                                                                aGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCTCGATGCCa                                                                                                    <  2:1085712/53‑1 (MQ=255)
                                                                                                                aGCGGTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCa                                                                                                    >  1:1085712/1‑53 (MQ=255)
                                                                                                                    gTTGATTGATCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACa          >  2:293551/1‑139 (MQ=255)
                                                                                                                        attgatCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                >  2:1057565/1‑81 (MQ=255)
                                                                                                                        attgatCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  1:1057565/81‑1 (MQ=255)
                                                                                                                            atCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCAt  >  2:702069/1‑139 (MQ=255)
                                                                                                                             tCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  >  2:1253717/1‑58 (MQ=255)
                                                                                                                             tCAGGTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  <  1:1253717/58‑1 (MQ=255)
                                                                                                                                 gTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                                                           <  2:664597/77‑1 (MQ=255)
                                                                                                                                 gTAGA‑GGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                                                           >  1:664597/1‑77 (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 < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑