Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,857,521	Δ1 bp	9.1%	coding (415/558 nt)	PP_4272 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,857,521	0	A	.	9.1%	167.2 / 11.1	55	coding (415/558 nt)	PP_4272	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base . (2/3);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCACCCATCCCCACGCCAAC  >  NC_002947/4857377‑4857666                                                                                                                                                  |                                                                                                                                                  gTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCTAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAgc                                                                                                                                               <  5:294755/149‑1 (MQ=255) gTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTTGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAgc                                                                                                                                               <  2:170079/149‑1 (MQ=255)    acacTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGTCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAgcggc                                                                                                                                            <  8:71500/149‑1 (MQ=255)      acTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTAGGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                                          <  8:101883/149‑5 (MQ=255)      acTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGgcgc                                                                                                                                          >  5:169733/1‑149 (MQ=255)                gTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGt                                                                                                                                 <  7:198051/148‑1 (MQ=255)                gTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGt                                                                                                                                 >  8:198051/1‑148 (MQ=255)                      ccGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑CCTGCATGAAGCACAGCGCCGCCACGGCTGTacccaca                                                                                                                          >  6:89392/1‑149 (MQ=255)                       cGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacac                                                                                                                         <  6:169733/149‑1 (MQ=255)                            tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                                          >  1:66906/1‑123 (MQ=255)                            tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgatgg                                                                                                                                          <  2:66906/127‑5 (MQ=255)                            tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgat                                                                                                                                            >  7:191381/1‑123 (MQ=255)                            tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGTGCTTCATGCAGC‑CAgcgat                                                                                                                                            <  8:191381/125‑3 (MQ=255)                            tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacacccac                                                                                                                     >  8:12830/1‑148 (MQ=255)                            tCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTacccacacccac                                                                                                                     >  5:146449/1‑148 (MQ=255)                                gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCAc                                                                                                                                       <  6:115831/126‑1 (MQ=255)                                gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCAc                                                                                                                                       >  5:115831/1‑126 (MQ=255)                                gAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAacta                                                                                                                <  1:157557/149‑1 (MQ=255)                                      ccccACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAg                                                                                                          <  8:145234/149‑1 (MQ=255)                                          acacCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTaa                                                                                                      <  7:266652/149‑1 (MQ=255)                                                 acacAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAagag                                                                                                 >  7:167119/1‑147 (MQ=255)                                                 acacAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTa                                                                                               <  8:167119/149‑1 (MQ=255)                                                      aaaTGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGa                                                                                          >  7:83555/1‑149 (MQ=255)                                                       aaTGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACgg                                                                                                                                     >  4:185113/1‑105 (MQ=255)                                                             tGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAc                                                                                   <  6:146449/149‑1 (MQ=255)                                                             tGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAc                                                                                   >  2:161266/1‑149 (MQ=255)                                                             tGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAc                                                                                   <  6:271665/149‑1 (MQ=255)                                                               gCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAAcc                                                                                  <  2:291858/148‑1 (MQ=255)                                                               gCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGAAGGTAAAAGAGTAGGCGCCCCCAACCa                                                                                 >  2:305861/1‑149 (MQ=255)                                                                ccAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACgg                                                                                                                                     >  7:254132/1‑96 (MQ=255)                                                                ccAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACgg                                                                                                                                     <  8:254132/96‑1 (MQ=255)                                                                       ctTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTccc                                                                         >  7:83840/1‑149 (MQ=255)                                                                           aGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTAc                                                                     >  2:90099/1‑149 (MQ=255)                                                                           aGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTAc                                                                     <  8:83555/149‑1 (MQ=255)                                                                               aGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATc                                                                 <  7:49865/149‑1 (MQ=255)                                                                                gTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATca                                                                <  8:83840/149‑1 (MQ=255)                                                                                    cgcgCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACaa                                                             <  1:161266/148‑1 (MQ=255)                                                                                                ggcgtcCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATcgc                                                >  2:74877/1‑149 (MQ=255)                                                                                                    tcCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGaaa                                            <  7:12830/149‑1 (MQ=255)                                                                                                       tGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGCCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCAACGATCACAACCAGCAGCATCGCGaaaacc                                         >  5:245851/1‑147 (MQ=255)                                                                                                            gAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGcc                                    >  1:339110/1‑149 (MQ=255)                                                                                                               aaGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACCACTACTGCAGGTAAAAGCGTAGGCGACCGCACCCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGt                                 >  6:367853/1‑149 (MQ=255)                                                                                                                aGAACCCTGCCCCCATCGCTG‑GCTGCATGAACCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAATGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTg                                <  5:89392/149‑1 (MQ=255)                                                                                                                    cccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTaaaa                                                                                                    <  3:4080/77‑1 (MQ=255)                                                                                                                    cccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTaaaa                                                                                                    >  4:4080/1‑77 (MQ=255)                                                                                                                     ccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTccc                                                                         >  7:234932/1‑103 (MQ=255)                                                                                                                     ccTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTccc                                                                         <  8:234932/103‑1 (MQ=255)                                                                                                                      cTGCCCCCATCGCTG‑CCTGCATGAAGCACAGCGGCGCCACGGCTTTACCCACACCCACAAGTACTGCAGGTAAAAGAGTAGGCTACCGCAACAAGCCCTCCCCTCCGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAg                          <  1:74877/149‑1 (MQ=255)                                                                                                                       tGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCACCCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGt                         >  7:165582/1‑149 (MQ=255)                                                                                                                           cccATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAAt                     >  6:38575/1‑149 (MQ=255)                                                                                                                           cccATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAAt                     >  6:284228/1‑149 (MQ=255)                                                                                                                                            tGAAGCACAGCCGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGAAAAAGAGTAGGCGACCGCAACCAGCCCACCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACCCCGCTGGCGTAGTGAATCACCCATACCCACGc      >  5:203489/1‑148 (MQ=255)                                                                                                                                               aGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCACCCATCCCCACGCCAAc  >  2:96743/1‑149 (MQ=255)                                                                                                                                                gCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCAc                  <  5:206692/132‑1 (MQ=255)                                                                                                                                                gCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCAc                  >  6:206692/1‑132 (MQ=255)                                                                                                                                                 cacaGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATca                                                                <  2:198941/85‑1 (MQ=255)                                                                                                                                                 cacaGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATca                                                                >  1:198941/1‑85 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTG‑GCTGCATGAAGCACAGCGGCGCCACGGCTGTACCCACACCCACAACTACTGCAGGTAAAAGAGTAGGCGACCGCAACCAGCCCTCCCCTACGATCACAACCAGCAGCATCGCGAAAAGCACGCCGGTGGCGTAGTGAATCACCCATCCCCACGCCAAC  >  NC_002947/4857377‑4857666

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