Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,580,975 G→T 9.3% intergenic (‑60/‑42) serB ← / → PP_4910 phosphoserine phosphatase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,580,9750GT 9.3% 127.8 / 3.9 43intergenic (‑60/‑42)serB/PP_4910phosphoserine phosphatase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base T (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

GCAGTGAGACCGGGACGGTCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCACCAACATT  >  NC_002947/5580827‑5581124
                                                                                                                                                    |                                                                                                                                                      
gCAGTGAGACCGGGACGGTCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGc                                                                                                                                                       <  2:195384/150‑1 (MQ=255)
         ccGGGACGGTCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTccc                                                                                                                                              >  1:36353/1‑150 (MQ=255)
          cGGGACGGTCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTcccc                                                                                                                                             <  2:540168/150‑1 (MQ=255)
          cGGGACGGTCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTcccc                                                                                                                                             >  1:442710/1‑150 (MQ=255)
                  tCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTc                                                                                                                                     >  1:177792/1‑150 (MQ=255)
                    cTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCAt                                                                                                                                   <  2:901007/150‑1 (MQ=255)
                         cGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCAt                                                                                                                              <  2:36353/150‑1 (MQ=255)
                              aTGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAaga                                                                                                                         <  2:799008/150‑1 (MQ=255)
                                gTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGc                                                                                                                       >  1:506183/1‑150 (MQ=255)
                                       aGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCTCTATACTCCCCCCACTTTCATGCCATTAAGAGc                                                                                                                       <  2:627792/143‑1 (MQ=255)
                                       aGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCTCTATACTCCCCCCACTTTCATGCCATTAAGAGc                                                                                                                       <  2:627338/143‑1 (MQ=255)
                                       aGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCTCTATACTCCCCCCACTTTCATGCCATTAAGAGc                                                                                                                       >  1:627792/1‑143 (MQ=255)
                                       aGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCTCTATACTCCCCCCACTTTCATGCCATTAAGAGc                                                                                                                       >  1:627338/1‑143 (MQ=255)
                                         gACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCtg                                                                                                              >  1:92246/1‑150 (MQ=255)
                                          aCGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCtgt                                                                                                             <  2:92246/150‑1 (MQ=255)
                                                                tCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCaaa                                                                                      >  2:160726/1‑150 (MQ=255)
                                                                      ttCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCAt                                                                                                                                   <  1:57482/100‑1 (MQ=255)
                                                                      ttCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCAt                                                                                                                                   >  2:57482/1‑100 (MQ=255)
                                                                      ttCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGAc                                                                                >  1:707080/1‑150 (MQ=255)
                                                                            gAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATttc                                                                          <  2:442710/150‑1 (MQ=255)
                                                                                      cATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCtgatga                                                                <  1:160726/150‑1 (MQ=255)
                                                                                              cGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATcgc                                                        >  1:641138/1‑150 (MQ=255)
                                                                                                aTTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCATCCTGATGATCTATCgcgc                                                      >  1:62130/1‑150 (MQ=255)
                                                                                                       aTGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGt                                               <  2:62130/150‑1 (MQ=255)
                                                                                                        tGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTc                                              <  2:641138/150‑1 (MQ=255)
                                                                                                           ccATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGc                                                                                                                                 >  1:492619/1‑65 (MQ=255)
                                                                                                           ccATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGc                                                                                                                                 <  2:492619/65‑1 (MQ=255)
                                                                                                           ccATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAAc                                           <  1:860055/150‑1 (MQ=255)
                                                                                                              tCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAAcgtc                                        <  2:707080/150‑1 (MQ=255)
                                                                                                                      cAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCGAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAAcgtcgt                                      >  2:410581/1‑144 (MQ=255)
                                                                                                                      cAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCGAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAAcgtcgt                                      <  1:410581/144‑1 (MQ=255)
                                                                                                                       aaCGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGc                               >  1:393388/1‑150 (MQ=255)
                                                                                                                       aaCGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGc                               >  1:1107259/1‑150 (MQ=255)
                                                                                                                          gATACTTTGCTCACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATCACG‑GCTGCGCCCGGTCGAGCCGGCGAACGCCAGTCAAACCAGACAATTCCTTCCTGATGAGCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTgg                            <  1:1177845/150‑1 (MQ=255)
                                                                                                                               tttGCCGACAAAACCCGGCTCGCTATACTTCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTg                       >  1:280881/1‑150 (MQ=255)
                                                                                                                                    cGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTgc                      >  2:1029842/1‑146 (MQ=255)
                                                                                                                                    cGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTgc                      <  1:1029842/146‑1 (MQ=255)
                                                                                                                                          aaCCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTg            <  2:506183/150‑1 (MQ=255)
                                                                                                                                             ccGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCAc         >  1:268512/1‑150 (MQ=255)
                                                                                                                                               ggCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGa                                                                                 <  2:238743/76‑1 (MQ=255)
                                                                                                                                               ggCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGa                                                                                 >  1:238743/1‑76 (MQ=255)
                                                                                                                                                gCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCACCaa      <  2:393388/150‑1 (MQ=255)
                                                                                                                                                  tCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCACCAACa    >  2:110469/1‑150 (MQ=255)
                                                                                                                                                    gCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCACCAACAtt  >  1:801106/1‑150 (MQ=255)
                                                                                                                                                    gCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCACCAACAtt  >  2:709585/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCAGTGAGACCGGGACGGTCCTCACCGGTGATGTTGATCAGGACGATTTCGCGCACAACCTGGCTCCGACTTCGATGAAAAATGCGCATTCTACCGATTTTCCATGACCATCGGGCGCCAACGATACTTTGCCGACAAAACCCGGCTCGCTATACTCCCCCCACTTTCATGCCATTAAGAGCCGAGCTCTGT‑GAACCGGCCCACGCCAGTCAAACCAGACAATTTCTTCCTGATGATCTATCGCGCCCTGAGTCAACGTCGTGTGCCGCTGGCACTGCGCATCGCCTGCACCAACATT  >  NC_002947/5580827‑5581124

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