Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,466,113 C→T 8.1% intergenic (‑159/‑145) betT‑III ← / → nhaA‑II choline/carnitine/betaine transporter BetTC/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,466,1130CT 8.1% 137.2 / 2.2 49intergenic (‑159/‑145)betT‑III/nhaA‑IIcholine/carnitine/betaine transporter BetTC/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base T (2/2);  total (22/27)
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.13e-01

CTCCTTGCAAGGCGATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTATG  >  NC_002947/4465965‑4466260
                                                                                                                                                    |                                                                                                                                                   
ctcCTTGCAAGGCGATGAGCACAGCGTGCCTGGGCACGACGACATTTACGAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGcc                                                                                                                                                     >  5:234583/1‑149 (MQ=255)
 tcCTTGCAAGGCGATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGcc                                                                                                                                                     >  2:325402/1‑148 (MQ=255)
   cTTGCAAGGCGATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCa                                                                                                                                                  >  4:52272/1‑149 (MQ=255)
         aGGCGATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAg                                                                                                                                            <  1:21501/149‑1 (MQ=255)
             gATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGaa                                                                                                                                        <  3:52272/149‑1 (MQ=255)
             gATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGaa                                                                                                                                        <  3:327434/149‑1 (MQ=255)
               tGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATg                                                                                                                                      >  3:27144/1‑149 (MQ=255)
                    acaGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGttt                                                                                                                                  <  1:325402/148‑1 (MQ=255)
                    acaGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTg                                                                                                                                 <  8:90371/149‑1 (MQ=255)
                            ccTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCtttt                                                                                                                         <  4:27144/149‑1 (MQ=255)
                            ccTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCttt                                                                                                                          <  7:3429/148‑1 (MQ=255)
                                       cgacATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGc                                                                                                              <  5:344113/149‑1 (MQ=255)
                                         agctttACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTATCGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGa                                                                                                            <  1:296746/146‑1 (MQ=255)
                                         acATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGa                                                                                                            <  6:234583/149‑1 (MQ=255)
                                           aTTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGt                                                                                                          >  8:230424/1‑149 (MQ=255)
                                            tttACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTg                                                                                                         <  3:110303/149‑1 (MQ=255)
                                               aCCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGt                                                                                                       >  4:170225/1‑148 (MQ=255)
                                                    ggCTAGCAGCTGGATTAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCg                                                                                                 <  3:203936/149‑1 (MQ=255)
                                                     gCTACCACCTGGATAAAATCTCGTTTTCCAACCTCCGTGTGCTCGAAAACAAACACAAAACCTAACGTCTGCATGAAAAAATCTCAAACCGTCGCCTCATTGGATTGAATTGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCgg                                                                                                <  7:214405/149‑1 (MQ=255)
                                                                gATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCg                                                                                                 <  4:102679/137‑1 (MQ=255)
                                                                gATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCg                                                                                                 >  3:102679/1‑137 (MQ=255)
                                                                gATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCg                                                                                      >  1:197502/1‑148 (MQ=255)
                                                                   aaaCTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACg                                                                                  >  2:188488/1‑149 (MQ=255)
                                                                     aCTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACggg                                                                                >  2:289050/1‑149 (MQ=255)
                                                                      ctctCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGt                                                                               <  6:20679/149‑1 (MQ=255)
                                                                       tctcGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGc                                                                                                              <  7:311053/117‑1 (MQ=255)
                                                                       tctcGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGc                                                                                                              >  8:311053/1‑117 (MQ=255)
                                                                         tcGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGc                                                                             >  1:335652/1‑148 (MQ=255)
                                                                            tttGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGc                                                                         <  1:188488/149‑1 (MQ=255)
                                                                                ccAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGc                                                                      <  3:170225/148‑1 (MQ=255)
                                                                                  aaGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACg                                                                   <  2:197502/149‑1 (MQ=255)
                                                                                               tCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAgg                                                        >  7:44805/1‑147 (MQ=255)
                                                                                                 gAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAgg                                                    <  1:289050/149‑1 (MQ=255)
                                                                                                 gAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAgg                                                    >  3:269930/1‑149 (MQ=255)
                                                                                                 gAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTCCACGGTGCAGTGAGGCAgg                                                    >  6:57786/1‑149 (MQ=255)
                                                                                                    aagatccAAAAACCCTAACGTCTGCAGGAAAAATTCTCAAACCGACGCCTCATTGGAGTTAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCTCACGGGTGCCAGCTGCACTGTGCAGTGAGGCAGGTTa                                                 <  3:157472/144‑1 (MQ=255)
                                                                                                             aaaCCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAgc                                        >  4:246751/1‑149 (MQ=255)
                                                                                                                  cTAACGTCTGCAGTAAAAAATCTCAAACCGACGCTTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAgg                                                        <  8:159128/128‑1 (MQ=255)
                                                                                                                  cTAACGTCTGCAGGAAAAAATCTCAAACCGACGCTTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAgg                                                        >  3:156759/1‑128 (MQ=255)
                                                                                                                  cTAACGTCTGCAGGAAAAAATCTCAAACCGACGCTTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAgg                                                        <  4:156759/128‑1 (MQ=255)
                                                                                                                  cTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGAAGGTTAAAGAACAGCGCAgg                                   >  3:156506/1‑149 (MQ=255)
                                                                                                                  cTAACGTCTGCAGGAAAAAATCTCAAACAGACGCTTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAgg                                                        >  7:159128/1‑128 (MQ=255)
                                                                                                                             aGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCGCGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCa                        >  1:299747/1‑149 (MQ=255)
                                                                                                                              ggAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCAt                       >  6:115365/1‑149 (MQ=255)
                                                                                                                              ggAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCAt                       >  8:12654/1‑149 (MQ=255)
                                                                                                                                 aaaaaTCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTATGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCg                    >  1:254157/1‑149 (MQ=255)
                                                                                                                                     aTCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAg                <  3:246751/149‑1 (MQ=255)
                                                                                                                                                aCGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGcc      >  1:297443/1‑148 (MQ=255)
                                                                                                                                                 cGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTa    <  7:230424/149‑1 (MQ=255)
                                                                                                                                                 cGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTa    <  2:335652/149‑1 (MQ=255)
                                                                                                                                                  gCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTAt   <  4:269930/149‑1 (MQ=255)
                                                                                                                                                   ccTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTCCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTATg  <  5:57786/149‑1 (MQ=255)
                                                                                                                                                    cTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTATg  >  2:290899/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCCTTGCAAGGCGATGAGCACAGCGTGCCTGGGCACGACGACATTTACCAAGGCTAGCAGCTGGATAAACTCTCGTTTGCCAAGCTCCGTGTGCTCGAAAACAACCACAAACCCTAACGTCTGCAGGAAAAAATCTCAAACCGACGCCTCATTGGAGTGAATGGTTTGGCGCTTTTTGTTCATCTGCGAGTGGTTATTCGGTAGACTACCGCACGGGTGCCAGCTGCACGGTGCAGTGAGGCAGGTTAAAGAACAGCGCAGGTCGGGCCGCCATGCGGGAGCACCTGTGCCTATG  >  NC_002947/4465965‑4466260

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