Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,563,828 A→C 17.1% intergenic (+94/+28) PP_1372 → / ← pitB hypothetical protein/phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,563,8280AC17.1% 98.6 / 15.9 46intergenic (+94/+28)PP_1372/pitBhypothetical protein/phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (24/14);  new base C (3/5);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

TTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGCAGGCCACT  >  NC_002947/1563687‑1563972
                                                                                                                                             |                                                                                                                                                 
ttGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCTCGGGCGcc                                                                                                                                            <  3:70182/149‑1 (MQ=255)
   gTGCGCGGGGTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTAGGAAAAAGGCGCCCGAGGGCGCCtt                                                                                                                                          >  6:370001/1‑148 (MQ=255)
     gcgcGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCTCGGGCGCCttttt                                                                                                                                       <  4:212638/149‑4 (MQ=255)
       gcgGGTTGAGGGGGTCGTTGTTAGGGCGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCTCGGGCGCCtttttca                                                                                                                                     <  5:269825/149‑6 (MQ=255)
          ggTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTcattcttt                                                                                                                                  <  7:280531/149‑1 (MQ=255)
                tGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGgtg                                                                                                                            >  7:318358/1‑149 (MQ=255)
                tGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGgtg                                                                                                                            >  7:300922/1‑149 (MQ=255)
                       gttgttAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACt                                                                                                                       >  6:255380/1‑147 (MQ=255)
                                 ggCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAAt                                                                                                               <  1:288226/145‑1 (MQ=255)
                                 ggCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAAt                                                                                                               >  2:288226/1‑145 (MQ=255)
                                   cAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCttt                                                                                                         <  1:375198/149‑1 (MQ=255)
                                            gTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCACGCTCAGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACgccag                                                                                               >  5:121548/1‑149 (MQ=255)
                                             tAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACgccagc                                                                                              <  8:188183/149‑1 (MQ=255)
                                               gTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACgccagc                                                                                              >  5:44836/1‑147 (MQ=255)
                                                tCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGa                                                                                                                                                   >  1:166194/1‑94 (MQ=255)
                                                tCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGa                                                                                                                                                   <  2:166194/94‑1 (MQ=255)
                                                      tGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAgg                                                                                     >  6:124555/1‑149 (MQ=255)
                                                                cATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGgccgc                                                                             >  3:289042/1‑147 (MQ=255)
                                                                cATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGgccgc                                                                             >  8:34794/1‑147 (MQ=255)
                                                                     ccTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGATATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGcc                                                                       >  1:58652/1‑148 (MQ=255)
                                                                       ttGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCAt                                                                    >  8:181368/1‑149 (MQ=255)
                                                                         ggTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATc                                                                   <  6:61508/148‑1 (MQ=255)
                                                                           tAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGag                                                                >  3:165537/1‑149 (MQ=255)
                                                                             gCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGagg                                                               >  2:259775/1‑148 (MQ=255)
                                                                               acaAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcag                                                            >  8:239060/1‑149 (MQ=255)
                                                                               acaAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcag                                                            >  6:113347/1‑149 (MQ=255)
                                                                                caAGCCCGCTCCGGCAAAGAACAAGGTGTCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcagg                                                           <  8:300922/149‑1 (MQ=255)
                                                                                            ggCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGCGGGCGCCTTCATTCAGTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCCGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCCGGCAGGGTCAGGAc                                               >  8:256331/1‑149 (MQ=255)
                                                                                               aaaGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCa                                            >  5:163392/1‑149 (MQ=255)
                                                                                                  gAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGc                                         <  5:255380/149‑1 (MQ=255)
                                                                                                    aCAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCca                                       <  4:289042/149‑1 (MQ=255)
                                                                                                        ggTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCcagcag                                   >  8:186148/1‑149 (MQ=255)
                                                                                                        ggTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCcagcag                                   >  6:298939/1‑149 (MQ=255)
                                                                                                             cGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATACATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATgg                              >  3:229269/1‑149 (MQ=255)
                                                                                                                 ggCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTTTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGTGTCAGGACCCATGCCAGCAGGATGGtttt                          <  6:121548/149‑1 (MQ=255)
                                                                                                                 ggCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGtttt                          <  5:113347/149‑1 (MQ=255)
                                                                                                                 ggCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGttt                           <  5:138862/148‑1 (MQ=255)
                                                                                                                    aTAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                        >  7:97803/1‑148 (MQ=255)
                                                                                                                     tAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTcaca                      <  8:97803/149‑1 (MQ=255)
                                                                                                                     tAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTcaca                      <  6:163392/149‑1 (MQ=255)
                                                                                                                              tGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCt             >  5:57782/1‑149 (MQ=255)
                                                                                                                                aaaaaGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCtt            <  7:181368/148‑1 (MQ=255)
                                                                                                                                aaaaaGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTg           >  3:335341/1‑149 (MQ=255)
                                                                                                                                  gaaGGCGCCCTCGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                        <  6:290812/133‑1 (MQ=255)
                                                                                                                                  gaaGGCGCCCTCGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                        >  5:290812/2‑134 (MQ=255)
                                                                                                                                  gaaGGCGCCCTCGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                        <  2:75656/133‑1 (MQ=255)
                                                                                                                                  gaaGGCGCCCTCGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                        >  1:75656/2‑134 (MQ=255)
                                                                                                                                          ccGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGAACGCCAGCCAGAACAGCCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGCAGGCCACt  >  1:275112/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
TTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGG‑ACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGCAGGCCACT  >  NC_002947/1563687‑1563972

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