Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,431,526 T→C 7.4% A64A (GCT→GCC PP_3926 → short‑chain dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,431,5260TC 7.4% 155.7 / 2.6 54A64A (GCT→GCCPP_3926short‑chain dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (19/31);  new base C (2/2);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGCAGCCACCAACCCGCAATTCTGCAATGTGCTGGACACCGACCCAGGG  >  NC_002947/4431381‑4431667
                                                                                                                                                   |                                                                                                                                             
cttgcGGCGCCTGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCCGCAAGGGGCCCTGGTGATCGTTTCCAGCCTCAAGCTTGCCGGGGGCCAGCAGTTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCtgt                                                                                                                                             <  5:5264/144‑1 (MQ=255)
       cgcCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACAGCCGTGGCCttgc                                                                                                                                     <  3:88610/149‑4 (MQ=255)
       cgcCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTCACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACAGCCGTGGCCttgc                                                                                                                                     >  4:88610/1‑146 (MQ=255)
        gcCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCca                                                                                                                                    >  2:277503/1‑149 (MQ=255)
            gCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACAtt                                                                                                                                >  1:159090/1‑149 (MQ=255)
             ccGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTg                                                                                                                               <  1:277503/149‑1 (MQ=255)
             ccGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGGGCCAGCCGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTg                                                                                                                               <  5:248163/149‑1 (MQ=255)
               gTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGt                                                                                                                             >  1:361430/1‑149 (MQ=255)
                tGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAATCTTGACGGGTGCCAGCAGGTGGCCGACGCCACTCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGt                                                                                                                             <  1:82157/149‑1 (MQ=255)
                tGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGt                                                                                                                             <  1:77541/149‑1 (MQ=255)
                tGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGt                                                                                                                             >  6:2717/1‑148 (MQ=255)
                tGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTg                                                                                                                            <  8:187452/149‑1 (MQ=255)
                    ctCGGCGAGGCCATCGCCCTCTTGCTCGCGCAGCAAGGGGCGCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCt                                                                                                                        <  8:60319/148‑1 (MQ=255)
                    aTCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCt                                                                                                                        <  6:38850/149‑1 (MQ=255)
                    aTCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCt                                                                                                                        >  4:32429/1‑149 (MQ=255)
                    aTCGGCGAGGCCATCGCCCACTTCCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCCCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACAGCCGTGGCCtt                                                                                                                                       <  7:296866/134‑2 (MQ=255)
                     ggggcttggCCATCGCCCACTTGCTCGCGCAGAAAGGGGCCCATGGTAGCGTTTCCAGCCGCACGCTTGACGGGGGCCAGCAGGGGCCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTg                                                                                                                       <  2:159090/142‑1 (MQ=255)
                       ggCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTgg                                                                                                                      <  6:225118/148‑1 (MQ=255)
                            ggCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGa                                                                                                                <  1:170931/149‑1 (MQ=255)
                             gCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGAt                                                                                                               >  8:17576/1‑149 (MQ=255)
                               cATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGTGGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTc                                                                                                             <  3:86061/149‑1 (MQ=255)
                               cATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGAtt                                                                                                              >  1:139813/1‑148 (MQ=255)
                                   gCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTcagca                                                                                                         <  1:25048/149‑1 (MQ=255)
                                    cccACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTcagcag                                                                                                        <  5:12448/149‑1 (MQ=255)
                                    cccACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTcagcag                                                                                                        >  6:12448/1‑149 (MQ=255)
                                     ccACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAgg                                                                                                       <  3:32429/149‑1 (MQ=255)
                                     ccACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAgg                                                                                                       >  4:364055/1‑149 (MQ=255)
                                         ttGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGtt                                                                                                   >  7:158735/1‑149 (MQ=255)
                                           ccgcgcgCAGCATGGGGCCCATGTGATCTTTTCCCCCCGCACGCTTGCGGGGGGCCCGCCGGTGCCCGACGCCATCTTTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTc                                                                                                  <  2:361430/145‑1 (MQ=255)
                                             tCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGa                                                                                                                >  5:179793/1‑132 (MQ=255)
                                             tCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGa                                                                                                                <  6:179793/132‑1 (MQ=255)
                                               gcgcAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGcc                                                                                               <  3:135879/147‑1 (MQ=255)
                                               gcgcAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCgg                                                                                             >  2:176582/1‑149 (MQ=255)
                                                 gcagcaAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCa                                                                                           >  5:82641/1‑149 (MQ=255)
                                                      aaGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTcgc                                                                                      <  3:239586/149‑1 (MQ=255)
                                                       aGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTcgcg                                                                                     <  6:179826/149‑1 (MQ=255)
                                                        ggggCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAg                                                                                                                          >  8:109508/1‑111 (MQ=255)
                                                        ggggCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAg                                                                                                                          <  7:109508/111‑1 (MQ=255)
                                                         gggCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGaa                                                                                   >  6:159770/1‑149 (MQ=255)
                                                          ggCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAAc                                                                                  <  3:207353/149‑1 (MQ=255)
                                                                 gtgATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGtt                                                                              >  7:135787/1‑146 (MQ=255)
                                                                 gtgATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGtt                                                                              <  8:135787/146‑1 (MQ=255)
                                                                  tgATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCggg                                                                          <  8:158735/149‑1 (MQ=255)
                                                                       gTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGGGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGAc                                                                      <  3:154281/148‑1 (MQ=255)
                                                                             aGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGt                                                               >  8:176949/1‑149 (MQ=255)
                                                                                 gCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTg                                                           <  3:173291/149‑1 (MQ=255)
                                                                                 gCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTg                                                           >  4:173291/1‑149 (MQ=255)
                                                                                 gCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTg                                                           <  3:364055/149‑1 (MQ=255)
                                                                                 gCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGCTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTg                                                           <  2:139813/149‑1 (MQ=255)
                                                                                   aaGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGt                                                         <  2:131881/149‑1 (MQ=255)
                                                                                               tGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGCAGc                                             <  5:2717/149‑1 (MQ=255)
                                                                                                            ccGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGCAGCCACCAACCCGCaa                                <  3:50864/149‑1 (MQ=255)
                                                                                                                   catcatTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGAAGCCACCAACCCGCAATTCTGc                          >  6:307139/1‑148 (MQ=255)
                                                                                                                                          ggCCACAGCCGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGCAGCCACCAACCCGCAATTCTGCAATGTGCTGGACAccccagc      >  4:272728/1‑142 (MQ=255)
                                                                                                                                            ccACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGCAGCCACCAACCCGCAATTCTGCAATGTGCTGGACACCGACCCAggg  <  1:176582/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCCCACTTGCTCGCGCAGCAAGGGGCCCATGTGATCGTTTCCAGCCGCAAGCTTGACGGGTGCCAGCAGGTGGCCGACGCCATCATTGCCGCCGGCGGCAAGGCCACGGCTGTGGCCTGCCACATTGGTGAGCTGGAACAGATTCAGCAGGTGTTCGCCGGCATTCGCGAACAGTTCGGGCGACTGGACGTGCTGGTCAACAATGCAGCCACCAACCCGCAATTCTGCAATGTGCTGGACACCGACCCAGGG  >  NC_002947/4431381‑4431667

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