Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,138,434 A→G intergenic (‑249/‑68) PP_3641 ← / → PP_3643 cytosine/purine/uracil/thiamine/allantoin permease family protein/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,138,4340AG97.8% 143.4 / ‑4.6 46intergenic (‑249/‑68)PP_3641/PP_3643cytosine/purine/uracil/thiamine/allantoin permease family protein/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/21);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAG  >  NC_002947/4138289‑4138573
                                                                                                                                                 |                                                                                                                                           
tGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCgg                                                                                                                                          >  7:212986/1‑149 (MQ=255)
tGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCgg                                                                                                                                          >  5:77124/1‑149 (MQ=255)
   tCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGta                                                                                                                                        >  1:118239/1‑148 (MQ=255)
       cATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCa                                                                                                                                   <  3:64527/149‑1 (MQ=255)
           aaaaaaCGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGATTGCAAAATAGGTATTGAGGTATACAAt                                                                                                                                  >  5:393339/1‑146 (MQ=255)
              aaaCGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGaca                                                                                                                            >  4:198464/1‑149 (MQ=255)
                    ccACTGATGCTCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCATTGCCACGTTTTTCGCGGCTTTGGTCCGTTCATCGCGAATACTCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCaa                                                                                                                      <  6:393339/149‑1 (MQ=255)
                    ccACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCaa                                                                                                                      >  4:359057/1‑149 (MQ=255)
                     cACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAAc                                                                                                                     >  5:375489/1‑149 (MQ=255)
                     cACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAAc                                                                                                                     >  6:49070/1‑149 (MQ=255)
                       cTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGa                                                                                                                   <  6:316341/149‑1 (MQ=255)
                          aTCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGATGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATa                                                                                                                 >  7:303002/1‑148 (MQ=255)
                             cGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATaaaaa                                                                                                             >  5:369195/1‑149 (MQ=255)
                             cGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATaaaaa                                                                                                             <  2:118239/149‑1 (MQ=255)
                              gCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAAc                                                                                                            >  3:247884/1‑149 (MQ=255)
                                               gACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGcc                                                                                           <  6:375489/149‑1 (MQ=255)
                                                     ccGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTAATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGACCTTCa                                                                                      >  7:289438/1‑148 (MQ=255)
                                                       gCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCagag                                                                                   <  3:359057/149‑1 (MQ=255)
                                                       gCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCagag                                                                                   <  3:198464/149‑1 (MQ=255)
                                                                 tCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCga                                                                         >  8:52278/1‑149 (MQ=255)
                                                                   aaGGCCACGTTTTTCGCGGCGTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTCCAGAGGTACTCCCgat                                                                        <  4:247884/148‑1 (MQ=255)
                                                                      gCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATc                                                                    >  4:203245/1‑149 (MQ=255)
                                                                        cACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGa                                                                  >  2:325113/1‑149 (MQ=255)
                                                                         aCGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGAc                                                                 <  3:203245/149‑1 (MQ=255)
                                                                         aCGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGAc                                                                 <  6:77124/149‑1 (MQ=255)
                                                                              tttCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCaa                                                             <  8:212986/148‑1 (MQ=255)
                                                                                  gcgGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTc                                                        >  4:237304/1‑149 (MQ=255)
                                                                                     gCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTa                                                      >  2:321219/1‑148 (MQ=255)
                                                                                     gCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTAc                                                     <  8:303002/149‑1 (MQ=255)
                                                                                          gggCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGaa                                                >  4:261442/1‑149 (MQ=255)
                                                                                          gggCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGa                                                 <  5:380386/148‑1 (MQ=255)
                                                                                          gggCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGa                                                 <  6:110083/148‑1 (MQ=255)
                                                                                          gggCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGa                                                 <  3:219500/148‑1 (MQ=255)
                                                                                            gCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAAc                                               <  3:237304/148‑1 (MQ=255)
                                                                                                ttCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCt                                          <  7:189419/149‑1 (MQ=255)
                                                                                                   aTCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACCTTCTggg                                       >  8:130942/1‑149 (MQ=255)
                                                                                                               cAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTccc                              >  6:252476/1‑146 (MQ=255)
                                                                                                               cAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTccc                              <  5:252476/146‑1 (MQ=255)
                                                                                                                   gTTGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGc                       <  7:52278/149‑1 (MQ=255)
                                                                                                                    ttGAAAATTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCg                      >  4:163483/1‑149 (MQ=255)
                                                                                                                       aaaaTTCCGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGtt                    <  1:321219/148‑1 (MQ=255)
                                                                                                                              cGCTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGtcat            <  6:311924/149‑1 (MQ=255)
                                                                                                                                cTTGCAAAATAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGtcatca          >  5:63249/1‑149 (MQ=255)
                                                                                                                                     aaaaTAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTc                                <  6:343565/122‑1 (MQ=255)
                                                                                                                                     aaaaTAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTc                                >  5:343565/1‑122 (MQ=255)
                                                                                                                                        aTAGGTATTGCGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCAGATGATCGACGCAACTGTCTACAACAACGATCTCGGCACATTCCCGTCCGGCGTTACCGTCATCACAACCCAg  >  5:159/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
TGATCCGCATGAAAAAACGGCCACTGATCCGCTTGCAAAACCCCTGCGACAAACCGCCTCACCGCTCAAGGCCACGTTTTTCGCGGCTTTGGGCCGTTCATCGCGAAAACGCAGGGTTGAAAATTCCGCTTGCAAAATAGGTATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAG  >  NC_002947/4138289‑4138573

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