Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,773,577 G→A A455V (GCC→GTC)  betT‑I ← choline transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,773,5770GA98.0% 169.8 / ‑4.1 51A455V (GCC→GTC) betT‑Icholine transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (33/17);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CAGCACCACCGAGAACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGGCATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCGGGAAGT  >  NC_002947/5773447‑5773709
                                                                                                                                  |                                                                                                                                    
cagcaccaCCGAGAACGGCAGGCCCGCCCCCACGACCATGGTCTGCATCGCTTCGAAGATACCGGCAAACAGCAGGCCGATGGACACCAGGGAGATCAACACGGACCAGAACACCACCATCCAGAGCGGGACATCt                                                                                                                                 >  2:228168/1‑136 (MQ=255)
cagcaccaCCGAGAACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCAACCACTGCGGGACATCt                                                                                                                                 >  1:194774/1‑136 (MQ=255)
cagcaccaCCGAGAACGGCAGGCCCGCCAGCACGACAATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCt                                                                                                                                 >  1:1058304/1‑136 (MQ=255)
  gcaccaCCGAGAACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTtc                                                                                                                               >  1:936790/1‑136 (MQ=255)
     ccaccGAGAACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTtcgtc                                                                                                                            >  1:1359853/1‑136 (MQ=255)
          gagaACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTtcg                                                                                                                              >  2:222997/1‑129 (MQ=255)
          gagaACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTtcg                                                                                                                              <  1:222997/129‑1 (MQ=255)
          gagaACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCt                                                                                                                       <  2:936790/136‑1 (MQ=255)
          gagaACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCA‑‑‑‑‑‑‑GGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGGCATCTTCGTCGACCTTgccgcc                                                                                                                <  1:268175/136‑1 (MQ=255)
           agaACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCtt                                                                                                                      >  2:1317140/1‑136 (MQ=255)
             aaCGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTgc                                                                                                                    >  1:858949/1‑136 (MQ=255)
             aaCGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTgc                                                                                                                    <  2:1281756/136‑1 (MQ=255)
             aaCGGCAGCCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTgc                                                                                                                    <  2:913700/136‑1 (MQ=255)
              aCGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCCCCATCCAGTGCGGGACATCTTCGTCGACCTTgcc                                                                                                                   >  2:631981/1‑136 (MQ=255)
               cGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTgccg                                                                                                                  <  1:325637/136‑1 (MQ=255)
                  cAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACGGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCt                                                                                                               <  1:899823/136‑1 (MQ=255)
                   aGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCtg                                                                                                              >  1:568526/1‑136 (MQ=255)
                              cACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAgg                                                                                                   <  1:1362862/136‑1 (MQ=255)
                              cACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAgg                                                                                                   <  2:1309069/136‑1 (MQ=255)
                              cACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAgg                                                                                                   <  1:1289593/136‑1 (MQ=255)
                                       ggTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCcatc                                                                                          >  2:1215975/1‑136 (MQ=255)
                                                ggCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACcgcg                                                                                 >  1:1197366/1‑136 (MQ=255)
                                                ggCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACcgcg                                                                                 >  1:419611/1‑136 (MQ=255)
                                                ggCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACcgcg                                                                                 >  1:984915/1‑136 (MQ=255)
                                                ggCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACcgcg                                                                                 >  1:1160927/1‑136 (MQ=255)
                                                ggCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCAAGACAGGTTGGCCATCATCACcgcc                                                                                 >  1:990198/1‑135 (MQ=255)
                                                        aGTTACCGGCAAACAGCAGCCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAggc                                                                         <  1:183916/136‑3 (MQ=255)
                                                               ggCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTcggccggc                                                                  >  1:1000083/1‑136 (MQ=255)
                                                               ggCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTcggccggc                                                                  <  1:1365169/136‑1 (MQ=255)
                                                                       gcagGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACAACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGaa                                                          >  2:708189/1‑136 (MQ=255)
                                                                          ggCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAg                                                       <  1:84298/136‑1 (MQ=255)
                                                                           gCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGc                                                      >  1:675396/1‑136 (MQ=255)
                                                                                  gTCACCAGGGTGATCACCACTGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGTCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAAcc                                               <  1:228168/136‑1 (MQ=255)
                                                                                        aGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCcgacga                                         >  1:726194/1‑136 (MQ=255)
                                                                                             gATCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACt                                    >  1:763351/1‑136 (MQ=255)
                                                                                              aTCACCACGTACCAGAACACCACCATCCTGTGCGGGACATCTTGGTCGACCTTGCCGCCTGTGCACGACTGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTg                                   <  1:646368/136‑1 (MQ=255)
                                                                                              aTCACCACGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTg                                   <  2:763351/136‑1 (MQ=255)
                                                                                                 accacGGACCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGccc                                >  2:302927/1‑135 (MQ=255)
                                                                                                        aCCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAAt                         >  1:134159/1‑136 (MQ=255)
                                                                                                        aCCAGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAAt                         >  1:720863/1‑136 (MQ=255)
                                                                                                           aGAACACCACCATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGATGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGAc                      >  1:142191/1‑136 (MQ=255)
                                                                                                                accaccATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAg                                                             <  1:106749/92‑1 (MQ=255)
                                                                                                                accaccATCCAGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAg                                                             >  2:106749/1‑92 (MQ=255)
                                                                                                                          aGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCgg       >  1:614482/1‑136 (MQ=255)
                                                                                                                          aGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCgg       >  1:958060/1‑136 (MQ=255)
                                                                                                                          aGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCgg       >  1:712924/1‑136 (MQ=255)
                                                                                                                          aGTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCAGGCCTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCgg       >  1:1280684/1‑136 (MQ=255)
                                                                                                                          aGTGCGGGACATCTTCGTCGACCTTACCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCgg       >  1:1174015/1‑136 (MQ=255)
                                                                                                                          aGTGCGGGACAACTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCAGCGTGAGGAACAGCACGAACCCGACGAACACAGCCACACCAATGACCATCTTGGCTTCCgc       >  2:159936/1‑135 (MQ=255)
                                                                                                                          aGTGCCGGACATCTTAGTCGACATTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGAGAGCAAAAGCACGAACCCGCCGAACACTGCCACACCAATGACGCTCTTGGCTGCCgg       >  2:139921/1‑136 (MQ=255)
                                                                                                                           gTGCGGGACATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCggg      <  1:1079611/136‑1 (MQ=255)
                                                                                                                               gggaCATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCGGGAAgt  >  2:1099149/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                    
CAGCACCACCGAGAACGGCAGGCCCGCCAGCACGACCATGGTCTGCATGGCTTCGAAGTTACCGGCAAACAGCAGGCCGATGGTCACCAGGGTGATCACCACGGACCAGAACACCACCATCCAGTGCGGGGCATCTTCGTCGACCTTGCCGCCTGTGCACGACAGGTTGGCCATCATCACCGCGCCGGAGTCGGCCGGCGTGAGGAACAGCACGAACCCGACGAACACTGCCACACCAATGACGATCTTGGCTGCCGGGAAGT  >  NC_002947/5773447‑5773709

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