Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,211,907 C→A 7.4% F13L (TTC→TTA)  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,9070CA 7.4% 171.0 / 3.8 54F13L (TTC→TTA) PP_1955cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base C (28/22);  new base A (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 = 4.68e-01

TGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGG  >  NC_002947/2211764‑2212051
                                                                                                                                                |                                                                                                                                                
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                                                                                         >  3:279560/1‑117 (MQ=255)
                                                                                     atTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgcgcgc                                                                                         <  4:279560/117‑1 (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                                  >  1:226626/1‑149 (MQ=255)
                                                                                                            ggAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                  >  3:327392/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)
                                                                                                                                                |                                                                                                                                                
TGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGA‑GGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGG  >  NC_002947/2211764‑2212051

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