Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,452	G→T	7.4%	S137* (TCA→TAA) ‡	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,452	0	G	T	7.4%	163.5 / 2.1	54	S137* (TCA→TAA) ‡	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base T (2/2);  total (28/26)
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.99e-01

GTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGT  >  NC_002947/1017306‑1017596                                                                                                                                                   |                                                                                                                                                 gTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgat                                                                                                                                                <  6:88462/149‑1 (MQ=255)       aTTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACACCTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAacac                                                                                                                                          >  6:303752/1‑149 (MQ=255)           tGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCa                                                                                                                                        >  7:258862/1‑147 (MQ=255)             aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCaaa                                                                                                                                      <  8:309186/147‑1 (MQ=255)             aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAg                                                                                                                                     <  7:159378/148‑1 (MQ=255)             aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCaaa                                                                                                                                      <  8:75569/147‑1 (MQ=255)                 ggTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTc                                                                                                                                >  6:92033/1‑149 (MQ=255)                          tGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACAttt                                                                                                                       >  2:48452/1‑149 (MQ=255)                                gagaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTttcttc                                                                                                                  <  1:150058/148‑1 (MQ=255)                                gagaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCg                                                                                                                 >  8:130378/1‑149 (MQ=255)                                   aTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGaa                                                                                                                                             >  7:124907/1‑118 (MQ=255)                                   aTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGaa                                                                                                                                             <  8:124907/118‑1 (MQ=255)                                         ggggCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATaaa                                                                                                        <  6:38900/149‑1 (MQ=255)                                              cGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTa                                                                                                   <  8:258862/149‑1 (MQ=255)                                              cGCTTTGCGGACCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGACTCGTGACATTTTCTTAGGTATATAAATTGTa                                                                                                   >  5:168827/1‑149 (MQ=255)                                                     cGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGAtttt                                                                                             <  1:48452/148‑1 (MQ=255)                                                       gCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTgg                                                                                          >  7:112465/1‑149 (MQ=255)                                                       gCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCaaa                                                                                                                                      >  2:272892/1‑105 (MQ=255)                                                       gCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAACGCAGGGCCCTGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCaaa                                                                                                                                      <  1:272892/105‑1 (MQ=255)                                                         ccTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTgggg                                                                                        >  5:321127/1‑149 (MQ=255)                                                             ttGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGa                                                                                      <  5:92033/147‑1 (MQ=255)                                                                cTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgtgt                                                                                  <  3:53445/148‑1 (MQ=255)                                                                                cAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAtttt                                                                  <  6:321127/148‑1 (MQ=255)                                                                                cAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAAAAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAttttt                                                                 >  5:14716/1‑149 (MQ=255)                                                                                          gtgCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                       >  1:223289/1‑149 (MQ=255)                                                                                           tgCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                       <  2:223289/148‑1 (MQ=255)                                                                                                ggCAGGGCCCGGATGGCCGGCTTGAGATCAATAAATCCAGGTTGGAGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGCGGAGTGTGCCAATGGCGCATTTTTCGCTGCGGGCCTGTAg                                                 <  8:112465/149‑1 (MQ=255)                                                                                                  cAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGat                                               >  2:58739/1‑149 (MQ=255)                                                                                                   aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATATATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                       <  4:101744/140‑1 (MQ=255)                                                                                                   aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATATATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                       >  3:101744/1‑140 (MQ=255)                                                                                                   aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGATTGTGACAATGGCGTATTTTTCGCTGCGGTCCTGTAGata                                              >  3:39831/1‑149 (MQ=255)                                                                                                     ggCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGGCAATGGCGTATTTTTCGCTGCGGGCCTGTAGatat                                             >  8:77919/1‑148 (MQ=255)                                                                                                             tGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCa                                    <  7:130378/149‑1 (MQ=255)                                                                                                                tgtgCTTTAGATTAATAAATTCAGTTTGGTGTTCTGATTAACACGAAAGGCCTAGTGACATTTTGTTCGGTATATTAATTTTATATTTTTGGGGGAGTGTTCCAATGTCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTAATa                                  <  6:168827/145‑1 (MQ=255)                                                                                                                          aTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTa                       <  5:303752/149‑1 (MQ=255)                                                                                                                          aTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTa                       >  2:18712/1‑149 (MQ=255)                                                                                                                            cAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTAt                      >  5:277132/1‑148 (MQ=255)                                                                                                                             aaTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATc                     >  3:329991/1‑148 (MQ=255)                                                                                                                              aTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCa                    >  8:103320/1‑148 (MQ=255)                                                                                                                               tAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGTTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCaaa                  <  1:58739/149‑1 (MQ=255)                                                                                                                               aaaaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCaaa                  <  6:14716/148‑1 (MQ=255)                                                                                                                                aaaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAat                 >  2:130883/1‑149 (MQ=255)                                                                                                                                 aaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCAt                                   >  7:218472/1‑130 (MQ=255)                                                                                                                                 aaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCAt                                   <  8:218472/130‑1 (MQ=255)                                                                                                                                  aTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAAt                             <  7:75916/135‑1 (MQ=255)                                                                                                                                  aTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAAt                             >  8:75916/1‑135 (MQ=255)                                                                                                                                  aTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatc                >  3:54344/1‑148 (MQ=255)                                                                                                                                      aGTTTGGTGTTCTGATGAACAACATAGGCCTCGTTAGATTTTCTTAGTTATATAAATTGTAGATTTTTGTGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatctatct           <  4:39831/149‑1 (MQ=255)                                                                                                                                        tttGGTGTTCATCAGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTg         >  8:194697/1‑149 (MQ=255)                                                                                                                                            gtgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatct               >  4:152416/1‑139 (MQ=255)                                                                                                                                            gtgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatct               <  3:152416/139‑1 (MQ=255)                                                                                                                                             tgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCaaaa    <  4:54344/149‑1 (MQ=255)                                                                                                                                                tCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGt  >  2:81410/1‑148 (MQ=255)                                                                                                                                                 cTGATGAACACCAAAGGCCTCGTGACA‑TTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCgg                                                         <  3:53220/91‑1 (MQ=255)                                                                                                                                                 cTGATGAACACCAAAGGCCTCGTGACA‑TTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCgg                                                         >  4:53220/1‑91 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGT  >  NC_002947/1017306‑1017596

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