Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,697	A→G	12.7%	F63L (TTC→CTC)	PP_1879 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,104,697	0	A	G	12.7%	133.4 / 10.8	55	F63L (TTC→CTC)	PP_1879	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (5/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

AATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGA  >  NC_002947/2104551‑2104841                                                                                                                                                   |                                                                                                                                                 aaTAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACt                                                                                                                                                <  5:269675/149‑1 (MQ=255) aaTAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGGCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAAc                                                                                                                                                 <  6:198872/148‑1 (MQ=255)     aGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAc                                                                                                                                            >  5:100264/1‑149 (MQ=255)                 ggTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAgggggg                                                                                                                                 >  1:18188/1‑148 (MQ=255)                  gTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGaa                                                                                                                               >  2:13845/1‑149 (MQ=255)                     ccAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAg                                                                                                                            <  1:13845/149‑1 (MQ=255)                     ccAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAg                                                                                                                            <  2:229276/149‑1 (MQ=255)                        atatGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgt                                                                                                                                        >  6:188480/1‑133 (MQ=255)                        atatGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgt                                                                                                                                        <  5:188480/134‑2 (MQ=255)                                  ctctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACt                                                                                                               >  1:248202/1‑149 (MQ=255)                                     ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGa                                                                                                            >  2:248517/1‑149 (MQ=255)                                       cAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGa                                                                                                            <  1:317628/147‑1 (MQ=255)                                              gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTgg                                                                                                   <  1:248517/149‑1 (MQ=255)                                                          gAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGCGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCt                                                                                        >  7:260861/1‑148 (MQ=255)                                                                  tcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCATGTTGGGTTTCCCCCCTTCAACATCg                                                                               >  6:58813/1‑149 (MQ=255)                                                                   cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                               <  6:100264/148‑1 (MQ=255)                                                                     tGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTc                                                                             >  5:139686/1‑148 (MQ=255)                                                                       tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgtg                                                                                                                                       >  2:235319/1‑86 (MQ=255)                                                                       tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCAACTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgtg                                                                                                                                       <  1:235319/88‑3 (MQ=255)                                                                          aCCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCcaccac                                                                       <  2:18188/149‑1 (MQ=255)                                                                             aGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTgg                                                                                                   >  3:208248/1‑118 (MQ=255)                                                                             aGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTgg                                                                                                   <  4:208248/118‑1 (MQ=255)                                                                              gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCtt                                                                   <  2:248202/149‑1 (MQ=255)                                                                              gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCtt                                                                   <  2:296729/149‑1 (MQ=255)                                                                               tACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTcc                                                                            >  2:241712/1‑139 (MQ=255)                                                                               tACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTcc                                                                            <  1:241712/139‑1 (MQ=255)                                                                                     aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACgg                                                            >  1:241624/1‑149 (MQ=255)                                                                                         cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgg                                                        <  2:241624/149‑1 (MQ=255)                                                                                           gTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCATGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgggg                                                      <  5:58813/149‑1 (MQ=255)                                                                                                   aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                     >  5:316703/1‑142 (MQ=255)                                                                                                   aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                     <  6:316703/142‑1 (MQ=255)                                                                                                    aaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACaa                                                >  7:55693/1‑146 (MQ=255)                                                                                                            ccTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTcc                                                                                             >  1:279632/1‑93 (MQ=255)                                                                                                            ccTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTcc                                                                                             <  2:279632/93‑1 (MQ=255)                                                                                                             cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                    >  3:216588/1‑149 (MQ=255)                                                                                                              tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTc                                   >  8:69490/1‑149 (MQ=255)                                                                                                                  cTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGt                               >  4:199973/1‑149 (MQ=255)                                                                                                                     tttCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGtt                              >  2:238609/1‑147 (MQ=255)                                                                                                                     tttCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCa                            <  8:260861/149‑1 (MQ=255)                                                                                                                          gggTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCCTCAACATCGTCCACCACGCTTGGTAGGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTGAGGGGt                       >  8:65307/1‑149 (MQ=255)                                                                                                                           ggTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGtt                      <  1:179025/149‑1 (MQ=255)                                                                                                                                gAATTTTCCCACGGGGGGAACTCCACCCGGAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                         <  2:154029/141‑1 (MQ=255)                                                                                                                                gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGGTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGtgatcg                  >  1:154029/1‑143 (MQ=255)                                                                                                                                gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                     >  4:302113/1‑145 (MQ=255)                                                                                                                                gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                     <  3:302113/145‑1 (MQ=255)                                                                                                                                      ccctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg           >  3:45669/6‑149 (MQ=255)                                                                                                                                        ctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCCCGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGTGTTCTGCGGt         >  6:243863/4‑149 (MQ=255)                                                                                                                                        ccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTg             >  6:155907/1‑145 (MQ=255)                                                                                                                                        ccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTg             <  5:155907/145‑1 (MQ=255)                                                                                                                                        ccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt         >  7:36885/1‑149 (MQ=255)                                                                                                                                         tgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt         >  8:260266/3‑148 (MQ=255)                                                                                                                                         cACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt         >  8:133837/1‑148 (MQ=255)                                                                                                                                              gTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagagag                                           <  4:312104/109‑1 (MQ=255)                                                                                                                                              gTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagagag                                           >  3:312104/1‑109 (MQ=255)                                                                                                                                               tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGa  <  6:139686/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGA  >  NC_002947/2104551‑2104841

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