Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,211,906 T→G 7.7% 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 7.7% 155.0 / 4.6 52F13C (TTC→TGC) ‡PP_1955cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base G (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GGCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTG  >  NC_002947/2211758‑2212050
                                                                                                                                                      |                                                                                                                                                 
gcgctctTCGGCAATCTAACTTAAATATAGATTTCCTCAGCGTAATGTAGAGATACTTTCTTGGTGGGATGCCGCGTTCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGTACGTCCTCAAGCGCCATCGGACtt                                                                                                                                                   <  2:98316/146‑1 (MQ=255)
    tctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAAt                                                                                                                                               >  4:47206/1‑149 (MQ=255)
    tctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAAt                                                                                                                                               >  5:78691/1‑149 (MQ=255)
         ggCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTCGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAt                                                                                                                                          <  3:112932/149‑1 (MQ=255)
           cACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGATTGAAGTCCGAtgg                                                                                                                                        <  2:304196/149‑2 (MQ=255)
             ctctAAGTTAAATCTAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgagt                                                                                                                                      <  6:222296/149‑1 (MQ=255)
              tctAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACGTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCGGAtgagtg                                                                                                                                     >  8:151248/1‑149 (MQ=255)
                tAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGTCGCCATCGGACTTCAATCCGAtgagt                                                                                                                                      <  1:4170/147‑1 (MQ=255)
                  aGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCa                                                                                                                                 <  1:332983/149‑1 (MQ=255)
                           taGATTTCCTCAGCGTAATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTc                                                                                                                         <  4:290955/148‑1 (MQ=255)
                                   cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGGCCTCAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGGGCAGAGTTTTCGTGATc                                                                                                                   >  1:256406/1‑147 (MQ=255)
                                   cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCa                                                                                                                 <  4:26231/148‑1 (MQ=255)
                                   cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAg                                                                                                                >  5:41988/1‑149 (MQ=255)
                                    tCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGc                                                                                                               >  1:365241/1‑149 (MQ=255)
                                     cAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGcc                                                                                                              >  3:286261/1‑149 (MQ=255)
                                               tAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGc                                                                                                     <  3:47206/148‑1 (MQ=255)
                                                 gTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAg                                                                                                  >  7:199052/1‑149 (MQ=255)
                                                      acttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgc                                                                                             >  5:211935/1‑149 (MQ=255)
                                                        ttagttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTAGATCGTCCTCAAGCGCCATCGGACTTCAATCAGATGAGTGAGAAGAGTTTTCGTGATCATGCATCAAACTGCCAGCGTga                                                                                             >  5:84358/1‑146 (MQ=255)
                                                            ttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGcgcgcg                                                                                        >  1:211159/1‑148 (MQ=255)
                                                              ggtggGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCgaa                                                                                      <  4:86868/148‑1 (MQ=255)
                                                              ggtggGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCgaa                                                                                      <  2:365241/148‑1 (MQ=255)
                                                                   gTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgag                                                                                                                                       >  8:117112/1‑94 (MQ=255)
                                                                   gTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgag                                                                                                                                       <  7:117112/94‑1 (MQ=255)
                                                                   gTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcc                                                                                <  4:286261/149‑1 (MQ=255)
                                                                      gCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcccc                                                                              <  6:78691/148‑1 (MQ=255)
                                                                                   aGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTAt                                                                <  6:171645/149‑1 (MQ=255)
                                                                                    ggAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTAt                                                                <  7:293674/148‑1 (MQ=255)
                                                                                     gAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATgc                                                              >  7:190888/1‑149 (MQ=255)
                                                                                           atTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccgc                                                         >  5:276498/1‑148 (MQ=255)
                                                                                            ttGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccgc                                                         >  3:248904/1‑147 (MQ=255)
                                                                                               ttCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGc                                                    <  2:211159/149‑1 (MQ=255)
                                                                                                   gagaATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCgtgt                                                <  6:41988/149‑1 (MQ=255)
                                                                                                      aaTAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCCGCCTCAATCTGCCAGCGTGCGCGCGAAGAGACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGa                                             >  2:370045/1‑149 (MQ=255)
                                                                                                          aCAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGGGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATgg                                          <  7:151248/148‑1 (MQ=255)
                                                                                                                aTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAAcgcg                                   >  1:204578/1‑149 (MQ=255)
                                                                                                                    aaaTACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGAGTGCGCGCGAAGAAACCCCTGTCTTTTTCTAt                                                                >  4:304077/1‑116 (MQ=255)
                                                                                                                    aaaTACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGAGTGCGCGCGAAGAAACCCCTGTCTTTTTCTAt                                                                <  3:304077/116‑1 (MQ=255)
                                                                                                                       tACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaa                            <  4:248904/149‑1 (MQ=255)
                                                                                                                             aTCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGcc                      <  8:199052/149‑1 (MQ=255)
                                                                                                                              tCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCg                     >  5:49133/1‑149 (MQ=255)
                                                                                                                                     cAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCCTTTTTCTATGCGCCGCTAGGCgtgt                                                <  1:5181/116‑1 (MQ=255)
                                                                                                                                     cAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCgtgt                                                >  2:5181/1‑115 (MQ=255)
                                                                                                                                        gcgcCATCTGCCTTCAATCCGATGATTGAGCAGATTTTTGGTGATCCATCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCTATAGGCGTGTGGATGGTAACGCGTCTAGAAGATGCCTAGCGCGtttt           <  6:84358/149‑1 (MQ=255)
                                                                                                                                        gcgcCATCGGACTTCAATCCGATGAGTGTGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt           <  2:256406/149‑1 (MQ=255)
                                                                                                                                        gcgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt           <  6:49133/149‑1 (MQ=255)
                                                                                                                                        gcgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt           >  5:2324/1‑149 (MQ=255)
                                                                                                                                            tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTg          >  6:226524/2‑147 (MQ=255)
                                                                                                                                            tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTg          <  5:226524/146‑1 (MQ=255)
                                                                                                                                             cATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGc       >  6:180253/1‑149 (MQ=255)
                                                                                                                                             cATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGc       >  8:48463/1‑149 (MQ=255)
                                                                                                                                                cGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAg    >  6:16781/1‑149 (MQ=255)
                                                                                                                                                  gACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTg  >  4:152040/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GGCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTG  >  NC_002947/2211758‑2212050

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: