Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,007,282 T→C 12.9% Q41Q (CAA→CAG PP_1790 ← acylneuraminate cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,007,2820TC12.9% 133.8 / 8.8 54Q41Q (CAA→CAGPP_1790acylneuraminate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base C (3/4);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACG  >  NC_002947/2007136‑2007422
                                                                                                                                                  |                                                                                                                                            
tGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGc                                                                                                                                           >  1:1512101/1‑150 (MQ=255)
      gTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaag                                                                                                                                     <  2:718325/150‑2 (MQ=255)
                       aGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaag                                                                                                                                     <  1:1162449/133‑2 (MQ=255)
                       aGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaag                                                                                                                                     >  2:1162449/1‑132 (MQ=255)
                            ggTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAAtt                                                                                                               >  2:1140307/1‑150 (MQ=255)
                              tcttctCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCg                                                                                                             >  1:943786/1‑150 (MQ=255)
                                ttctCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTg                                                                                                           <  1:1140307/150‑1 (MQ=255)
                                ttctCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTg                                                                                                           <  2:1512101/150‑1 (MQ=255)
                                     gtcgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGtt                                                                                                                                             >  2:1203098/1‑111 (MQ=255)
                                     gtcgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGtt                                                                                                                                             <  1:1203098/111‑1 (MQ=255)
                                     gtcgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGgccgc                                                                                                      >  2:1063940/1‑150 (MQ=255)
                                      tcgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCtt                                                                                                                                         <  1:1325098/114‑1 (MQ=255)
                                      tcgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCtt                                                                                                                                         >  2:1325098/1‑114 (MQ=255)
                                        gCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCaa                                                                                                   >  2:1333546/1‑150 (MQ=255)
                                           tATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAAc                                                                                                                                     <  2:747956/113‑1 (MQ=255)
                                           tATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAAc                                                                                                                                     >  1:747956/1‑113 (MQ=255)
                                              gAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaagca                                                                                                                                   >  1:769006/1‑109 (MQ=255)
                                              gAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaagca                                                                                                                                   <  2:769006/112‑4 (MQ=255)
                                                   ttGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTtata                                                                                        <  1:250087/150‑1 (MQ=255)
                                                       ttttATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGttt                                                                                            <  2:1839300/142‑1 (MQ=255)
                                                       ttttATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGttt                                                                                            >  1:1839300/1‑142 (MQ=255)
                                                          tattatTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGc                                                                                                                                           <  1:630569/92‑1 (MQ=255)
                                                          tattatTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGc                                                                                                                                           >  2:630569/1‑92 (MQ=255)
                                                            ttattTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTc                                                                               >  1:552502/1‑150 (MQ=255)
                                                               ttttGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGa                                                                            <  1:1063940/150‑1 (MQ=255)
                                                               ttttGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGa                                                                            >  1:474227/1‑150 (MQ=255)
                                                                     aGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGt                                                                      <  2:943786/150‑1 (MQ=255)
                                                                      gCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGtt                                                                     >  2:927816/1‑150 (MQ=255)
                                                                       cGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAg                                                                                                                      >  1:550886/1‑100 (MQ=255)
                                                                       cGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAg                                                                                                                      <  2:550886/100‑1 (MQ=255)
                                                                          tGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCt                                                                 <  1:648134/150‑1 (MQ=255)
                                                                               cGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCtt                                                            <  1:1255695/150‑1 (MQ=255)
                                                                                   gTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCATCTGCTTGAAGCAGTTGCTTAAACGTGATGTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCg                                                        <  2:942687/150‑1 (MQ=255)
                                                                                       aGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcag                                                    <  2:589205/150‑1 (MQ=255)
                                                                                              tttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCTGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGaa                                             >  2:1632531/1‑150 (MQ=255)
                                                                                                   tcgtTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGc                                                                                                                     >  1:1844019/1‑73 (MQ=255)
                                                                                                   tcgtTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGc                                                                                                                     <  2:1844019/73‑1 (MQ=255)
                                                                                                               aGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCa                            >  1:1183384/1‑150 (MQ=255)
                                                                                                                 tACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCata                          <  2:1183384/150‑1 (MQ=255)
                                                                                                                    gATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATattt                       >  1:1638995/1‑150 (MQ=255)
                                                                                                                      tGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatttca                     <  2:474227/150‑1 (MQ=255)
                                                                                                                       ggAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatttcat                    >  2:851844/1‑150 (MQ=255)
                                                                                                                          aTCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAac                                          >  1:973547/1‑125 (MQ=255)
                                                                                                                          aTCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAac                                          <  2:973547/125‑1 (MQ=255)
                                                                                                                             aaCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCg                                                                              >  2:1922111/1‑86 (MQ=255)
                                                                                                                             aaCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCg                                                                              <  1:1922111/86‑1 (MQ=255)
                                                                                                                             aaCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTTGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatttca                     <  2:792169/143‑1 (MQ=255)
                                                                                                                             aaCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAAGGCAGAAAACACGTCGTTTTTTCATatttca                     >  1:792169/1‑143 (MQ=255)
                                                                                                                              aCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGttt                                 >  1:900423/1‑130 (MQ=255)
                                                                                                                              aCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGttt                                 <  2:900423/130‑1 (MQ=255)
                                                                                                                                         ttGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACg  >  2:977355/1‑150 (MQ=255)
                                                                                                                                         ttGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTAGTCGGAATGGGTTCGCTGCCTTTTTGGCAGGGCAGAAAAGACGTCGTTTTTTCGTATTTCTTTTCAGTGGCCCATAAACg  >  2:276133/1‑150 (MQ=255)
                                                                                                                                         ttGAAGCAACTGCTTCAACTTGATTTCAATCAAGcc                                                                                                                    <  2:1264284/36‑1 (MQ=21)
                                                                                                                                         ttGAAGCAACTGCTTCAACTTGATTTCAATCAAGcc                                                                                                                    >  1:1264284/1‑36 (MQ=21)
                                                                                                                                            aaGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACAcgt                                      >  2:53158/1‑111 (MQ=255)
                                                                                                                                            aaGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACAcgt                                      <  1:53158/111‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACG  >  NC_002947/2007136‑2007422

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: