Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,198	Δ1 bp	8.9%	coding (772/1320 nt)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,054,198	0	A	.	8.9%	169.0 / 10.9	56	coding (772/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base . (3/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGG  >  NC_002947/4054056‑4054339                                                                                                                                                 |                                                                                                                                              cgACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTgg                                                                                                                                         <  2:88795/149‑1 (MQ=255)  gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGt                                                                                                                                        <  6:34841/149‑1 (MQ=255)       ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGc                                                                                                                                    <  1:255613/148‑1 (MQ=255)        tgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGAGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCt                                                                                                                                  <  6:34068/149‑1 (MQ=255)          ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGa                                                                                                                                <  6:110933/149‑1 (MQ=255)          ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGa                                                                                                                                <  8:245648/149‑1 (MQ=255)             gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATG‑CCCTGGTGGGcgaggg                                                                                                                               <  7:204177/147‑6 (MQ=255)                 cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTc                                                                                                                         <  3:128092/149‑1 (MQ=255)                  tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCTACCAGGGTCATG‑CCCTGGTggg                                                                                                                                     <  1:277775/136‑1 (MQ=255)                  tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCTACCAGGGTCATG‑CCCTGGTggg                                                                                                                                     >  2:277775/1‑136 (MQ=255)                     gCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAcc                                                                                                                              <  5:292271/140‑1 (MQ=255)                     gCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAcc                                                                                                                              >  6:292271/1‑140 (MQ=255)                         gCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaa                                                                                                                  >  6:193769/1‑148 (MQ=255)                                     ggACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgcgc                                                                                                     >  5:180472/1‑149 (MQ=255)                                       aCGCGGAAATTGCCAAGCCCGGTCGCGAGCATGGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCAGCCTGACCCAGTCGTTCGAAAACGGc                                                                                                            >  1:43068/1‑141 (MQ=255)                                       aCGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCgcgcgc                                                                                                   <  4:303342/149‑1 (MQ=255)                                       aCGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCgcgcgc                                                                                                   <  3:349826/149‑1 (MQ=255)                                       aCGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCAGCCTGACCCAGTCGTTCGAAAACGGc                                                                                                            <  2:43068/140‑1 (MQ=255)                                         gcgGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                 >  1:139590/1‑149 (MQ=255)                                         gcgGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                 >  6:316468/1‑149 (MQ=255)                                          cgGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                                <  6:41521/149‑1 (MQ=255)                                              aaTTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGTCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAgg                                                                                            <  5:193769/149‑1 (MQ=255)                                              aaTTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAcc                                                                                                                              <  5:191384/115‑1 (MQ=255)                                              aaTTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAcc                                                                                                                              >  6:191384/1‑115 (MQ=255)                                                   ccAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAc                                                                                       <  5:256665/149‑1 (MQ=255)                                                    cAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                      >  7:43861/1‑149 (MQ=255)                                                    cAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                      >  8:18204/1‑149 (MQ=255)                                                        cccGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGt                                                                                  <  1:229599/149‑1 (MQ=255)                                                         ccGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTc                                                                                 <  2:139590/149‑1 (MQ=255)                                                             tCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCgaga                                                                             <  1:206022/149‑1 (MQ=255)                                                                 gAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                         >  3:89326/1‑149 (MQ=255)                                                                        gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                   >  1:153241/1‑149 (MQ=255)                                                                         tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcgc                                                                  >  4:26284/1‑149 (MQ=255)                                                                           aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcgg                                                                >  5:180630/1‑149 (MQ=255)                                                                             cATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                               >  5:359233/1‑148 (MQ=255)                                                                              aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTGCCGTGCCCTCGCCCACCAGGG‑CATGACCCTCGACGGCCTGACCCAGTCGTACGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcga                                                             >  8:253783/1‑149 (MQ=255)                                                                                                    gtggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                        <  1:75806/148‑1 (MQ=255)                                                                                                     tggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTgg                                      <  6:180630/149‑1 (MQ=255)                                                                                                       gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGAt                                    <  5:222476/149‑1 (MQ=255)                                                                                                              cggggACTTCCGTGCCCGCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACTGCGTCGCGCGCTTCCATGACAACCTGGTCGAGAACATCAATCgcgg                                                                <  4:37741/111‑1 (MQ=255)                                                                                                              cGTGGACTTACGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGCCCCAGTCGTGCGAAAACGGAGTCGCGCGCTTCCAGGACAACCGGGTCGAGAACATCAATCgcgg                                                                >  3:37741/1‑114 (MQ=255)                                                                                                               gTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                              >  3:235884/1‑67 (MQ=255)                                                                                                               gTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                              <  4:235884/67‑1 (MQ=255)                                                                                                                tGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGAt                           <  3:26284/149‑1 (MQ=255)                                                                                                                   aCTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGcc                        <  7:18204/149‑1 (MQ=255)                                                                                                                    cTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                       <  6:359233/149‑1 (MQ=255)                                                                                                                         gTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                               <  3:28968/120‑1 (MQ=255)                                                                                                                         gTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGCTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc                  <  3:230528/149‑1 (MQ=255)                                                                                                                         gTGCCCTCGCCCACCAGGG‑CATGACACTGGTCGGCCTGACCCAGTCGATCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                               >  4:28968/1‑120 (MQ=255)                                                                                                                             ccTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAAc              >  3:241145/1‑149 (MQ=255)                                                                                                                             ccTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAAc              >  7:284044/1‑149 (MQ=255)                                                                                                                                   cgaccaGGGTCATG‑CCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACg        >  2:128337/3‑149 (MQ=255)                                                                                                                                   cgaccaGGGTCATG‑CCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGACGCGCGCATCCAGGACAACCTGGTCGAGAACAACAATCGAGGCGACGACCACTACCTGGACCTGCAGCATGCGGCCGATGCCTACATCGAACGCaa          >  4:112288/3‑147 (MQ=255)                                                                                                                                   ccaccaGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                          >  6:156754/1‑99 (MQ=255)                                                                                                                                   ccaccaGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                          <  5:156754/99‑1 (MQ=255)                                                                                                                                         ggg‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTgg  <  4:241145/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              CGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCAT‑GTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGG  >  NC_002947/4054056‑4054339

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