Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,118,338	Δ1 bp	7.3%	coding (1028/2256 nt)	PP_3622 →	isoquinoline 1‑oxidoreductase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,118,338	0	A	.	7.3%	167.5 / 7.0	55	coding (1028/2256 nt)	PP_3622	isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base . (2/2);  total (27/28)
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.02e-01

GCGTGCAGTGGAGGCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGCCCAA  >  NC_002947/4118195‑4118485                                                                                                                                                |                                                                                                                                                       gCGTGCAGTGGAGGCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGGCGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAggg                                                                                                                                                    >  6:19035/1‑148 (MQ=255)         tGGAGGCGATTCAGGTGGACTGTAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGcc                                                                                                                                           >  7:57378/1‑149 (MQ=255)         tGGAGGCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGcc                                                                                                                                           >  1:358866/1‑149 (MQ=255)           gAGGCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAACAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCgg                                                                                                                                         >  1:173845/1‑149 (MQ=255)            aGGCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGgc                                                                                                                                        <  5:347516/149‑1 (MQ=255)              gCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGgcgc                                                                                                                                      <  8:301120/149‑1 (MQ=255)                 aTTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGTCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑tt                                                                                                                                 <  3:323500/149‑1 (MQ=255)                  ttCAGGTCGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTg                                                                                                                                <  8:324320/149‑1 (MQ=255)                     aGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGc                                                                                                                              >  7:28928/1‑148 (MQ=255)                     aGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGc                                                                                                                              <  6:199420/148‑1 (MQ=255)                     aGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGc                                                                                                                              <  3:361823/148‑1 (MQ=255)                             tGGAAAGAACCGGCTGCCGACAGCAAAGTGCCACCGATGCCGGACGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCACGAGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCAAACGCCa                                                                                                                      >  5:62037/1‑148 (MQ=255)                               gAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATAtt                                                                                                                                              <  8:265033/124‑1 (MQ=255)                               gAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATAtt                                                                                                                                              >  7:265033/1‑124 (MQ=255)                                    aaCCGGCGGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAg                                                                                                             <  8:299034/149‑1 (MQ=255)                                             ccGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAgg                                                                                                      <  8:218071/147‑1 (MQ=255)                                             ccGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGc                                                                                                     <  7:27677/148‑1 (MQ=255)                                             ccGACAGCAAAGTGAGACCGATGCGGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGGGACGGTGCAAGGGCGGTCCCGTGATGAAGAGTAGGAAGGGGATTTTGCAGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAg                                                                                                       <  5:329458/146‑1 (MQ=255)                                                 cAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGc                                                                                                 >  4:223283/1‑148 (MQ=255)                                                     aaaGTGCGACCGATGCCGGCCGATTTTTCCAGAGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACca                                                                                            >  4:45677/1‑149 (MQ=255)                                                        gTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACaa                                                                                         >  1:84188/1‑149 (MQ=255)                                                         tGCGACCGATGCCGTCCGTTTTTTCCAGCGAGGCCTGGCGCGAACTCCTCGCGACGGTGCAAGGGCCGTCCGGTTATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACgacatc                                                                                        <  6:62037/149‑3 (MQ=255)                                                         tGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAAc                                                                                        <  8:28928/149‑1 (MQ=255)                                                               cGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAg                                                                                                       <  8:134701/128‑1 (MQ=255)                                                               cGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAg                                                                                                       >  7:134701/1‑128 (MQ=255)                                                                     cggccgATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACCACCAGTACCTGacc                                                                            >  7:245743/1‑147 (MQ=255)                                                                              tttCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCa                                                                    <  2:84188/148‑1 (MQ=255)                                                                                   aGCGACGCCTGGCGCGAACAACTCCCGAATGTCCAACGGCCGACCCCTGATGAAGAGAAAGAAGGGGATATTGCCGTCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGACGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTg                                                                >  5:180351/1‑147 (MQ=255)                                                                                    gCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAgg                                                                                                            >  5:335559/1‑102 (MQ=255)                                                                                    gCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAgg                                                                                                            <  6:335559/102‑1 (MQ=255)                                                                                       aCGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCa                                                          <  2:358866/149‑1 (MQ=255)                                                                                                  gAACACCTCGCGACGCTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTACAGGCCAGCTACCACAACCAGTACCTGAACCATGCACAGCTGGAGCCACCGTCCGCCCt                                               >  8:289016/1‑149 (MQ=255)                                                                                                       ccTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCCGGTTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGc                                            >  1:293733/1‑149 (MQ=255)                                                                                                       ccTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGccc                                          <  4:333412/149‑1 (MQ=255)                                                                                                         tCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGACCAGCTGGAGCCACCGTCCGCCCTCGCCCGc                                        >  5:192917/1‑149 (MQ=255)                                                                                                                  tGCAAGGGCCGTCCCGTGATTAAGAGAAAGAAGGGGATATTGCCGGCGCC‑‑TTGGCCACCGCAAAGAC‑CAAGGTCGAGGCCAGCTACCACAACACGTACCTGACCCAGGCGCAGCTGGAGCCACCGTCCGCACTCGCCCGCTTCAAACCg                               >  8:341943/1‑149 (MQ=255)                                                                                                                  tGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACcac                                                                                           >  1:321165/1‑89 (MQ=255)                                                                                                                  tGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACcac                                                                                           <  2:321165/89‑1 (MQ=255)                                                                                                                      aGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCTCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCGAGCTACCACAACCAGTACCTGAACCATGCCCTGCTGGAGCCAGCGTCCGCCCTCGCCGGCTTCAACCCGGAc                            <  3:11060/148‑1 (MQ=255)                                                                                                                      aGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGACGGCGCC‑‑TTGGCCAACGCCAAGAC‑AAAGGTCGAGGCCAGCTACCACAACAAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTAGCCCGCTTAAACCCGGACg                           >  7:5515/1‑149 (MQ=255)                                                                                                                       gggCCGTCCCGTGATGAAGAGAACGAAGGTGATATTGCCGGCGCC‑‑TTGGCCAACGCCCAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTAACTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCGTCACCCCAGACga                          >  4:11060/1‑148 (MQ=255)                                                                                                                       gggCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACgg                          <  2:173845/149‑1 (MQ=255)                                                                                                                                cGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGa                 >  3:363232/1‑149 (MQ=255)                                                                                                                                       aaGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGACGCAAGGTCGAGGCCAGCTACCACAACCAGAACCTGAACCATGCCCCGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCa                                    >  1:82241/1‑124 (MQ=255)                                                                                                                                       aaGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCa                                    <  2:82241/123‑1 (MQ=255)                                                                                                                                          agaACG‑AGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGc      >  3:203655/1‑149 (MQ=255)                                                                                                                                          agaACG‑AGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGc      <  4:203655/149‑1 (MQ=255)                                                                                                                                          agaACG‑AGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACAACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGc      >  5:170252/1‑149 (MQ=255)                                                                                                                                           gaACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGc      <  4:172628/149‑1 (MQ=255)                                                                                                                                           gaACG‑AGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGc      <  6:170252/148‑1 (MQ=255)                                                                                                                                            aaCGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCt            >  5:196611/1‑142 (MQ=255)                                                                                                                                            aaCGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCt            <  6:196611/142‑1 (MQ=255)                                                                                                                                             aCGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGccc    <  4:363232/149‑1 (MQ=255)                                                                                                                                              cGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCAGAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGccc    >  8:345944/1‑148 (MQ=255)                                                                                                                                               gAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGCCCaa  <  5:19035/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                       GCGTGCAGTGGAGGCGATTCAGGTGGACTGGAAAGAACCGGCTGCCGACAGCAAAGTGCGACCGATGCCGGCCGATTTTTCCAGCGACGCCTGGCGCGAACACCTCGCGACGGTGCAAGGGCCGTCCCGTGATGAAGAGAACGAAGGGGATATTGCCGGCGCC‑‑TTGGCCAACGCCAAGAC‑CAAGGTCGAGGCCAGCTACCACAACCAGTACCTGAACCATGCCCAGCTGGAGCCACCGTCCGCCCTCGCCCGCTTCAACCCGGACGGCAGCCTCGACGTCTGGCTGCCCAA  >  NC_002947/4118195‑4118485

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