Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,155	A→T	10.2%	V263E (GTG→GAG)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,181,155	0	A	T	10.2%	125.1 / 6.1	49	V263E (GTG→GAG)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base T (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCC  >  NC_002947/5181012‑5181300                                                                                                                                                 |                                                                                                                                                  gCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACgg                                                                                                                                                >  7:71195/1‑147 (MQ=255)   ttGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGcgc                                                                                                                                             <  4:241789/149‑3 (MQ=255)         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                                                                                                                                       >  4:7191/1‑92 (MQ=25)                                                                 gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGaaa                                                                                                                                       <  3:7191/92‑1 (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                                                    <  7:279485/144‑1 (MQ=255)                                                                                                gCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtagtag                                                    >  8:279485/1‑144 (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                                                               >  2:190646/1‑106 (MQ=255)                                                                                                                           attTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCt                                                               <  1:190646/106‑1 (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)                                                                                                                                                 |                                                                                                                                                  GCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCC  >  NC_002947/5181012‑5181300

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