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)
             cccAAAAAGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTa                                                                                                                 >  1:290458/1‑139 (MQ=255)
                  aaaGAAGAAGAGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGcc                                                                                                            <  1:441645/139‑1 (MQ=255)
                     gaagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGTAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGa                                                                                                         <  2:93586/139‑1 (MQ=255)
                     gaagaagaGTCGAAAGAACTGTGTGCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGa                                                                                                         <  2:290458/139‑1 (MQ=255)
                                            tgCGCAGGTAGAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  <  2:567329/139‑1 (MQ=255)
                                                      gAAGCTTTGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAgctgct                                                                        >  1:1496726/1‑139 (MQ=255)
                                                           tttGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGcgcg                                                                   <  2:1496726/139‑1 (MQ=255)
                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGTCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                 >  1:944231/1‑139 (MQ=255)
                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                 >  2:1169497/1‑139 (MQ=255)
                                                             tGGAGATTATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAt                                                                 >  2:1009224/1‑139 (MQ=255)
                                                              ggAGATTATCGTCACTGCAATGCTTCGCAATATGTCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  2:944231/139‑1 (MQ=255)
                                                                    tATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                          >  1:141625/1‑139 (MQ=255)
                                                                    tATCGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaa                                                          >  2:818895/1‑139 (MQ=255)
                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                          <  1:258313/104‑1 (MQ=255)
                                                                       cGTCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTga                                                                                          >  2:258313/1‑104 (MQ=255)
                                                                         tCACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                     <  1:818895/139‑1 (MQ=255)
                                                                          cACTGCAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGt                                                    >  1:235985/1‑139 (MQ=255)
                                                                               cAATGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTg                                               <  1:1169497/139‑1 (MQ=255)
                                                                                aaTGCTTCGCAATATGGTGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTCAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         >  1:956392/1‑128 (MQ=255)
                                                                                aaTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         <  2:956392/128‑1 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacg                                                                                        <  1:136149/96‑1 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTgacg                                                                                        >  2:136149/1‑96 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             <  2:235985/139‑1 (MQ=255)
                                                                                 aTGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             <  1:1009224/139‑1 (MQ=255)
                                                                                  tGCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAg                                            >  2:1190432/1‑139 (MQ=255)
                                                                                   gCTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGc                                           >  1:975903/1‑139 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                                                               <  1:1529588/102‑1 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACgg                                                                               >  2:1529588/1‑102 (MQ=255)
                                                                                    cTTCGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCa                                          <  2:141625/139‑1 (MQ=255)
                                                                                      aagagacagggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGTGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATc                                        <  2:398122/130‑1 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                     >  1:1237576/1‑125 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAaagaag                                                     <  2:1237576/125‑1 (MQ=255)
                                                                                       cGCAATATGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATcc                                       >  2:243483/1‑139 (MQ=255)
                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  <  2:786124/89‑1 (MQ=255)
                                                                                              tGGCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  >  1:786124/1‑89 (MQ=255)
                                                                                               ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGt                               <  1:243483/139‑1 (MQ=255)
                                                                                               ggCGCAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGt                               >  1:1353380/1‑139 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                >  1:57250/1‑104 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  2:57250/104‑1 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             >  1:510725/1‑123 (MQ=255)
                                                                                                 cgcAAAATGACCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGaa                                             <  2:510725/123‑1 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGAGGGGGGCTGTACGAGGTAAAGCACGATGACAGAATTCCTGAAGAGGATACGGAGCTGATGGGCGAATAGGTAAAGAAGATATTGAAGCATCCTAGTCAGTaactggctctt                    >  1:398122/1‑130 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         <  1:993203/102‑1 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTaaa                                                         >  2:993203/1‑102 (MQ=255)
                                                                                                          aCCAACAGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAAt                    <  2:761915/139‑1 (MQ=255)
                                                                                                                aGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCTCGATGCCa                                                                                                    <  2:1085712/53‑1 (MQ=255)
                                                                                                                aGCGGTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCa                                                                                                    >  1:1085712/1‑53 (MQ=255)
                                                                                                                    gTTGATTGATCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACa          >  2:293551/1‑139 (MQ=255)
                                                                                                                        attgatCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                >  2:1057565/1‑81 (MQ=255)
                                                                                                                        attgatCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGAtt                                                                <  1:1057565/81‑1 (MQ=255)
                                                                                                                            atCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTAAAGAAGTTATTGAAGCATCCTCGTCAGTAAAAAGTTAATCTTTTCAACAGCTGTCAt  >  2:702069/1‑139 (MQ=255)
                                                                                                                             tCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  >  2:1253717/1‑58 (MQ=255)
                                                                                                                             tCAGGTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATa                                                                                  <  1:1253717/58‑1 (MQ=255)
                                                                                                                                 gTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                                                           <  2:664597/77‑1 (MQ=255)
                                                                                                                                 gTAGAGGGGCGCTGTACGAGGTAAAGCCCGATGCCAGCATTCCTGACGACGATACGGAGCTGCTGCGCGATTACGTa                                                           >  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: