Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,735 T→A 18.2% K59I (AAA→ATA) ‡ 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,7350TA18.2% 119.6 / 16.9 54K59I (AAA→ATA) ‡PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base A (6/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.00e+00

CTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAG  >  NC_002947/5518596‑5518881
                                                                                                                                            |                                                                                                                                                  
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                                                                                                                                             >  1:65478/1‑138 (MQ=255)
         gCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTc                                                                                                                                             <  2:65478/138‑1 (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                                                                                                                                <  3:40788/134‑1 (MQ=255)
                          ccTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAgg                                                                                                                                >  4:40788/1‑134 (MQ=255)
                           cTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGc                                                                                                                  >  7:213605/1‑147 (MQ=255)
                           cTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGc                                                                                                                  <  8:213605/147‑1 (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)
                                                                                                                                            |                                                                                                                                                  
CTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAG  >  NC_002947/5518596‑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: