Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,776,798 T→G 8.1% T353T (ACT→ACG)  betA‑II → choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,776,7980TG 8.1% 153.0 / 3.7 57T353T (ACT→ACG) betA‑IIcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base G (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAAG  >  NC_002947/5776654‑5776945
                                                                                                                                                |                                                                                                                                                   
ccTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAATAGATTTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACGCGTcc                                                                                                                                                >  1:1060472/1‑150 (MQ=255)
   gCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACGAGTCCGGa                                                                                                                                             <  1:561096/150‑1 (MQ=255)
              gCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGc                                                                                                                                  <  2:138226/150‑1 (MQ=255)
              gCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGc                                                                                                                                  <  2:1574811/150‑1 (MQ=255)
                  ccAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGAGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGaa                                                                                                                              <  1:197136/150‑1 (MQ=255)
                     gCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACAt                                                                                                                           >  1:1405022/1‑150 (MQ=255)
                        ggTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCATCTAGCTCTAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGACCCTCCa                                                                                                                        >  1:1755628/1‑150 (MQ=255)
                        ggTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCa                                                                                                                        <  2:664645/150‑1 (MQ=255)
                           ttCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTa                                                                                                                     <  2:254658/150‑1 (MQ=255)
                               cTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTAccgc                                                                                                                 >  1:860469/1‑148 (MQ=255)
                               cTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCCAACATCCAGTACCAc                                                                                                                 <  2:465610/150‑1 (MQ=255)
                                tataCCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACt                                                                                                                >  2:1537918/1‑150 (MQ=255)
                                 ataCCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACtt                                                                                                               >  1:93907/1‑150 (MQ=255)
                                   aCCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTcc                                                                                                             >  1:243104/1‑150 (MQ=255)
                                     ccGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTg                                                                                                           >  2:246559/1‑150 (MQ=255)
                                           ctgctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGtt                                                                                                     <  1:413831/150‑1 (MQ=255)
                                              ctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       <  2:236781/113‑1 (MQ=255)
                                              ctgTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCaa                                                                                                                                       >  1:236781/1‑113 (MQ=255)
                                                      gAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTa                                                                                          >  1:566042/1‑150 (MQ=255)
                                                        aCCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACa                                                                                        >  1:444835/1‑150 (MQ=255)
                                                               ggCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTc                                                                                 >  2:1475548/1‑150 (MQ=255)
                                                                    tCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATg                                                                            <  1:1537918/150‑1 (MQ=255)
                                                                            gAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTGCAACGGCTCCAATGGCGTGAAg                                                                    <  2:1670767/150‑1 (MQ=255)
                                                                                  aTGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACaa                                                              >  1:46002/1‑149 (MQ=255)
                                                                                   tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       >  1:1127524/1‑92 (MQ=255)
                                                                                   tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       <  2:1127524/92‑1 (MQ=255)
                                                                                   tGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGGTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             >  2:548321/1‑150 (MQ=255)
                                                                                       cAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCtt                                                         >  1:1659872/1‑150 (MQ=255)
                                                                                        aaTGACACCGTGATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTc                                                        >  2:1368450/1‑150 (MQ=255)
                                                                                          tGGCGCCGATACCGGTGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATg                                                             <  1:201568/143‑1 (MQ=255)
                                                                                           ggCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAg                                                     <  2:247359/150‑1 (MQ=255)
                                                                                            gCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAgg                                                    <  2:1774737/150‑1 (MQ=255)
                                                                                              aCCGGTATCGGCGCCAGCCACCAGTTCGAGGCTGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgcg                                                  >  2:419490/1‑150 (MQ=255)
                                                                                              aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       <  2:27434/81‑1 (MQ=255)
                                                                                              aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAg                                                                                                                       >  1:27434/1‑81 (MQ=255)
                                                                                              aCCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGGTTACCTACAACGGCTCCAATGGCGCGACGGCACACGGCTTCCAcgcg                                                  >  2:208626/1‑150 (MQ=255)
                                                                                               ccGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgcgc                                                 >  2:1634775/1‑150 (MQ=255)
                                                                                                  gTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGcaca                                              >  2:228837/1‑150 (MQ=255)
                                                                                                     tCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATgg                                           <  1:1600864/150‑1 (MQ=255)
                                                                                                      gggCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTACCTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATggg                                          <  1:1301387/149‑1 (MQ=255)
                                                                                                        gcgcCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCt                                        <  1:1924023/150‑1 (MQ=255)
                                                                                                            cAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCAt                                    <  2:1182026/150‑1 (MQ=255)
                                                                                                            cAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCAt                                    <  1:1475548/150‑1 (MQ=255)
                                                                                                                aaCCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATgcgc                                <  2:983473/150‑1 (MQ=255)
                                                                                                                aaCCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATgcgc                                >  2:1487672/1‑150 (MQ=255)
                                                                                                                  ccAGTTCGAGGCCGGGGGCTCCATCCAGCCCCGTACAGGCTTCAAGTGGCCGAACCCGCGGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTc                              <  1:1454281/150‑1 (MQ=255)
                                                                                                                   cAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCg                             <  1:1634775/150‑1 (MQ=255)
                                                                                                                     ggtCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGcc                           <  2:305213/148‑1 (MQ=255)
                                                                                                                        cGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGc                        <  2:1405022/150‑1 (MQ=255)
                                                                                                                         gAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGGTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGcc                       >  1:1101159/1‑150 (MQ=255)
                                                                                                                               ggggggTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGc                 <  2:243104/150‑1 (MQ=255)
                                                                                                                                 ggggTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGCCg               >  2:35364/1‑150 (MQ=255)
                                                                                                                                        ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc                                                   <  1:309053/106‑1 (MQ=255)
                                                                                                                                        ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGgc                                                   >  2:309053/2‑107 (MQ=255)
                                                                                                                                        ctCCGGACGAGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAg        >  1:164366/2‑150 (MQ=255)
                                                                                                                                         tCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAgg       >  1:1319643/1‑150 (MQ=255)
                                                                                                                                              aCTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAAg  <  2:1659872/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCTGCAGTACGCCTGCACCCAGCCGGTTTCGCTATACCCGTCGCTGCTGTGGTGGAACCAGCCGGCGATCGGTGCCGAGTGGATGTTCAATGGCACCGGTATCGGCGCCAGCAACCAGTTCGAGGCCGGGGGGTTCATCCGGACTCGTCCGGAGTTCAAGTGGCCGAACATCCAGTACCACTTCCTGCCGGTTGCGATTAACTACAACGGCTCCAATGGCGTGAAGGAACATGGCTTCCAGGCGCACATGGGCTCCATGCGCTCGCCTGCCCGTGGCCGCATCCAGGTCAAG  >  NC_002947/5776654‑5776945

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