Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,382,535	G→T	8.3%	L256I (CTC→ATC)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,382,535	0	G	T	8.3%	146.2 / 2.1	47	L256I (CTC→ATC)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base T (2/2);  total (23/24)
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.94e-01

AGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAC  >  NC_002947/1382396‑1382683                                                                                                                                            |                                                                                                                                                     aGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                               >  4:308766/1‑147 (MQ=255) aGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                               >  7:71634/1‑147 (MQ=255)   gAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGa                                                                                                                                           <  5:211010/149‑1 (MQ=255)     aCTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATc                                                                                                                                         <  6:76689/149‑1 (MQ=255)      cTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGATGACCTCGAAGAtga                                                                                                                                        <  1:191963/149‑3 (MQ=255)         ctgctgCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGc                                                                                                                                      >  5:111287/1‑148 (MQ=255)            ctgcAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTc                                                                                                                                    <  8:71634/147‑1 (MQ=255)                aGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGAACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGc                                                                                                                              >  1:161452/1‑149 (MQ=255)                                       cTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTc                                                                                                                           <  4:77124/129‑1 (MQ=255)                                          gATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTgg                                                                                                    >  7:93081/1‑149 (MQ=255)                                             cAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGgcgc                                                                                                 >  2:238242/1‑149 (MQ=255)                                                gCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGa                                                                                              >  5:284540/1‑149 (MQ=255)                                                 ctcgGTTCCATTTCACTGGCTTCCTCTTCACCGCCGAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGTCATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGAc                                                                                             <  7:81699/149‑1 (MQ=255)                                                        ccATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGa                                                                                              <  3:45128/141‑1 (MQ=255)                                                        ccATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGa                                                                                              >  4:45128/1‑141 (MQ=255)                                                         cATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTc                                                                                                                                    >  4:51686/1‑102 (MQ=255)                                                         cATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGc                                                                                     <  6:111287/149‑1 (MQ=255)                                                           tttCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                                   <  1:238242/149‑1 (MQ=255)                                                                      tCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCg                                                                        >  4:273431/1‑149 (MQ=255)                                                                       ccTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTc                                                                         <  8:93081/147‑1 (MQ=255)                                                                        ctctTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTa                                                                      <  5:171219/149‑1 (MQ=255)                                                                          ctTCAAGGCCCAGTTCGGACAGGCGCTTGACACGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGTTGGAGCTCACCAACAGGCTGGCATCGCCTCGTaac                                                                    >  4:100295/1‑147 (MQ=255)                                                                           ttCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGAt                                                                                                                                          <  4:40001/78‑1 (MQ=255)                                                                             cAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTAcc                                                                    <  6:284540/146‑1 (MQ=255)                                                                              aaGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCa                                                                 <  3:308766/148‑1 (MQ=255)                                                                                       ttCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGt                                                       <  3:100295/149‑1 (MQ=255)                                                                                         cGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTc                                                      <  8:4316/148‑1 (MQ=255)                                                                                           gACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAg                                                   >  1:322578/1‑149 (MQ=255)                                                                                                gcgcTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAAc                                                                                         <  1:86750/106‑1 (MQ=255)                                                                                                gcgcTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAAc                                                                                         >  2:86750/1‑106 (MQ=255)                                                                                                       aCCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtg                                       >  2:11421/1‑149 (MQ=255)                                                                                                        cccGGTGCACCCACTGGCCGACATCATCA‑‑‑TCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcagg                                   <  4:241317/149‑1 (MQ=255)                                                                                                        cccGGTGCACCCACTGGCCGACATCATCA‑‑‑TCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcagg                                   >  3:241317/1‑149 (MQ=255)                                                                                                                aCCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcagg                                   <  5:38689/144‑1 (MQ=255)                                                                                                                aCCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcagg                                   >  6:38689/1‑144 (MQ=255)                                                                                                                  ccACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCTAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTACCGa                            >  5:77524/1‑149 (MQ=255)                                                                                                                     cTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTTCTTGCAGGTGCAGGGTGCCGACGt                         <  5:146422/149‑1 (MQ=255)                                                                                                                     cTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGt                         <  2:322578/149‑1 (MQ=255)                                                                                                                       ggCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTgg                       >  1:265677/1‑149 (MQ=255)                                                                                                                         ccGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                                   >  4:204840/1‑87 (MQ=255)                                                                                                                           gACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGa                   <  1:11421/149‑1 (MQ=255)                                                                                                                                  tgatcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgc                                            >  6:233869/3‑117 (MQ=255)                                                                                                                                  tgatcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgc                                            <  5:233869/115‑1 (MQ=255)                                                                                                                                  tgatcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACAAACAGGCTGGCATAGCCTCGTACCGCATTGTCCTCGTAAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTg            >  5:340028/3‑149 (MQ=255)                                                                                                                                     tcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTgg                       >  2:281938/1‑135 (MQ=255)                                                                                                                                     tcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTgg                       <  1:281938/135‑1 (MQ=255)                                                                                                                                     tcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGGTc         >  3:195541/1‑149 (MQ=255)                                                                                                                                            ggTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAc  >  4:131881/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                     AGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAC  >  NC_002947/1382396‑1382683

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