Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 146,127 T→A 9.0% intergenic (+32/‑116) 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,1270TA 9.0% 123.3 / 2.8 43intergenic (+32/‑116)gltP/PP_0138glutamate/aspartate‑proton DAACS transporter/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base A (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 = 9.93e-01

ATTGGTCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTTCCTG  >  NC_002947/145978‑146275
                                                                                                                                                     |                                                                                                                                                    
aTTGGTCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACt                                                                                                                                                      >  2:353197/1‑150 (MQ=255)
 ttGGTCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACtt                                                                                                                                                     <  2:126424/150‑1 (MQ=255)
      cATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTTTCgg                                                                                                                                                <  2:866847/150‑1 (MQ=255)
       aTTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACAAGTCggg                                                                                                                                               >  1:300523/1‑150 (MQ=255)
        ttGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgggg                                                                                                                                              >  1:575632/1‑150 (MQ=255)
        ttGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCgggg                                                                                                                                              >  1:237085/1‑150 (MQ=255)
         tGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGc                                                                                                                                             >  2:126459/1‑150 (MQ=255)
              gCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCttttt                                                                                                                                        <  2:237085/150‑1 (MQ=255)
              gCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACAAGTCGGGGCtcttg                                                                                                                                        <  1:895745/150‑2 (MQ=255)
              gCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACAAGTCGGGGCtcttg                                                                                                                                        <  2:300523/150‑2 (MQ=255)
                tGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTga                                                                                                                                      <  2:826318/150‑2 (MQ=255)
                 gggAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGtt                                                                                                                                     <  1:353197/150‑1 (MQ=255)
                     aGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGCCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCt                                                                                                                                 <  1:686533/150‑1 (MQ=255)
                          tGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTa                                                                                                                                    >  1:915112/1‑142 (MQ=255)
                          tGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTa                                                                                                                                    <  2:915112/142‑1 (MQ=255)
                            tACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAGGCGCTGAGTCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGAt                                                                                                                          >  1:335786/1‑150 (MQ=255)
                             acgacgCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCg                                                                                                                                                 <  2:966348/126‑1 (MQ=255)
                             acgacgCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCg                                                                                                                                                 >  1:966348/1‑126 (MQ=255)
                                    ggTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAACCCGCCACGCGCCGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGcc                                                                                                                  <  2:625747/150‑1 (MQ=255)
                                                gCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAg                                                                                                      <  1:1011053/150‑1 (MQ=255)
                                                    tactacGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGTGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCGTACAAAGGCTc                                                                                                  >  1:489065/1‑150 (MQ=255)
                                                    tactacGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTc                                                                                                  <  1:785387/150‑1 (MQ=255)
                                                       tacGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCACGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTcgcg                                                                                               >  2:602817/1‑150 (MQ=255)
                                                          gCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCACCGCCTAGTCCGGGTTCTGCTTGATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATc                                                                                            >  1:271098/1‑150 (MQ=255)
                                                             tCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGCATGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGcc                                                                                         >  1:1002096/1‑150 (MQ=255)
                                                                    tGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTcaca                                                                                  <  1:126459/150‑1 (MQ=255)
                                                                     ggCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAg                                                                                 >  1:586976/1‑150 (MQ=255)
                                                                      gCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAgg                                                                                >  2:660687/1‑150 (MQ=255)
                                                                           gaAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCg                                                                                                                              >  1:313916/1‑99 (MQ=255)
                                                                           gaAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCg                                                                                                                              <  2:313916/99‑1 (MQ=255)
                                                                           gaAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCg                                                                           <  2:38368/150‑1 (MQ=255)
                                                                                         tAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTc                                                             >  1:357823/1‑150 (MQ=255)
                                                                                                         cGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGAGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAgggg                                             <  2:586976/150‑1 (MQ=255)
                                                                                                             aaGCGCTTAGCCTGATGTCTGTTTGACCAAGAGCGCCGACTTGTCGTGGCTTTTTGTTATCTGCAAGATCTTGAGCCGCGCCTACAAAGGCTCGGGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAAc                                         <  2:489065/150‑1 (MQ=255)
                                                                                                                 gcTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATAc                                     <  1:1127053/150‑1 (MQ=255)
                                                                                                                    gAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGa                                  <  1:374533/150‑1 (MQ=255)
                                                                                                                     aGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGAt                                 <  2:575632/150‑1 (MQ=255)
                                                                                                                       ccTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGc                               <  2:1002096/150‑1 (MQ=255)
                                                                                                                            tGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGTGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGACAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGTGGTATCATCCGGATGCTCAAc                          >  2:90585/1‑150 (MQ=255)
                                                                                                                             ggCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACg                         >  2:1072612/1‑150 (MQ=255)
                                                                                                                                  tttGACCAAGAGCACCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCtt                    >  2:50183/1‑150 (MQ=255)
                                                                                                                                    tGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTg                  >  2:694217/1‑150 (MQ=255)
                                                                                                                                       ccAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCt               >  1:1128992/1‑150 (MQ=255)
                                                                                                                                       caaaaagCCCCGACAAGTCGGGGCTCTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCt               >  1:785968/6‑150 (MQ=255)
                                                                                                                                                    cTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTTCCtg  <  1:1072612/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTGGTCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTTCCTG  >  NC_002947/145978‑146275

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