Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,229,227 G→A 12.2% intergenic (‑121/+55) PP_4606 ← / ← PP_4607 ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,229,2270GA12.2% 99.4 / 6.3 41intergenic (‑121/+55)PP_4606/PP_4607ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (18/18);  new base A (2/3);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

GGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCG  >  NC_002947/5229081‑5229376
                                                                                                                                                  |                                                                                                                                                      
gggACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCgg                                                                                                                                                     <  2:458960/150‑1 (MQ=255)
     atctcgcgGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTaaa                                                                                                                                                <  1:608726/146‑1 (MQ=255)
               gTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTa                                                                                                                                      >  2:568782/1‑150 (MQ=255)
                tGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTAc                                                                                                                                     >  2:849920/1‑150 (MQ=255)
                    cTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgtgg                                                                                                                                <  1:849920/150‑1 (MQ=255)
                              gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTAc                                                                                                                         >  2:362552/1‑147 (MQ=39)
                              gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGGTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTAc                                                                                                                         <  1:362552/147‑1 (MQ=39)
                               cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                     <  2:850385/150‑1 (MQ=255)
                               cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                     <  2:291179/150‑1 (MQ=255)
                                  aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTATGGGCGGGTTTACCCGCGaa                                                                                                                  <  2:76662/150‑1 (MQ=38)
                                    cTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAaga                                                                                                                >  1:314279/1‑150 (MQ=255)
                                        gAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCt                                                                                                            <  1:568782/150‑1 (MQ=255)
                                          gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTgg                                                                                                          <  2:66705/150‑1 (MQ=255)
                                          gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCgc                                                                                                                     <  2:200493/139‑2 (MQ=37)
                                          gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCgc                                                                                                                     >  1:200493/1‑138 (MQ=37)
                                             ggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGcac                                                                                                       >  1:242219/1‑150 (MQ=255)
                                                    attGCGAATCAATCGCACTTTCACTCATTACACGGCTGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgggg                                                                                                                                <  2:390221/118‑1 (MQ=255)
                                                    attGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgggg                                                                                                                                >  1:390221/1‑118 (MQ=255)
                                                           atcaatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGt                                                                                                                             <  2:716803/114‑1 (MQ=255)
                                                           atcaatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGt                                                                                                                             >  1:716803/1‑114 (MQ=255)
                                                                tcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                    <  1:773586/150‑1 (MQ=255)
                                                                  gCACTTTCACTCATTACACGGATGTCACGTGCTGTTACCTACCATAATAGCGGAAAGATGTTTGCTGTGGCGGACTCTGCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGt                                                                                  <  1:406072/150‑1 (MQ=255)
                                                                     cTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCGCGaa                                                                                                                  >  1:720223/1‑115 (MQ=25)
                                                                     cTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCGCGaa                                                                                                                  <  2:720223/115‑1 (MQ=25)
                                                                           cTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCa                                                                                      >  2:590886/1‑137 (MQ=255)
                                                                           cTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCa                                                                                      <  1:590886/137‑1 (MQ=255)
                                                                               ttACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAac                                                                     >  1:504124/1‑150 (MQ=255)
                                                                                 acacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAacac                                                                   >  1:202693/1‑150 (MQ=255)
                                                                                   acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTAACCGTGAAGAGGATGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACgg                                                                 >  1:57714/1‑150 (MQ=255)
                                                                                                gCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAAGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaaa                                                    >  1:540559/1‑150 (MQ=255)
                                                                                                   gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTg                                                                                                           >  2:96458/1‑92 (MQ=255)
                                                                                                   gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                        <  2:316166/127‑1 (MQ=255)
                                                                                                   gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                        >  1:316166/1‑127 (MQ=255)
                                                                                                   gTTACCTCACCTGTGCCGGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAGCCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTg                                                                                                           <  1:96458/92‑1 (MQ=255)
                                                                                                          taCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                    >  1:166228/2‑92 (MQ=255)
                                                                                                          cACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                    <  2:166228/92‑1 (MQ=255)
                                                                                                           aCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCt                                                                                                            <  2:386834/83‑1 (MQ=255)
                                                                                                           aCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCt                                                                                                            >  1:386834/1‑83 (MQ=255)
                                                                                                                      gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa                              >  2:772670/1‑150 (MQ=255)
                                                                                                                        aaTATTTTTGCTGTGGCGGCCTCTTCGCGGGAAAAGCCGCCCCTACGGGGGGTTTAACCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                            >  1:112370/1‑150 (MQ=255)
                                                                                                                                           ccTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCg         <  2:202693/150‑1 (MQ=255)
                                                                                                                                                  gcgGGGAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCAccc  >  2:541376/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GGGACGCAACGCGGGGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCG  >  NC_002947/5229081‑5229376

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