Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,211,909	A→G	7.7%	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_002947	2,211,909	0	A	G	7.7%	143.4 / 2.6	52	N14S (AAT→AGT)	PP_1955	cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base A (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 = 7.71e-01

TCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTG  >  NC_002947/2211762‑2212050                                                                                                                                                      |                                                                                                                                               tctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAAt                                                                                                                                               >  4:47206/1‑149 (MQ=255) tctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAAt                                                                                                                                               >  5:78691/1‑149 (MQ=255)      ggCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTCGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGAt                                                                                                                                          <  3:112932/149‑1 (MQ=255)        cACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGATTGAAGTCCGAtgg                                                                                                                                        <  2:304196/149‑2 (MQ=255)          ctctAAGTTAAATCTAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGAtgagt                                                                                                                                      <  6:222296/149‑1 (MQ=255)           tctAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACGTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCGGAtgagtg                                                                                                                                     >  8:151248/1‑149 (MQ=255)             tAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGTCGCCATCGGACTTCAATCCGAtgagt                                                                                                                                      <  1:4170/147‑1 (MQ=255)               aGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCa                                                                                                                                 <  1:332983/149‑1 (MQ=255)                        taGATTTCCTCAGCGTAATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTc                                                                                                                         <  4:290955/148‑1 (MQ=255)                                cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGGCCTCAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGGGCAGAGTTTTCGTGATc                                                                                                                   >  1:256406/1‑147 (MQ=255)                                cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCa                                                                                                                 <  4:26231/148‑1 (MQ=255)                                cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAg                                                                                                                >  5:41988/1‑149 (MQ=255)                                 tCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGc                                                                                                               >  1:365241/1‑149 (MQ=255)                                  cAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGcc                                                                                                              >  3:286261/1‑149 (MQ=255)                                            tAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGc                                                                                                     <  3:47206/148‑1 (MQ=255)                                              gTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAg                                                                                                  >  7:199052/1‑149 (MQ=255)                                                   acttacttGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTgc                                                                                             >  5:211935/1‑149 (MQ=255)                                                     ttagttGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTAGATCGTCCTCAAG‑CGCCATCGGACTTCAATCAGATGAGTGAGAAGAGTTTTCGTGATCATGCATCAAACTGCCAGCGTga                                                                                             >  5:84358/1‑146 (MQ=255)                                                         ttGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGcgcgcg                                                                                        >  1:211159/1‑148 (MQ=255)                                                           ggtggGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCgaa                                                                                      <  4:86868/148‑1 (MQ=255)                                                           ggtggGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCgaa                                                                                      <  2:365241/148‑1 (MQ=255)                                                                gTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGAtgag                                                                                                                                       <  7:117112/94‑1 (MQ=255)                                                                gTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGAtgag                                                                                                                                       >  8:117112/1‑94 (MQ=255)                                                                gTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcc                                                                                <  4:286261/149‑1 (MQ=255)                                                                   gCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAAcccc                                                                              <  6:78691/148‑1 (MQ=255)                                                                                aGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTAt                                                                <  6:171645/149‑1 (MQ=255)                                                                                 ggAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTAt                                                                <  7:293674/148‑1 (MQ=255)                                                                                  gAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATgc                                                              >  7:190888/1‑149 (MQ=255)                                                                                        atTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGcgccgc                                                         >  5:276498/1‑148 (MQ=255)                                                                                         ttGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGcgccgc                                                         >  3:248904/1‑147 (MQ=255)                                                                                            ttCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGc                                                    <  2:211159/149‑1 (MQ=255)                                                                                                gagaATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCgtgt                                                <  6:41988/149‑1 (MQ=255)                                                                                                   aaTAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCCGCCTCAATCTGCCAGCGTGCGCGCGAAGAGACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGa                                             >  2:370045/1‑149 (MQ=255)                                                                                                       aCAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGGGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATgg                                          <  7:151248/148‑1 (MQ=255)                                                                                                             aTGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAAcgcg                                   >  1:204578/1‑149 (MQ=255)                                                                                                                 aaaTACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGAGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTAt                                                                >  4:304077/1‑116 (MQ=255)                                                                                                                 aaaTACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGAGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTAt                                                                <  3:304077/116‑1 (MQ=255)                                                                                                                    tACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaa                            <  4:248904/149‑1 (MQ=255)                                                                                                                          aTCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGcc                      <  8:199052/149‑1 (MQ=255)                                                                                                                           tCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCg                     >  5:49133/1‑149 (MQ=255)                                                                                                                                  cAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGTCCTTTTTCTATGCGCCGCTAGGCgtgt                                                <  1:5181/116‑1 (MQ=255)                                                                                                                                  cAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCgtgt                                                >  2:5181/1‑115 (MQ=255)                                                                                                                                     g‑cgcCATCTGCCTTCAATCCGATGATTGAGCAGATTTTTGGTGATCCATCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCTATAGGCGTGTGGATGGTAACGCGTCTAGAAGATGCCTAGCGCGtttt           <  6:84358/149‑1 (MQ=255)                                                                                                                                     g‑cgcCATCGGACTTCAATCCGATGAGTGTGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt           <  2:256406/149‑1 (MQ=255)                                                                                                                                     g‑cgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt           <  6:49133/149‑1 (MQ=255)                                                                                                                                     g‑cgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt           >  5:2324/1‑149 (MQ=255)                                                                                                                                         tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTg          >  6:226524/2‑147 (MQ=255)                                                                                                                                         tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTg          <  5:226524/146‑1 (MQ=255)                                                                                                                                          cATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGc       >  6:180253/1‑149 (MQ=255)                                                                                                                                          cATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGc       >  8:48463/1‑149 (MQ=255)                                                                                                                                             cGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAg    >  6:16781/1‑149 (MQ=255)                                                                                                                                               gACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTg  >  4:152040/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                               TCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGT‑GCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAG‑CGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCGCGCGAAGAAACCCCTGT‑CTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTG  >  NC_002947/2211762‑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: ‑