Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,211,909 A→G 7.4% N14S (AAT→AGT)  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,9090AG 7.4% 149.7 / 2.7 54N14S (AAT→AGT) PP_1955cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base G (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

TGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGGAAAC  >  NC_002947/2211764‑2212056
                                                                                                                                                  |                                                                                                                                                   
tGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATcc                                                                                                                                                  >  7:320914/1‑149 (MQ=255)
  cGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATc                                                                                                                                                   <  2:112870/146‑1 (MQ=255)
  cGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGa                                                                                                                                                >  5:179905/1‑149 (MQ=255)
     cACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGATTGAAGTCCGAtgg                                                                                                                                             <  1:21328/149‑2 (MQ=255)
      aCTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgag                                                                                                                                            <  2:246398/149‑1 (MQ=255)
            aGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGATTGAAGTCCGAtggcgct                                                                                                                                         <  2:175923/146‑6 (MQ=255)
            aGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCa                                                                                                                                      <  3:118780/149‑1 (MQ=255)
                 aaTATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                               <  1:16566/135‑1 (MQ=255)
                 aaTATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                               >  2:16566/1‑135 (MQ=255)
                     taGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCg                                                                                                                             <  6:179905/149‑1 (MQ=255)
                        aTTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGa                                                                                                                          >  4:278946/1‑149 (MQ=255)
                             cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                               <  2:243070/123‑1 (MQ=255)
                             cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCAGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                               >  1:243070/1‑123 (MQ=255)
                              tCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGAGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATCAGTGAGCAGAGTTTTCGTGATCCAg                                                                                                                     >  7:42051/1‑148 (MQ=255)
                                aGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGAGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATCAGTGAGCAGAGTTTTCGTGATCCAg                                                                                                                     <  8:42051/146‑1 (MQ=255)
                                  cGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCACAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGAACCAGCCTc                                                                                                                 >  8:226730/1‑148 (MQ=255)
                                        gTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCg                                                                                                                             >  2:111076/1‑130 (MQ=255)
                                        gTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGCAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCg                                                                                                                             <  1:111076/130‑1 (MQ=255)
                                             gATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGACCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCg                                                                                                     >  5:148173/1‑149 (MQ=255)
                                              aTACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGt                                                                                                    >  5:158921/1‑149 (MQ=255)
                                              aTACTTACTTGGTGGGTAGCCGCGTGCATGAGGGAAATAATTGTTCGGAGAATAACAAAAATGAAAATACTTGAACGTCCTCAACCGCCATCGCAATTCAATCCGATGAGTGAGAAGAGTTTTCATGATCCAGAATAAAACTGCCACCg                                                                                                     >  1:210661/1‑149 (MQ=255)
                                                    acttGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgcgc                                                                                                <  2:156986/147‑1 (MQ=255)
                                                          tggGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCgaaga                                                                                         >  7:201431/1‑148 (MQ=255)
                                                            tgTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGaaa                                                                                       <  7:233471/147‑1 (MQ=255)
                                                             gTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcc                                                                                     <  3:278946/149‑1 (MQ=255)
                                                                gCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCt                                                                                  <  6:208733/149‑1 (MQ=255)
                                                                  cgcgTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGt                                                                                >  3:226139/1‑149 (MQ=255)
                                                                           gA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGaaa                                                                                       <  2:290233/133‑1 (MQ=255)
                                                                           gA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGaaa                                                                                       >  1:290233/1‑133 (MQ=255)
                                                                            a‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCt                                                                        <  4:188279/147‑1 (MQ=255)
                                                                              ggAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATg                                                                     >  2:330637/1‑149 (MQ=255)
                                                                                    aatTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATgcgc                                                                  <  4:117595/146‑1 (MQ=255)
                                                                                    aatTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATgcgc                                                                  >  3:117595/1‑146 (MQ=255)
                                                                                     atTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgcgcgc                                                                                              <  4:279560/117‑1 (MQ=255)
                                                                                     atTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgcgcgc                                                                                              >  3:279560/1‑117 (MQ=255)
                                                                                      ttGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGTGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccgc                                                               >  1:245283/1‑147 (MQ=255)
                                                                                             gagaATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCgtgt                                                      >  7:61422/1‑149 (MQ=255)
                                                                                             gagaATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCgtgt                                                      <  1:330637/149‑1 (MQ=255)
                                                                                                   aaCAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATgg                                                >  5:131220/1‑149 (MQ=255)
                                                                                                      aaaaaTGGAAATACTTGATCGTCCTCAATCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTaa                                             <  6:148173/149‑1 (MQ=255)
                                                                                                       aaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAAc                                            >  3:253250/1‑149 (MQ=255)
                                                                                                          aTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAAcgcg                                         >  4:144661/1‑149 (MQ=255)
                                                                                                            ggAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                       >  3:327392/1‑149 (MQ=255)
                                                                                                            ggAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                       >  1:226626/1‑149 (MQ=255)
                                                                                                                 tACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaa                                  <  6:158921/149‑1 (MQ=255)
                                                                                                                       aTCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGcc                            >  5:47108/1‑149 (MQ=255)
                                                                                                                        tCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGATTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCCAGAAGATGcc                            >  3:230460/1‑148 (MQ=255)
                                                                                                                                aaGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGcgcg                     <  4:226139/147‑1 (MQ=255)
                                                                                                                                 aGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAGCCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGttt                  >  5:331277/1‑149 (MQ=255)
                                                                                                                                     tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt                 >  8:173167/2‑146 (MQ=255)
                                                                                                                                     tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAg              >  1:173635/2‑149 (MQ=255)
                                                                                                                                        tCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGGTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGa           <  2:226626/149‑1 (MQ=255)
                                                                                                                                        tCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGc             >  3:241891/1‑147 (MQ=255)
                                                                                                                                            aCTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGGCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTgg       <  2:245283/149‑1 (MQ=255)
                                                                                                                                                 aaTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAGCCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGGAAAc  <  6:331277/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGGAAAC  >  NC_002947/2211764‑2212056

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