Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,110,826	A→G	26.7%	L868L (CTA→CTG)	PP_3616 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,110,826	0	A	G	26.7%	16.1 / 4.9	15	L868L (CTA→CTG)	PP_3616	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (4/7);  new base G (2/2);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACGCCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCTATCACACATAACGAAGTAGCCACCCCCCAACCAACCGTAGA  >  NC_002947/4110683‑4110957                                                                                                                                                |                                                                                                                                    cGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACGCCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATccgcc                                                                                                                               <  1:436339/150‑1 (MQ=255)                        cAATGATGGAGAGGCGCACGCCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCCCTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCAt                                                                                                        <  2:194434/150‑1 (MQ=255)                         aaTGATGGAGAGGCGCACGCCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTGTCCGCCGAGCAGGCCCTATCCCTCGGCATa                                                                                                       >  1:196162/1‑150 (MQ=255)                                                   cAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCCCAATGAACATTACGGGTTCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCa                                                                             >  1:30614/1‑150 (MQ=255)                                                           aCTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAAcc                                                                     <  1:69341/150‑1 (MQ=255)                                                                            cctgcTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTGTCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTgg                                                                >  2:428153/1‑138 (MQ=255)                                                                            cctgcTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTGTCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTgg                                                                <  1:428153/138‑1 (MQ=255)                                                                                     aCCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCt                                           >  2:139741/1‑150 (MQ=255)                                                                                     aCCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACAACCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCt                                           <  1:380018/150‑1 (MQ=255)                                                                                       cTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGAt                                                                                    >  1:443812/1‑107 (MQ=255)                                                                                       cTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGAt                                                                                    <  2:443812/107‑1 (MQ=255)                                                                                        tCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACGCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCTATc                                        >  2:330004/1‑150 (MQ=255)                                                                                            tGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCTATcacac                                    <  1:139741/150‑1 (MQ=255)                                                                                                              cTTGGCACAATGAACATTACGGGCCCAAGCTACTGTCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCTATCACACATAACGAAGTAGCCAccc                  <  2:196162/150‑1 (MQ=255)                                                                                                                              ttACGAGCCCAAGATACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCTATCACACATAACGAAGTAGCCACCCCCCAACCAACCGTaga  <  1:195820/150‑1 (MQ=255)                                                                                                                                           ctactaTCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCg                                                                                      >  2:420875/1‑53 (MQ=255)                                                                                                                                           ctactaTCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCg                                                                                      <  1:420875/53‑1 (MQ=255)                                                                                                                                                |                                                                                                                                    CGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACGCCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAATTAAGGCTTGGCACAATGAACATTACGGGCCCAAGCTACTATCCGCCGAGCAGGCCCTATCCCTCGGCATAACTAAAGCCGACACCCGATGGCATCAGAATAACCACTGGTCGCCGCCAAATTGGACAGCTATCACACATAACGAAGTAGCCACCCCCCAACCAACCGTAGA  >  NC_002947/4110683‑4110957

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