Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,156,328	A→C	9.6%	S392A (TCT→GCT)	PP_2772 ←	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,156,328	0	A	C	9.6%	100.0 / 4.2	42	S392A (TCT→GCT)	PP_2772	monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (17/21);  new base C (2/2);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CTTCGCTGGGGGCGACCTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGCGAGGCGTGGGTCGCGCGTGCGCCAGTGA  >  NC_002947/3156179‑3156471                                                                                                                                                      |                                                                                                                                                cTTCGCTGGGGGCGACCTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGa                                                                                                                                                 >  1:495766/1‑150 (MQ=255)         ggggCGACCTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGgggtgg                                                                                                                                         <  1:1100271/150‑1 (MQ=255)             cGACCTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGCTGGGGTGGGTcg                                                                                                                                     >  2:871540/1‑150 (MQ=255)             cGACCTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGATCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGCTGGGGTGGGTcg                                                                                                                                     >  2:871622/1‑150 (MQ=255)                ccTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCgt                                                                                                                                  <  2:747384/150‑1 (MQ=255)                  tGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCgtgt                                                                                                                                >  2:935896/1‑150 (MQ=255)                   gTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGtgtg                                                                                                                               <  1:268194/150‑1 (MQ=255)                    ttGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTgg                                                                                                                              >  2:1118304/1‑150 (MQ=255)                     tGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGCTGGGGTGGGTCGCGTGTGGc                                                                                                                             <  1:871540/150‑1 (MQ=255)                     tGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTgga                                                                                                                             <  2:857732/150‑2 (MQ=255)                     tGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGATCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGCTGGGGTGGGTCGCGTGTGGc                                                                                                                             <  1:871622/150‑1 (MQ=255)                               cGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCa                                                                                                                   <  1:447816/150‑1 (MQ=255)                               cGAGCAGCGCCTGCGATCTAGCGTCTCCCAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCa                                                                                                                   <  1:280741/150‑1 (MQ=255)                                          tGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCt                                                                                                        >  1:529851/1‑150 (MQ=255)                                                 tAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTgg                                                                                                                              <  1:60533/121‑1 (MQ=255)                                                 tAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTgg                                                                                                                              >  2:60533/1‑121 (MQ=255)                                                 tAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAg                                                                                                 >  2:1047484/1‑150 (MQ=255)                                                  aGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTTCCAAATTGCCCACTCACCCAGCTCGTGGAGt                                                                                                <  1:577962/150‑1 (MQ=255)                                                   gCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGtt                                                                                               >  2:223338/1‑150 (MQ=255)                                                      tgttCAAAGACCTGTGCGGTTGTCTGGATCGATGTCGTCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGTGATTCCGATGGGTAGAAAGATGGGGTGGGTCGCGGGTGTCGAAGTGTCCACTCACCTAGCTCTTGGCGTTTGt                                                                                            >  2:520585/3‑150 (MQ=255)                                                                 cTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCggg                                                                                 >  2:310380/1‑150 (MQ=255)                                                                   gtgCGGTTGCCCGGGTCGATGGCGCCAGTGACGGCGAAGGGGGGAAGGTCAGCAATTGCGTGGGTATCCCTATGGGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTc                                                                               <  1:976660/150‑1 (MQ=255)                                                                             ccGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGttt                                                                                              <  2:567809/125‑1 (MQ=255)                                                                             ccGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGttt                                                                                              >  1:567809/1‑125 (MQ=255)                                                                                gATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTTGGTCGCGGGTGTCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTTTAGTCTAGCGTGTCGTGCCGGCTCTgg                                                                  >  2:790096/1‑150 (MQ=255)                                                                                gATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTgg                                                                  <  2:452887/150‑1 (MQ=255)                                                                                       gTCGACAGTGACGGCGAAGCGGGATGTGTTATCAATTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTACCCACTCACCCGGCTCGTGCATTTTGTAGTCCAGCGTGTCGTGCTGGCTCTGGGTACTCa                                                           <  2:529851/150‑1 (MQ=255)                                                                                               tGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTgcgc                                                   <  1:750726/150‑1 (MQ=255)                                                                                                 aCGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCa                                                 <  1:696958/150‑1 (MQ=255)                                                                                                  cGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAg                                                <  1:223338/150‑1 (MQ=255)                                                                                                   ggCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGACGGCTCTGGGTAAGCAGGTTGCGCCAgg                                               >  2:906271/1‑149 (MQ=255)                                                                                                    gCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTa                                              <  2:495766/150‑1 (MQ=255)                                                                                                     cGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAg                                             <  1:1047484/150‑1 (MQ=255)                                                                                                       aaGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCg                                           >  2:167639/1‑150 (MQ=255)                                                                                                             ggAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAg                                     <  1:1118304/150‑1 (MQ=255)                                                                                                                    cAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGc                              <  2:957628/150‑1 (MQ=255)                                                                                                                    cAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGc                              >  2:186493/1‑150 (MQ=255)                                                                                                                      gCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCATGTTGCGCCAGTAGCGATCGCGGCGCCGCGa                            >  2:159982/1‑150 (MQ=255)                                                                                                                              aTGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGGGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGAGAGGCGTggg                    >  2:97068/1‑150 (MQ=255)                                                                                                                                    aTTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGCGAGGCGTGGGTcgcgc              >  1:178243/1‑150 (MQ=255)                                                                                                                                            gAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATc                                        >  1:196545/1‑116 (MQ=255)                                                                                                                                            gAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATc                                        <  2:196545/116‑1 (MQ=255)                                                                                                                                                aGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGCGAGGCGTGGGTCGCGCGTGCGCCAGTGa  <  1:310380/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                CTTCGCTGGGGGCGACCTGTTGCAGGTGGGCGAGCAGCGCCTGCGATCTAGCGTCTTCAAAGACCTGTGCGGTTGCCCGGATCGATGTCGCCAGTGACGGCGAAGGGGGGAGGGTCAGCAACTGCATGGCGATTCCTATGAGTAGAAAGATGGGGTGGGTCGCGTGTGGCTCAGTGCCCACTCACCCAGCTCGTGGAGTTTGTAGTCCAGCGGGTCGTGCAGGCTCTGGGTACGCAGGTTGCGCCAGTAGCGATCGAGGCGCAGCGAGGCGTGGGTCGCGCGTGCGCCAGTGA  >  NC_002947/3156179‑3156471

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