Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,500,613 G→A 8.2% intergenic (+17/+8) PP_2193 → / ← PP_2194 ferric siderophore receptor/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,500,6130GA 8.2% 139.3 / 3.1 49intergenic (+17/+8)PP_2193/PP_2194ferric siderophore receptor/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGCCCCT  >  NC_002947/2500465‑2500751
                                                                                                                                                    |                                                                                                                                          
gACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACgg                                                                                                                                           >  2:69015/1‑150 (MQ=255)
     gCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTg                                                                                                                                      >  2:1129115/1‑150 (MQ=255)
         tACAGCACGGTATGGGGACAGCGTAGGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCa                                                                                                                                  >  1:465716/1‑150 (MQ=255)
          aCAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAg                                                                                                                                 <  1:69015/150‑1 (MQ=255)
           cAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAgg                                                                                                                                <  2:309293/150‑1 (MQ=255)
                    aTGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAGGcc                                                                                                                              <  2:121745/143‑2 (MQ=255)
                    aTGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAGGcc                                                                                                                              >  1:121745/1‑142 (MQ=255)
                       gggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCa                                                                                                                    >  1:42879/1‑150 (MQ=255)
                           cAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAg                                                                                                                                 >  2:910108/1‑133 (MQ=255)
                           cAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAg                                                                                                                                 <  1:910108/133‑1 (MQ=255)
                                tGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCt                                                                                                           >  1:848617/1‑150 (MQ=255)
                                     cGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCa                                                                                                      >  1:364384/1‑150 (MQ=255)
                                             tATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATg                                                                                              >  1:391514/1‑150 (MQ=255)
                                              atacGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTc                                                                                                                                   >  1:463465/4‑112 (MQ=255)
                                              atacGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTc                                                                                                                                   <  2:463465/109‑1 (MQ=255)
                                                    cGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGcc                                                                                       >  2:631767/1‑150 (MQ=255)
                                                        cTACGGGCGGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTGCTGATTGTGGTGGCCTGAACGGGCCTGTTAAGGCAAGAGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGc                                                                                   >  2:159097/1‑150 (MQ=255)
                                                        cTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGc                                                                                   <  1:772923/150‑1 (MQ=255)
                                                         tACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGcc                                                                                  >  2:601684/1‑150 (MQ=255)
                                                                 ggCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACgccgcc                                                                                                                     <  2:1180262/107‑1 (MQ=255)
                                                                 ggCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACgccgcc                                                                                                                     >  1:1180262/1‑107 (MQ=255)
                                                                 ggCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCt                                                                          <  1:890793/150‑1 (MQ=255)
                                                                   cGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCtgt                                                                        <  1:631767/150‑1 (MQ=255)
                                                                          gTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACgg                                                                                                                                           <  2:452105/76‑1 (MQ=255)
                                                                          gTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACgg                                                                                                                                           >  1:452105/1‑76 (MQ=255)
                                                                          gTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATgg                                                                                             <  2:558984/122‑1 (MQ=255)
                                                                          gTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATgg                                                                                             >  1:558984/1‑122 (MQ=255)
                                                                           tATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTg                                                                <  2:364384/150‑1 (MQ=255)
                                                                               aCTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCTTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCgg                                                            <  2:465716/150‑1 (MQ=255)
                                                                                   ccAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCa                                                        >  1:403491/1‑150 (MQ=255)
                                                                                    cAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAg                                                       <  2:403491/150‑1 (MQ=255)
                                                                                     aGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAgg                                                      <  1:1129115/150‑1 (MQ=255)
                                                                                           gTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAAc                                                <  1:360195/150‑1 (MQ=255)
                                                                                           gTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAAc                                                <  2:42879/150‑1 (MQ=255)
                                                                                                 ggACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCAcc                                          <  2:848617/150‑1 (MQ=255)
                                                                                                     cTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCaa                                      >  1:7029/1‑150 (MQ=255)
                                                                                                       ataCGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCaaaa                                    >  1:1022629/1‑150 (MQ=255)
                                                                                                          cGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAACTCa                                 >  1:620559/1‑150 (MQ=255)
                                                                                                          cGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCa                                 >  1:813909/1‑150 (MQ=255)
                                                                                                           gCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCaa                                <  1:601684/150‑1 (MQ=255)
                                                                                                            cTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAg                               <  2:813909/150‑1 (MQ=255)
                                                                                                                  tGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAg                                                                                                                                 >  1:829645/1‑46 (MQ=25)
                                                                                                                  tGGAGCATGGCTTACTGATTGTGGTGGCCTGAACAGGCCCGTTCAg                                                                                                                                 <  2:829645/46‑1 (MQ=25)
                                                                                                                     aGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTc                      <  2:7029/150‑1 (MQ=255)
                                                                                                                                    ttGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCgg                                                            <  1:735851/97‑1 (MQ=255)
                                                                                                                                    ttGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCgg                                                            >  2:735851/1‑97 (MQ=255)
                                                                                                                                     tgtgGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGc      <  1:258559/150‑1 (MQ=255)
                                                                                                                                         gtggCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAGATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCACATAACAGCCCCt  >  2:438153/1‑150 (MQ=255)
                                                                                                                                              cTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCaaa                                                                                                    >  2:97359/1‑47 (MQ=255)
                                                                                                                                              cTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCaaa                                                                                                    <  1:97359/47‑1 (MQ=255)
                                                                                                                                              cTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGc                                                                                        >  1:199936/1‑59 (MQ=255)
                                                                                                                                              cTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGc                                                                                        <  2:199936/59‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCGCAAATACATGGGCTGCCTCGCCAACCCCCTGTGCCCGGTGCCGGGCCAGGGTGAACGGCACCGCAAAATCAAGCTCAGGGTCCAGCGCCCATAACAGCCCCT  >  NC_002947/2500465‑2500751

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