Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,680,476 A→T 12.5% A274A (GCT→GCA) ‡ PP_2350 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,680,4760AT12.5% 83.0 / 5.0 32A274A (GCT→GCA) ‡PP_2350LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (11/17);  new base T (2/2);  total (13/19)
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.67e-01

GGTGGGGGTGATCACCCAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTGAACTCCGGTGCCTGCGGTACCCGCTCGAACACGCCGCTGTCCAGATA  >  NC_002947/2680361‑2680610
                                                                                                                   |                                                                                                                                      
ggtggGGGTGATCACCCAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGtt                                                                                                      <  1:349647/150‑1 (MQ=255)
    ggggTGATCACCCAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAg                                                                                                  <  2:827820/150‑1 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCTGGTAGCTGCCGgcgc                                                                                                                         >  2:636329/1‑116 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCTGGTAGCTGCCGgcgc                                                                                                                         >  2:636324/1‑116 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCTGGTAGCTGCCGgcgc                                                                                                                         <  1:636329/116‑1 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCTGGTAGCTGCCGgcgc                                                                                                                         <  1:636324/116‑1 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGt                                                                                       >  1:112011/1‑150 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGt                                                                                       >  1:75041/1‑150 (MQ=255)
               ccAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGt                                                                                       <  1:161574/150‑1 (MQ=255)
                cAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTACACGAGAGTCGAACTTGCCGGCATCTACCAGCAGGAGAAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTc                                                                                      >  1:444710/1‑150 (MQ=255)
                  gaccgaccTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTcgc                                                                                    <  1:464462/150‑1 (MQ=255)
                       accTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTg                                                                               >  2:218044/1‑150 (MQ=255)
                          tGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGt                                                                            >  1:722878/1‑150 (MQ=255)
                          tGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGt                                                                            <  2:722878/150‑1 (MQ=255)
                          tGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGACAGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGt                                                                            <  1:262288/150‑1 (MQ=255)
                                        aaGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGaa                                                              >  2:642544/1‑150 (MQ=255)
                                             tCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCTTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTCTTCCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCg                                                         <  2:151613/150‑1 (MQ=255)
                                             tCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCg                                                         <  1:642544/150‑1 (MQ=255)
                                              cAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCgg                                                        <  1:218044/150‑1 (MQ=255)
                                               aGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCggg                                                       <  2:112011/150‑1 (MQ=255)
                                                     tCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTg                                                 <  2:444710/150‑1 (MQ=255)
                                                     tCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGGACACCAGGAAGGTCGGGTAGGTg                                                 >  1:284254/1‑150 (MQ=255)
                                                              tCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTGAACTCCGGt                                        <  2:75041/150‑1 (MQ=255)
                                                                gAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTGAACTCCGGTGc                                      <  2:589205/150‑1 (MQ=255)
                                                                     ttCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAg                                                                                                  <  2:97180/85‑1 (MQ=255)
                                                                     ttCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAg                                                                                                  >  1:97180/1‑85 (MQ=255)
                                                                     ntCAGATAACCGGGTCCCAGCGTTGTGANCAGTCGCTGGCGCCGGCANCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTGAACTCCGGTGCCTGCg                                 >  1:109654/2‑150 (MQ=255)
                                                                        aGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAg                                                    >  1:190736/1‑128 (MQ=255)
                                                                        aGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAg                                                    <  2:190736/128‑1 (MQ=255)
                                                                                     ccAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCa                                                                                                   >  2:48717/1‑68 (MQ=255)
                                                                                     ccAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCa                                                                                                   <  1:48717/68‑1 (MQ=255)
                                                                                                    atCGCTGGCGCCGGCAGCTAACAGCCGGAGCAGGATCCCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTGAACTCCGGTGCCTGCGGTACCCGCTCGAACACGCCGCTGTCCAGATa  <  1:393800/149‑1 (MQ=255)
                                                                                                                   |                                                                                                                                      
GGTGGGGGTGATCACCCAGACCGACCTGATCGCCGGCCTCAAGCGTCAGTTGCTCAGCAAGATCGAGGCTTCAGATAACCGGGTCCCAGCGTTGTGACCAGTCGCTGGCGCCGGCAGCTACCAGCCGGCGCAGGATCGCGAAGGCCTGTTGCAGGGCTTCGCTGTCGCGTTTGCGTGGGTACACCAGGAAGGTCGGGTAGGTGAACTCCGGTGCCTGCGGTACCCGCTCGAACACGCCGCTGTCCAGATA  >  NC_002947/2680361‑2680610

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