Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,734 T→A 16.8% K59N (AAA→AAT) ‡ PP_4851 ← phosphate starvation‑inducible protein PsiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,518,7340TA16.8% 102.7 / 16.9 54K59N (AAA→AAT) ‡PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base A (5/4);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

CAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAG  >  NC_002947/5518586‑5518881
                                                                                                                                                     |                                                                                                                                                   
cAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCt                                                                                                                                                     <  5:9942/149‑1 (MQ=255)
          cttcatcAGGCAATCGCTCATGAAGGCCTTGAGCTCATCCACCTTCAGGGGGTTGGTGGTGGCGTCAGCATTGCAGGTCTTAATTTTTTCATGCTGGGGTTGGGGCGCATCAATTAGGAGGCAGGTGCACATGCCCGAGATGCGTTCaa                                                                                                                                           >  8:113879/1‑148 (MQ=255)
           ttcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATc                                                                                                                                          >  4:110504/1‑149 (MQ=255)
             cttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAtga                                                                                                                                         <  2:164324/148‑3 (MQ=255)
             cttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAtg                                                                                                                                          <  7:237028/147‑2 (MQ=255)
                 aGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCtt                                                                                                                                    <  3:163214/149‑1 (MQ=255)
                   gCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTc                                                                                                                                             <  2:65478/138‑1 (MQ=255)
                   gCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTc                                                                                                                                             >  1:65478/1‑138 (MQ=255)
                          cTCATGAAGGCCTTGCGCTCATCCGCCTTCAGGGATTTGGTGGTGGCGTCAGCATTTCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCtt                                                                                                                           <  3:197794/149‑1 (MQ=255)
                          cTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTGTTGCTGGGTTTGGGGGACATCCTGTTTCATGACGGTGCTCATGATCGCTTTTCGTTCATCGCCCTTGAGGGccgg                                                                                                                           >  5:152957/1‑147 (MQ=255)
                             aTGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTg                                                                                                                          <  3:110504/147‑1 (MQ=255)
                                aaGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGt                                                                                                                     <  3:168524/149‑1 (MQ=255)
                                    ccTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAgg                                                                                                                                >  4:40788/1‑134 (MQ=255)
                                    ccTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAgg                                                                                                                                <  3:40788/134‑1 (MQ=255)
                                     cTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGc                                                                                                                  <  8:213605/147‑1 (MQ=255)
                                     cTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGc                                                                                                                  >  7:213605/1‑147 (MQ=255)
                                                  cccTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGGACATCCTTTTTCAGGAAGGTGCTCATGAACGCTTTGCGTTCATCGACCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGt                                                                                                   >  1:155790/1‑149 (MQ=255)
                                                    cTTCAGGGCTTTGGTTGTTGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTCGGGCACATCGTTTTTCAGGCTGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGttt                                                                                                 <  3:196501/149‑1 (MQ=255)
                                                                     tggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCtgctgc                                                                                <  1:115079/149‑1 (MQ=255)
                                                                       gCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTCAGGGCCTTGGCGGTGGCGTCGGCCTTGCAGGTTTTCATCTTTTCCTCCTgcgc                                                                              >  7:250050/1‑149 (MQ=255)
                                                                         gTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCgg                                                                            >  8:80514/1‑149 (MQ=255)
                                                                              cATTGCAGGTCTTCATATTTTCTTGCTGGGTTTCGGGCACATCCTTTTTCAGGCAGGTCCTCATGACCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGCTTTTCATCTTTTCCTGCTGCGCGGTATTg                                                                       >  2:53513/1‑149 (MQ=255)
                                                                                  gCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTAATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAg                                                                   >  5:73858/1‑149 (MQ=255)
                                                                                      gTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGa                                                                >  8:89647/1‑148 (MQ=255)
                                                                                         ttCATTTTTTCTTGCTTTGTTTGGTGCACATCCTTTTTCAGGCATGTGCTCATGAACGCTTTGCTTTCATCGCCCTTGAGGGCCTTGGCGGTGTCGTCTGCGTTGCAGTTTTTCATCTTTTCctgct                                                                                  <  5:176950/127‑1 (MQ=255)
                                                                                         ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGccc                                                                                                                                      <  7:179037/75‑1 (MQ=255)
                                                                                         ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGccc                                                                                                                                      >  8:179037/1‑75 (MQ=255)
                                                                                         ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAcc                                                            >  5:78073/1‑149 (MQ=255)
                                                                                         ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAAAGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTAatgct                                                                                  >  6:176950/1‑127 (MQ=255)
                                                                                            aTTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTg                                                         >  6:174595/1‑149 (MQ=255)
                                                                                            aTTTTTTCATGCTGGGTTTCGGGCACATCCTTTTTCAGGCACGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCtt                                                          <  5:51940/148‑1 (MQ=255)
                                                                                             ttttttCTTGCTGGGTTTGGTGCACATTCTTTTTCAGGCATGTGCTCATGATCGGTTTTCGTTCATAGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTg                                                         <  2:155790/148‑1 (MQ=255)
                                                                                             ttttttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                        <  6:78073/149‑1 (MQ=255)
                                                                                             ttttttCTTGCTGGGGGTGGGGCACATCCTTTTTCAGGCAGGGGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                        <  3:175932/149‑1 (MQ=255)
                                                                                               ttttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTggg                                                       <  2:130137/148‑1 (MQ=255)
                                                                                                      gCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCGTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGa                                                 >  1:59761/1‑148 (MQ=255)
                                                                                                      gCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTg                                                  <  5:174595/146‑1 (MQ=255)
                                                                                                          ggTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                        >  5:98632/1‑136 (MQ=255)
                                                                                                          ggTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                        <  6:98632/136‑1 (MQ=255)
                                                                                                                      aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgcca                               <  2:59761/149‑1 (MQ=255)
                                                                                                                          tttttCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAACTTGGGCACTGAACAACACCGCCAGCGCCAGTaa                           >  8:39415/1‑149 (MQ=255)
                                                                                                                             ttCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAAct                         <  7:80514/148‑2 (MQ=255)
                                                                                                                                aggcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACAcc                                         >  2:109370/1‑129 (MQ=255)
                                                                                                                                aggcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACAcc                                         <  1:109370/129‑1 (MQ=255)
                                                                                                                                    aggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCtt                                                          <  5:93453/108‑1 (MQ=255)
                                                                                                                                    aggTGCTCATGAAAGCTTTGCGTTCATCGCACTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCtt                                                          >  6:93453/1‑108 (MQ=255)
                                                                                                                                          cgaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCAc                                                    >  2:74755/3‑108 (MQ=255)
                                                                                                                                          cgaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCAc                                                    <  1:74755/106‑1 (MQ=255)
                                                                                                                                           gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt          >  5:8691/2‑149 (MQ=255)
                                                                                                                                           gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt          >  5:4865/2‑149 (MQ=255)
                                                                                                                                           gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt          >  3:75773/2‑149 (MQ=255)
                                                                                                                                              gAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGTGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCa       >  5:108749/1‑149 (MQ=255)
                                                                                                                                              gAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCa       >  7:150433/1‑149 (MQ=255)
                                                                                                                                               aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCAt      <  3:223104/149‑1 (MQ=255)
                                                                                                                                                   cTTTGCGTTCATCGCCCTTGAGGGCCTTGGTGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAg  <  6:108749/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAG  >  NC_002947/5518586‑5518881

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: