Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,053,476	A→T	10.3%	Q17L (CAG→CTG)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,053,476	0	A	T	10.3%	119.0 / 5.9	49	Q17L (CAG→CTG)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (19/23);  new base T (3/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGG  >  NC_002947/4053338‑4053623                                                                                                                                           |                                                                                                                                                    gCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCgc                                                                                                                                            <  2:293663/148‑2 (MQ=255) gCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTg                                                                                                                                           >  2:54907/1‑149 (MQ=255) gCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTg                                                                                                                                           <  3:258431/149‑1 (MQ=255)  ccATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTgg                                                                                                                                          >  5:311958/1‑149 (MQ=255)  ccATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTgg                                                                                                                                          <  3:214253/149‑1 (MQ=255)  ccATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCCGCCAGGCCGGCGTgg                                                                                                                                          >  4:230624/1‑149 (MQ=255)             gCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGaa                                                                                                                               <  5:222251/149‑1 (MQ=255)              cATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAac                                                                                                                              <  2:125202/149‑1 (MQ=255)                  ccTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCt                                                                                                                          >  5:241485/1‑149 (MQ=255)                  ccTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCt                                                                                                                          <  5:247340/149‑1 (MQ=255)                  ccTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCt                                                                                                                          <  5:309877/149‑1 (MQ=255)                               cAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTg                                                                                                               <  6:241485/147‑1 (MQ=255)                               cAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGc                                                                                                             <  5:168134/149‑1 (MQ=255)                                    aaCGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGcc                                                                                                        >  6:156831/1‑149 (MQ=255)                                        gcggTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAc                                                                                                    >  3:102300/1‑149 (MQ=255)                                         cggTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAcc                                                                                                   <  5:156831/149‑1 (MQ=255)                                         cggTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAcc                                                                                                   >  1:235837/1‑149 (MQ=255)                                         cggTCCTCCAAACCGAACTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGATCAAGCTTGGCGTGCCGCAcc                                                                                                   <  6:4095/149‑1 (MQ=255)                                            tcctccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTgg                                                                                                <  1:54907/149‑1 (MQ=255)                                                ccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCt                                                                                                                 >  6:173496/1‑128 (MQ=255)                                                ccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCt                                                                                                                 <  5:173496/128‑1 (MQ=255)                                                ccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCACGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACt                                                                                            >  6:3657/1‑149 (MQ=255)                                                    aCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAg                                                                                        <  6:311958/149‑1 (MQ=255)                                                     ccGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAgc                                                                                       <  2:235837/149‑1 (MQ=255)                                                        aCCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCa                                                                                    <  4:102300/149‑1 (MQ=255)                                                           tGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAAcc                                                                                 >  2:129154/1‑149 (MQ=255)                                                             aCAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCg                                                                                >  8:217945/1‑148 (MQ=255)                                                                 aaaCCGCCTGAACAAGGGCTCGGACATCACACTGAACAACCGCGAAATAGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAACATGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCaact                                                                           >  2:219019/1‑146 (MQ=255)                                                                            aCAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCCGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTgg                                                                >  5:190523/1‑149 (MQ=255)                                                                                gggCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTa                                                            >  6:128308/1‑149 (MQ=255)                                                                                      ggACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCgg                                                        <  2:244258/147‑1 (MQ=255)                                                                                            gACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTgg                                                 >  4:297699/1‑148 (MQ=255)                                                                                             acacTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTggg                                                <  1:57552/148‑1 (MQ=255)                                                                                             acacTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGa                                               <  3:297699/149‑1 (MQ=255)                                                                                              cacTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCCGCGTGGCCATGAGCGAACACCAGAGCAAGCTTTGCGTGCCGCACCTGGTAATGGAGCGCAACCGCATAGCGCAGGCCTGGCGTACCCGCCGCTGGGa                                               >  4:135585/1‑148 (MQ=255)                                                                                                cTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCCGCCCGGCCGGCCTGGCCCAGAGCGCACACCATAGCAAACTTGGCGGGGCGCCCCGGGGACTGGGGCGCCACCGCAAAGCGCACGCCAGGCGTAACCGCCCCAggggc                                              >  8:109196/1‑145 (MQ=255)                                                                                                         cTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCCGGGCGGCGTGGGCAAGAGGGCAAAACGTAGCAAACTTTGGGTGCCGCCCCTGGTAATGGGGCGCAACCGCAAAGCGCAGGCCTGGCGTTACCGCCGCTGGGGCTCGCGGGtggg                                   >  6:38652/1‑146 (MQ=255)                                                                                                                       ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAg                                                                      <  1:212365/100‑1 (MQ=255)                                                                                                                       ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTg                                                  <  1:133504/120‑1 (MQ=255)                                                                                                                       ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCt                     <  5:38652/149‑1 (MQ=255)                                                                                                                       ccTCGTCGTCGGCGCCGGCCACGCCCGCGTGGGCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCCCCGGGTACTGGAGCGCAACCCCATCGCGCAg                                                                      >  2:212365/1‑100 (MQ=255)                                                                                                                                  acgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGc          >  4:126753/2‑149 (MQ=255)                                                                                                                                  acgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGc          <  7:229607/148‑1 (MQ=255)                                                                                                                                  acgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGc          >  8:229607/2‑149 (MQ=255)                                                                                                                                  acgcCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCAAGACCGc          >  8:237683/2‑149 (MQ=255)                                                                                                                                   cgcCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCt         >  6:216179/1‑149 (MQ=255)                                                                                                                                        gCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGTCGTACCGGCCGCTGGGACTCGCTGTTCGCCAAcg                              <  4:11763/123‑1 (MQ=255)                                                                                                                                        gCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAAcg                              >  3:11763/1‑123 (MQ=255)                                                                                                                                          cAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCggg  <  7:109196/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    GCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGGGACTCGCTGGTCGCCAACGGCCCGGCCTGGCACGACCGCTTCCCGGG  >  NC_002947/4053338‑4053623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑