Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,803 G→C 13.0% intergenic (+69/+73) PP_2589 → / ← PP_2590 aldehyde dehydrogenase family protein/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,958,8030GC13.0% 65.8 / 5.3 31intergenic (+69/+73)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (15/11);  new base C (2/2);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

GATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACC  >  NC_002947/2958666‑2958950
                                                                                                                                         |                                                                                                                                                   
gATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACAGCGTCGTCTGCTTCCCGGTCATGCCCGCTCCCACAGGTTCAGCACAATTCTTATGGTTTTGCt                                                                                                                                          >  5:110456/1‑149 (MQ=255)
      tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAcc                                                                                                                                     <  7:83895/148‑1 (MQ=25)
        aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                  <  3:241086/149‑1 (MQ=255)
          tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                >  1:251972/1‑149 (MQ=255)
            cATTCCTTCGACACCTATTCGCAACTCAAGACCACCTGGCACCAGTTGCGGGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCATCACAAGTCTGAGGGTTGGGCTGGACCTGTgg                                                                                                                                <  8:32190/147‑1 (MQ=38)
              ttcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGc                                                                                                                            >  8:225599/1‑149 (MQ=255)
                              tCGCAACTCACGACCACGTGGTACCAGTTGCGCGGCTTATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAGCACAAGTCTGTGGGTTGTGCTGGACCTGTGGGAGGGGGCGCGCCCGCGaa                                                                                                             <  2:56259/148‑1 (MQ=39)
                              tCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTAGCGGGCATTCCCCCACCCACAGGTTAAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcccg                                                                                                               >  1:56259/1‑146 (MQ=255)
                                  aaCTAAAGACATCGTGTTACTTGTTGGGGGGTTTATGTTTTGATCTTTGACCGCGTCGTCTGCTTCGCGGGAATGCCCGCTCCCACATGTTCTGCATATGTCTGTGGGTTGTGCTGGACGTGTGGTAGCGGGCGCGCCCGCGAAGAgg                                                                                                         <  6:110456/148‑1 (MQ=255)
                                                    aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGTGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                      >  5:188961/1‑149 (MQ=255)
                                                                gCTGATGTATTGATCATTGCCCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCCCAGGTGCCCCACAAGTCTGAGGGGTGGGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                          <  7:225599/149‑1 (MQ=38)
                                                                           gATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACg                                                                >  2:197369/1‑148 (MQ=255)
                                                                                                          aTGCCCGCTCCCACAGGTTCAGCTCAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGGGGGACCCGGTTAGAAGTCCACGCTAGCGGAAATGTGATAGGTAGGAGGAGcccc                                >  3:136072/1‑149 (MQ=38)
                                                                                                             cccGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                                 <  1:78382/129‑1 (MQ=255)
                                                                                                             cccGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                                 >  2:78382/1‑129 (MQ=255)
                                                                                                                 cTCCCACAGGTTCACCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAcccc                         <  1:197369/149‑1 (MQ=255)
                                                                                                                    ccACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                      >  5:128347/1‑149 (MQ=25)
                                                                                                                     cacaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                     >  5:242867/1‑149 (MQ=255)
                                                                                                                      acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGt                    >  3:24262/1‑149 (MQ=255)
                                                                                                                      acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCCCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGt                    >  1:49037/1‑149 (MQ=255)
                                                                                                                      acaGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                     >  6:274472/1‑148 (MQ=37)
                                                                                                                        aGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCt                  <  4:142539/149‑1 (MQ=255)
                                                                                                                        aGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCt                  >  3:142539/1‑149 (MQ=255)
                                                                                                                          gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                <  6:128347/146‑1 (MQ=25)
                                                                                                                                tacaAGTGTGTGGGTTGTGCTGGAGCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccc                                 <  4:67858/125‑1 (MQ=39)
                                                                                                                                cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgcgc                                                                                            <  2:194646/67‑1 (MQ=255)
                                                                                                                                cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgcgc                                                                                            >  1:194646/1‑67 (MQ=255)
                                                                                                                                cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccc                                 >  3:67858/1‑126 (MQ=255)
                                                                                                                                cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                <  4:26725/143‑1 (MQ=255)
                                                                                                                                cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                >  3:26725/1‑143 (MQ=255)
                                                                                                                                        tGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAcc  >  7:214804/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
GATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACC  >  NC_002947/2958666‑2958950

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