Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,048,529	A→T	11.2%	Q195L (CAG→CTG)	PP_2661 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,048,529	0	A	T	11.2%	115.0 / 4.1	43	Q195L (CAG→CTG)	PP_2661	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/21);  new base T (3/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGC  >  NC_002947/3048381‑3048672                                                                                                                                                     |                                                                                                                                                ggTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCa                                                                                                                                                 <  6:307641/149‑1 (MQ=255)       gACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACt                                                                                                                                           <  2:136226/149‑1 (MQ=255)                  ccGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTg                                                                                                                                  >  4:6592/1‑147 (MQ=255)                   cgcgTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGaa                                                                                                                                <  2:272960/148‑1 (MQ=255)                    gcgTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGc                                                                                                                              >  6:105548/1‑149 (MQ=255)                        ggtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAgg                                                                                                                          <  1:31138/149‑1 (MQ=255)                         gtggtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCAGGAAGCCAgt                                                                                                                          <  7:36005/148‑2 (MQ=255)                            gtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTg                                                                                                                                          >  1:54378/1‑129 (MQ=255)                            gtggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTg                                                                                                                                          <  2:54378/129‑1 (MQ=255)                             tggtgTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGAtt                                                                                                                     >  2:108968/1‑149 (MQ=255)                                      cGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCt                                                                                                                >  5:98535/1‑145 (MQ=255)                                       gAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGa                                                                                                            >  3:60553/1‑148 (MQ=255)                                            aTGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGaa                                                                                                                                >  5:231019/1‑123 (MQ=255)                                            aTGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGaa                                                                                                                                <  6:231019/123‑1 (MQ=255)                                                 cACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTc                                                                                                                    >  2:143071/1‑130 (MQ=255)                                                 cACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTc                                                                                                                    <  1:143071/130‑1 (MQ=255)                                                      ggCGCAGAAGGTGAAGGACGCGGTTGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGgc                                                                                             <  3:6592/148‑1 (MQ=255)                                                           aGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCttt                                                                                       <  6:21803/149‑1 (MQ=255)                                                           aGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCttt                                                                                       <  8:242170/149‑1 (MQ=255)                                                                   aaGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCt                                                                               <  4:60553/149‑1 (MQ=255)                                                                         gcgGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCTGGCACTGGCTTCCTGAAGCCAGGAGATTCCGc                                                                                                                 <  5:308578/109‑1 (MQ=255)                                                                                 ggATGTGCAGGCACGGCTGGTCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGccc                                                                 >  4:212318/1‑149 (MQ=255)                                                                                        cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGa                                                                                                                                 <  1:150116/78‑1 (MQ=255)                                                                                        cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGa                                                                                                                                 >  2:150116/1‑78 (MQ=255)                                                                                           gCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCGGAAGCCAGGAGATTCCGCTCAGCTAATGTAATTGGGGCTGCTTTGCAGCCATTCGCCGGCGAGCCCGCGAAGGGcc                                                       >  8:211103/1‑149 (MQ=255)                                                                                                 ctggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGTCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGgctgct                                                                                         <  8:151510/109‑1 (MQ=255)                                                                                                    gcatggcAGCAACCGCCCCGTTCTGAACTGCCCGCTGGCCAGCGTGCGCAGGCACTGGCGTCCTGACGCCAGGCGATGCCGCTCAGATAATGTCATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggc                                                <  7:211103/144‑1 (MQ=255)                                                                                                      ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcc                                            >  3:191714/1‑149 (MQ=255)                                                                                                      ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCAAAGCAGccc                                                                                >  7:105662/1‑113 (MQ=255)                                                                                                      ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCAAAGCAGccc                                                                                <  8:105662/113‑1 (MQ=255)                                                                                                          cagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCgg                                        <  6:98535/149‑1 (MQ=255)                                                                                                                   cACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCaa                               >  4:19167/1‑149 (MQ=255)                                                                                                                     ccGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCAGGTTCATTCAATc                             >  3:15959/1‑149 (MQ=255)                                                                                                                        ttATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                          <  3:212318/149‑1 (MQ=255)                                                                                                                          aTGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAt                         <  3:19167/148‑1 (MQ=255)                                                                                                                           tGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                        <  3:76179/148‑1 (MQ=255)                                                                                                                                  agccagtGCCTGCGAACGCTGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCg                >  8:217173/7‑149 (MQ=255)                                                                                                                                   ccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                        <  8:277711/140‑1 (MQ=255)                                                                                                                                   ccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                        >  7:277711/1‑140 (MQ=255)                                                                                                                                     cagtGCCTGCGAACGCTGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt             >  6:32446/4‑149 (MQ=255)                                                                                                                                      agtGCCTGCGAACGCTGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTg            >  8:273061/3‑149 (MQ=255)                                                                                                                                           cAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCACCGGTTCATTCAATCGTATTGATACCCGCCTGCCCCt       >  7:129150/1‑149 (MQ=255)                                                                                                                                              cgaacGCTGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTc      <  7:273061/143‑1 (MQ=255)                                                                                                                                              cGTTCGCAGGCACTGGCTTGCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGCCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcct    <  1:234650/149‑1 (MQ=255)                                                                                                                                                ttggtaGGCACTGGCTTCTTGAAGCTAGGAGATTCCGTTCAGATAATGTAATTGGGGCTGCTTTGCAGCCGTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcctgc  <  4:15959/144‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GGTGTCGACCACGGTTTCCGCGTGGTGGTGGTGTGCCCGAAGGATGAGCACTGGGCGCAGAAGGTGAAGGACGCGGTGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGC  >  NC_002947/3048381‑3048672

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