Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,638,461 A→T 7.8% K393* (AAA→TAA)  PP_4103 → low‑affinity inorganic phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,638,4610AT 7.8% 146.1 / 2.7 51K393* (AAA→TAA) PP_4103low‑affinity inorganic phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base T (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TGGTACCGGCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGC  >  NC_002947/4638318‑4638608
                                                                                                                                               |                                                                                                                                                    
tGGTACCGGCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAac                                                                                                                                                >  1:149706/1‑149 (MQ=255)
  gTACCGGCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAac                                                                                                                                                <  6:307339/147‑1 (MQ=255)
  gTACCGGCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAac                                                                                                                                                <  7:33417/147‑1 (MQ=255)
  gTACCGGCGTTGGCGGCGCTGACCGGGATGAGCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAac                                                                                                                                                <  4:112999/147‑1 (MQ=255)
     ccGGCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAaccacc                                                                                                                                            <  4:106280/148‑1 (MQ=255)
       ggCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTgg                                                                                                                                         <  2:149706/149‑1 (MQ=255)
          gTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGa                                                                                                                                      <  8:172577/149‑1 (MQ=255)
              gcggcgCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGTCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACAtt                                                                                                                                  <  6:102993/149‑1 (MQ=255)
                ggcgCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTc                                                                                                                                 <  5:350497/148‑1 (MQ=255)
                ggcgCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTc                                                                                                                                 >  6:350497/1‑148 (MQ=255)
                       aCCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGa                                                                                                                         >  2:356097/1‑149 (MQ=255)
                         cGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGaca                                                                                                                       >  7:339592/1‑149 (MQ=255)
                         cGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGaca                                                                                                                       >  7:241790/1‑149 (MQ=255)
                                          agGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAg                                                                                                      <  8:272084/149‑1 (MQ=255)
                                           ggTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGATAAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCa                                                                                                       <  1:281156/148‑1 (MQ=255)
                                              cAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAgttgt                                                                                                  >  2:208799/1‑149 (MQ=255)
                                                       tGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGccc                                                                                         <  8:241790/149‑1 (MQ=255)
                                                       tGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGcc                                                                                          <  8:301176/148‑1 (MQ=255)
                                                             gctcAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCt                                                                                     <  1:208799/147‑1 (MQ=255)
                                                             gctcAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTg                                                                                    >  5:377008/1‑148 (MQ=255)
                                                             gctcAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGCGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTgg                                                                                   <  5:263182/149‑1 (MQ=255)
                                                                 aaGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGgacga                                                                               <  3:174434/149‑1 (MQ=255)
                                                                 aaGCGAGTACCGGCTGAAGCCATGGCCAACGTGAGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATTGAGA‑ACAAGCACCTGGTGACATTCGATGCCGACACCAGCAGACAAGTGCAGTTGTTCAAGGACAAGCTGgact                                                                                >  7:248846/1‑147 (MQ=255)
                                                                    cgcgTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGc                                                                            <  6:377008/149‑1 (MQ=255)
                                                                     gcgTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGcc                                                                           <  8:245719/149‑1 (MQ=255)
                                                                             gCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGT‑TAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCGAAGTGCAGTTGTTCTAGGCCAAGCTGTACGACGCCACCCTCta                                                                   >  6:297064/1‑149 (MQ=255)
                                                                             gCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACAAGTGAAGCAATCCGCCTCATGGAGT‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGACAAGTGCAGTTGTTCAAGGCCCAGATGGACGACGCCACCCGCta                                                                   >  6:291891/1‑149 (MQ=255)
                                                                              cTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGGGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCtat                                                                  <  3:54850/149‑1 (MQ=255)
                                                                                    ccATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCg                                                            <  8:339592/149‑1 (MQ=255)
                                                                                        ggCCAACTTGCGCAACGACATGTATCTGACCAGTTAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTTTTCAAGGCCCATCTGGACGACGCCAACCGATATATTACGTtgt                                                        >  6:362201/1‑149 (MQ=255)
                                                                                         gCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTtgtg                                                       <  2:168292/149‑1 (MQ=255)
                                                                                         cccAACGTGCGCATCGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCAGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTtgtg                                                       >  3:336135/2‑149 (MQ=255)
                                                                                                    cAACGACATGTATCTGACCAGTCAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTgg                                            <  6:151093/149‑1 (MQ=255)
                                                                                                      aCGACATGTATCTGACAATTGAAGCCATCCGCCTCATGGAGA‑AAAACCACATGGTTACAATCGATGCCGACACCAGAATACAAGTGCATTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTTGGTGAAGGTGGCt                                          <  8:248846/149‑1 (MQ=255)
                                                                                                          cATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCAGAAGCCAAGTGCAGTTGTTCAAGGCACAGCTGGACGACGCCACCCGCTATATTACGTTGTTGGTGAAGGTGGCTGTAg                                      >  7:264637/1‑149 (MQ=255)
                                                                                                           aTGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCATGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGaa                                                 >  1:122904/1‑137 (MQ=255)
                                                                                                           aTGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGaa                                                 <  2:122904/137‑1 (MQ=255)
                                                                                                                        gTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAg                                                <  3:6579/125‑1 (MQ=255)
                                                                                                                        gTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAg                                                >  4:6579/1‑125 (MQ=255)
                                                                                                                        gTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCgctgggg                        >  1:380199/1‑149 (MQ=255)
                                                                                                                         tGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCgctggggc                       <  1:356097/149‑1 (MQ=255)
                                                                                                                           aaGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGt                                                    >  8:72494/1‑118 (MQ=255)
                                                                                                                           aaGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGt                                                    <  7:72494/118‑1 (MQ=255)
                                                                                                                            aGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCgctggggctgg                    >  3:328022/1‑149 (MQ=255)
                                                                                                                            aGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCgctggggctgg                    >  8:66989/1‑149 (MQ=255)
                                                                                                                               cATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGa                 >  7:322793/1‑149 (MQ=255)
                                                                                                                                aTCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGAc                >  7:41403/1‑149 (MQ=255)
                                                                                                                                aTCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCAGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGAc                >  4:64116/1‑149 (MQ=255)
                                                                                                                                 tCCGCCTCATGTAGT‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGGTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGAcc               <  5:297064/149‑1 (MQ=255)
                                                                                                                                 tCCGCCTCAAGGAGT‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGAGCAGTTGTTCAAGGCCCAGCTCGAAGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGAcc               <  5:291891/149‑1 (MQ=255)
                                                                                                                                    gCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATg            >  8:413/1‑149 (MQ=255)
                                                                                                                                         atgcaca‑aaatCCACGTGGTGACATTCGATGCCTAGACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCCCCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGTGGCTGGCGACCATGGTtg        <  4:336135/142‑1 (MQ=255)
                                                                                                                                              ga‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTtggctggc  >  5:244685/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGGTACCGGCGTTGGCGGCGCTGACCGGGATGATCGGCGAAGAGGTCAAAGGCTATGGCTCGCTCAAGCGCGTACCGGCTGAAGCCATGGCCAACGTGCGCAACGACATGTATCTGACCAGTGAAGCCATCCGCCTCATGGAGA‑AAAACCACCTGGTGACATTCGATGCCGACACCCGCAGCCAAGTGCAGTTGTTCAAGGCCCAGCTGGACGACGCCACCCGCTATATTCCGTTGTGGGTGAAGGTGGCTGTAGCCATCGCGCTGGGGCTGGGGACCATGGTTGGCTGGC  >  NC_002947/4638318‑4638608

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: ‑