Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 146,128 T→A 13.7% intergenic (+33/‑115) gltP → / → PP_0138 glutamate/aspartate‑proton DAACS transporter/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947146,1280TA13.7% 131.5 / 13.1 51intergenic (+33/‑115)gltP/PP_0138glutamate/aspartate‑proton DAACS transporter/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base A (4/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

TCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTT  >  NC_002947/145983‑146271
                                                                                                                                                 |                                                                                                                                                
tCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCg                                                                                                                                              >  2:516713/1‑150 (MQ=255)
 cATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgg                                                                                                                                             <  2:464848/150‑1 (MQ=255)
 cATTGCGCGGGGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgg                                                                                                                                             >  2:756182/1‑150 (MQ=255)
   ttGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAGACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgggg                                                                                                                                           <  2:494599/150‑1 (MQ=255)
       tCGGTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACAAGTCGGGGCtct                                                                                                                                       <  1:46933/150‑3 (MQ=255)
       tCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCttt                                                                                                                                       >  1:634965/1‑150 (MQ=255)
       tCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCttt                                                                                                                                       <  2:332758/150‑1 (MQ=255)
             ggAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGTAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTa                                                                                                                                 <  1:735508/150‑1 (MQ=255)
              gAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGTTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTCTTAt                                                                                                                                <  2:332859/150‑1 (MQ=255)
              gAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTAt                                                                                                                                <  1:516713/150‑1 (MQ=255)
                    aggTACGCGGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCCAGAAGAAGGGAGCTGCAGTGGTTGGCGAAACCGCCAAGCGCTGGGCCTGATGGCTGTTTCACCAAGCCCCCCGACCTGTGGGGGCTTTTTGTTATCTGCga                                                                                                                          <  2:365672/148‑1 (MQ=255)
                           acgCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTg                                                                                                                   >  2:13270/1‑150 (MQ=255)
                               ggTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGcc                                                                                                               <  1:13270/150‑1 (MQ=255)
                                gTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgggg                                                                                                                                           <  1:459934/121‑1 (MQ=255)
                                gTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgggg                                                                                                                                           >  2:459934/1‑121 (MQ=255)
                                  cAAGGGCGAGCAGTACTACGCCTCTCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAAAGCCCCGACAAGTCGGGGCtctt                                                                                                                                      >  2:257682/1‑121 (MQ=255)
                                  cAAGGGCGAGCAGTACTACGCCTCTCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAAAGCCCCGACAAGTCGGGGCtctt                                                                                                                                      <  1:257682/124‑1 (MQ=255)
                                     gcgagggCAGTACTACGCTTTTCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTTTTTGACCAAGAGCCCCGACTTGTCGGGGTTTTTTGTTATCTTCGAGATCTTGAGCCGCGCCt                                                                                                         >  2:76277/7‑150 (MQ=255)
                                               tactacGCCTCGCTGATGGTCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTc                                                                                               >  2:635898/1‑150 (MQ=255)
                                                  tacGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTcgcg                                                                                           >  2:220187/1‑150 (MQ=255)
                                                                 gCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCTAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAgg                                                                            <  1:7672/150‑1 (MQ=255)
                                                                      gaAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGc                                                                                      >  2:335008/1‑135 (MQ=255)
                                                                      gaAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGc                                                                                      <  1:335008/135‑1 (MQ=255)
                                                                       aaGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGa                                                                      <  1:787458/150‑1 (MQ=255)
                                                                           aGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGATTACGATCGCCAGTCACAGGCTGCGATACg                                                                  >  1:399786/1‑150 (MQ=255)
                                                                            gggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGc                                                                 >  1:540489/1‑150 (MQ=255)
                                                                            gggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGc                                                                 >  2:316306/1‑150 (MQ=255)
                                                                             ggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACAAGTCGGGGCt                                                                                                                                         <  2:214842/78‑1 (MQ=255)
                                                                             ggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACAAGTCGGGGCt                                                                                                                                         >  1:214842/1‑78 (MQ=255)
                                                                                     aGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCg                                                        <  2:540489/150‑1 (MQ=255)
                                                                                         gTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCa                                                    >  2:454436/1‑150 (MQ=255)
                                                                                              cGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAg                                                                                                   <  2:572618/99‑1 (MQ=255)
                                                                                              cGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAg                                                                                                   >  1:572618/1‑99 (MQ=255)
                                                                                               gAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACATGGCTGTGCGTGCGGCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTc                                              >  2:419035/1‑150 (MQ=255)
                                                                                                     gCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAggggg                                        >  2:221603/1‑150 (MQ=255)
                                                                                                      ncAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCGTGAGCCGCGCCTACAAAAGCATGCCATCACCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGAGGGa                                       >  2:585403/2‑150 (MQ=255)
                                                                                                          gcgcTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTc                                                                                               >  2:319285/1‑91 (MQ=255)
                                                                                                          gcgcTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTc                                                                                               <  1:319285/91‑1 (MQ=255)
                                                                                                            gcTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATAc                                 <  1:756182/150‑1 (MQ=255)
                                                                                                                aGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGAt                             >  2:721338/1‑150 (MQ=255)
                                                                                                                        ggCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACg                     <  2:634965/150‑1 (MQ=255)
                                                                                                                         gCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTAACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACgg                    >  2:307167/1‑150 (MQ=255)
                                                                                                                             tttGACCACGAGCCCTGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACATGAGCGCTCAATAAGCTGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCtt                <  1:454436/150‑1 (MQ=255)
                                                                                                                              ttGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCttt               <  1:316306/150‑1 (MQ=255)
                                                                                                                                 aCCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACtttt                                               >  1:565965/1‑115 (MQ=255)
                                                                                                                                 aCCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACtttt                                               <  2:565965/115‑1 (MQ=255)
                                                                                                                                  caaaaagCCCCGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCt           >  1:283884/6‑150 (MQ=255)
                                                                                                                                   aaaaagCCCCGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTc          >  1:141328/5‑150 (MQ=255)
                                                                                                                                           ccGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCttttt  >  2:582680/1‑150 (MQ=255)
                                                                                                                                            cGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCgg                                                       >  2:5013/1‑96 (MQ=255)
                                                                                                                                            cGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCgg                                                       <  1:5013/96‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTT  >  NC_002947/145983‑146271

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