Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,211,904 C→T 9.6% D12D (GAC→GAT)  PP_1955 → cytochrome P450 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,211,9040CT 9.6% 139.5 / 5.2 52D12D (GAC→GAT) PP_1955cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base T (3/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGGCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAG  >  NC_002947/2211756‑2212044
                                                                                                                                                    |                                                                                                                                            
gtggCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGAc                                                                                                                                              >  4:116409/1‑149 (MQ=255)
      tctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAAt                                                                                                                                        >  5:259718/1‑149 (MQ=255)
           ggCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCTAt                                                                                                                                   <  3:297354/149‑1 (MQ=255)
           ggCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                   <  3:775/149‑1 (MQ=255)
           ggCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                   <  3:297351/149‑1 (MQ=255)
             cACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtga                                                                                                                                 <  2:122473/149‑1 (MQ=255)
                    aGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGATTGAAGTCCGAtggcgctt                                                                                                                            <  7:288163/147‑7 (MQ=255)
                      ttAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCaga                                                                                                                        <  6:259718/149‑1 (MQ=255)
                      ctAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCag                                                                                                                         <  1:297877/147‑1 (MQ=255)
                      ctAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCag                                                                                                                         >  2:297877/2‑148 (MQ=255)
                         aaTCTAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCTTCGGACTTCAATCCGATGAGTGAGCAGAGt                                                                                                                      <  1:222064/148‑1 (MQ=255)
                                 tttCCTCAGCTTCATTTATTTATACTTACTTTGTTGGTAGCCGCTTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGa                                                                                                              <  2:123998/148‑1 (MQ=255)
                                 tttCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGAt                                                                                                             <  3:116409/149‑1 (MQ=255)
                                     cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGtttt                                                                                                                   <  8:75476/139‑1 (MQ=255)
                                     cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGtttt                                                                                                                   >  7:75476/1‑139 (MQ=255)
                                     cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAg                                                                                                         >  1:149251/1‑149 (MQ=255)
                                     cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAACCGCCATCTGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAg                                                                                                         >  5:187874/1‑149 (MQ=255)
                                                         cttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCg                                                                                                                 <  7:103984/121‑1 (MQ=255)
                                                         cttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCg                                                                                                                 >  8:103984/1‑121 (MQ=255)
                                                           tacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgcgcg                                                                                   >  6:131975/1‑149 (MQ=255)
                                                              ttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGcgcgcg                                                                                 >  4:168902/1‑148 (MQ=255)
                                                                    ggTATCCGCGTGCCTGAGGAAATAATTGTTCTGAGAATAACAAATTTGGAAATACTTGATCGTCCTCAAGCGCCTTCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGaaa                                                                           <  6:187874/148‑1 (MQ=255)
                                                                    ggTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAc                                                                          >  5:158036/1‑149 (MQ=255)
                                                                    ggTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAACGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAc                                                                          <  2:223016/149‑1 (MQ=255)
                                                                     gTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAc                                                                          <  3:168902/148‑1 (MQ=255)
                                                                      tAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAccc                                                                        >  4:165489/1‑149 (MQ=255)
                                                                            cgTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCt                                                                  >  6:19216/1‑149 (MQ=255)
                                                                                  tGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTc                                                             >  4:313515/1‑148 (MQ=255)
                                                                                   gAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGa                                                                                                              <  4:276743/98‑1 (MQ=255)
                                                                                   gAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGa                                                                                                              >  3:276743/1‑98 (MQ=255)
                                                                                         ataatTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGcgcc                                                     >  8:196738/1‑149 (MQ=255)
                                                                                                  cGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCg                                             <  6:158036/148‑1 (MQ=255)
                                                                                                   ggAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcccc                                                                       >  6:25438/1‑121 (MQ=255)
                                                                                                   ggAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcccc                                                                       <  5:25438/121‑1 (MQ=255)
                                                                                                         ttaCAAAAATGGAAATTCTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGATGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGa                                       <  5:240567/145‑1 (MQ=255)
                                                                                                         tAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTCAGCAGAGTTTTCGTGATCCAGCCTCAAGCTGCCAGCGTGCGCGCGAAGAAACCCCTGACTTTTTCTATGCGCCGATAGGCGAGTGGa                                       >  6:240567/1‑147 (MQ=255)
                                                                                                           aCAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATgg                                    <  5:19216/148‑1 (MQ=255)
                                                                                                             aaaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAc                                                                          >  5:184294/1‑108 (MQ=255)
                                                                                                             aaaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAc                                                                          <  6:184294/108‑1 (MQ=255)
                                                                                                              aaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAAc                                <  3:313515/149‑1 (MQ=255)
                                                                                                                    gAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                           <  5:131975/148‑1 (MQ=255)
                                                                                                                     aaaTACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCga                         <  2:149251/149‑1 (MQ=255)
                                                                                                                         aCTTGATCGTCCTCAAGCGCCATCGGTCTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCgag                        >  7:105594/1‑146 (MQ=255)
                                                                                                                               tCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGcc                >  6:92979/1‑148 (MQ=255)
                                                                                                                                     tCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGcgcg         >  3:294898/1‑149 (MQ=255)
                                                                                                                                        tcactcATCGGATTGAAGTCCGATGGGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGAt                   >  4:241595/6‑136 (MQ=255)
                                                                                                                                        tcactcATCGGATTGAAGTCCGATGGGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGAt                   <  3:241595/131‑1 (MQ=255)
                                                                                                                                         gcgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt     >  3:233894/1‑149 (MQ=255)
                                                                                                                                          cgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTg    <  7:196738/149‑1 (MQ=255)
                                                                                                                                          actcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt     >  3:233658/4‑148 (MQ=255)
                                                                                                                                            tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAg  >  3:27384/2‑149 (MQ=255)
                                                                                                                                             cATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaaga                    >  2:203536/1‑130 (MQ=255)
                                                                                                                                             cATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaaga                    <  1:203536/130‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GTGGCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAG  >  NC_002947/2211756‑2212044

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