Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,160	Δ1 bp	7.5%	coding (783/957 nt)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,160	0	A	.	7.5%	158.8 / 7.4	53	coding (783/957 nt)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base . (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

GGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCGGGGTGA  >  NC_002947/5181020‑5181308                                                                                                                                              |                                                                                                                                                     ggggTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTCCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                               <  2:269255/148‑1 (MQ=255) gggTTTGGGGGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                               <  8:167459/148‑1 (MQ=255)    tttGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGAAatgaa                                                                                                                                           <  6:183091/149‑4 (MQ=38)    tttGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAAtata                                                                                                                                           <  5:45333/149‑1 (MQ=255)       gggTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGa                                                                                                                                         <  1:160793/148‑1 (MQ=255)        ggTGCTCGTGACCCAGTACAGTCAACGTATCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGaa                                                                                                                                        <  5:249118/148‑1 (MQ=255)          tGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAAt                                                                                                                                      >  5:186709/1‑148 (MQ=255)                   cccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                                            >  5:44787/1‑149 (MQ=255)                    ccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                                            >  4:167812/1‑148 (MQ=255)                              cAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                                  <  7:57022/148‑1 (MQ=255)                                aGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCa                                                                                                               >  8:219316/1‑149 (MQ=255)                                 ggTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTGGATCACCGCATTCATTTCGGC‑GCAGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                              >  5:114727/1‑149 (MQ=255)                                                   cccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATg                                                                                             <  6:186709/148‑1 (MQ=255)                                                    ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                                  >  6:15320/1‑126 (MQ=255)                                                    ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                                  <  5:15320/126‑1 (MQ=255)                                                    ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGc                                                                                           >  1:169086/1‑149 (MQ=255)                                                         gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGaaa                                                                                                                                               <  3:7191/92‑1 (MQ=25)                                                         gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGaaa                                                                                                                                               >  4:7191/1‑92 (MQ=25)                                                             tCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTc                                                                                  <  7:219316/149‑1 (MQ=255)                                                               tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATTCTGCCGTACATGGCGTTGTAGTCg                                                                                 >  3:38003/1‑148 (MQ=255)                                                                 tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCgg                                                                                <  2:301670/147‑1 (MQ=255)                                                                   tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACAt                                                                                              >  1:185643/1‑131 (MQ=255)                                                                   tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACAt                                                                                              <  2:185643/131‑1 (MQ=255)                                                                     gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGGTGCCGTAGATGGCGTGGTAGTCGGCTAAgg                                                                          >  4:26224/1‑149 (MQ=255)                                                                     gCCCGGGTTCTTGCCTTCGGCTGACAGGTGCTCGCTCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCGGGCGAAATATAGAAATACTGCAACAGACGGATGATTGCACCAATGCAGCGGTAAAGGGCGCTGGAGGCGGCAAAgg                                                                          >  6:155928/1‑149 (MQ=255)                                                                               ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                                                >  4:22109/1‑149 (MQ=255)                                                                                 gCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtag                                                               <  3:167812/148‑1 (MQ=255)                                                                                        gCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtagtag                                                            >  8:279485/1‑144 (MQ=255)                                                                                        gCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtagtag                                                            <  7:279485/144‑1 (MQ=255)                                                                                        gCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCaaa                                                       >  6:228168/1‑149 (MQ=255)                                                                                          tGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagt                                                              <  3:40119/140‑1 (MQ=255)                                                                                          tGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagt                                                              >  4:40119/1‑140 (MQ=255)                                                                                            aCATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGccc                                                   <  3:119019/149‑1 (MQ=255)                                                                                                          acctgattcattTCGGC‑GCCGAGCAGCAGCACGTCAGCCGAAATATAGAAATAGAGCAAGAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTAGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCATGGAGGCCACGa                                      <  4:38003/143‑1 (MQ=255)                                                                                                          aCCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGa                                      <  3:22109/148‑1 (MQ=255)                                                                                                           ccGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATc                                    >  5:178990/1‑149 (MQ=255)                                                                                                                ttcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac                               >  8:53715/1‑149 (MQ=255)                                                                                                                   attTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCt                                                                       <  1:190646/106‑1 (MQ=255)                                                                                                                   attTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCt                                                                       >  2:190646/1‑106 (MQ=255)                                                                                                                   attTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacacca                             >  2:291722/1‑148 (MQ=255)                                                                                                                      tCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACCGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGGCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACTATCCACACCACaa                          <  5:155928/148‑1 (MQ=255)                                                                                                                       cGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGATAAGCCCAGGGAGGCCACGATCCACACCACAAc                         <  8:71195/148‑1 (MQ=255)                                                                                                                         gctgcCGTGCTGCTGCTCGGC‑GCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc                      >  8:170349/4‑149 (MQ=255)                                                                                                                         gc‑gcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc                      >  6:274472/1‑149 (MQ=255)                                                                                                                         gc‑gcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc                      >  8:127895/1‑149 (MQ=255)                                                                                                                            gcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc                      >  3:159930/1‑147 (MQ=255)                                                                                                                            gcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCa                    <  5:228168/149‑1 (MQ=255)                                                                                                                                 gcagcagcaCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCAcc               <  6:178990/149‑1 (MQ=255)                                                                                                                                      agcaCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACGACCGCCAGCACCGAGcc          >  6:120178/1‑149 (MQ=255)                                                                                                                                        caCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGa                                      <  7:35267/119‑1 (MQ=255)                                                                                                                                        caCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGa                                      >  8:35267/1‑119 (MQ=255)                                                                                                                                           ggCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACGACCGCCAGCACCGAGCCCgggg     <  5:120178/149‑1 (MQ=255)                                                                                                                                           ggCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCgggg     <  4:159930/149‑1 (MQ=255)                                                                                                                                              aGCCGAAATATAGAAATACAGCAACAGCACGAGGATTGCACCAATGCGGCCGTAGATGGCGTTGTAGTCGGCAAAGGTCTTATCTTATTATGAATATCCCTGGGAGTCCTAGATCTACACCACAACCGCCAGCACCGAGCCCGGGGtga  <  6:114727/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     GGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCGGGGTGA  >  NC_002947/5181020‑5181308

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