Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,211,906 T→G 9.5% F13C (TTC→TGC) ‡ 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,9060TG 9.5% 155.1 / 6.4 52F13C (TTC→TGC) ‡PP_1955cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base T (22/25);  new base G (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.90e-01

TCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGG  >  NC_002947/2211762‑2212052
                                                                                                                                                |                                                                                                                                                  
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                                                                                                                      >  3:276743/1‑98 (MQ=255)
                                                                             gAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGa                                                                                                                      <  4:276743/98‑1 (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                           <  3:241595/131‑1 (MQ=255)
                                                                                                                                  tcactcATCGGATTGAAGTCCGATGGGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGAt                           >  4:241595/6‑136 (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)
                                                                                                                                               ttCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTggg  >  3:13657/1‑148 (MQ=255)
                                                                                                                                               ttCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTggg  >  8:104406/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGG  >  NC_002947/2211762‑2212052

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