Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,181,155	A→T	11.5%	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	11.5%	115.6 / 4.2	43	V263E (GTG→GAG)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base T (3/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCG  >  NC_002947/5181018‑5181296                                                                                                                                           |                                                                                                                                               tCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGa                                                                                                                                     >  6:76383/1‑148 (MQ=255)          ggTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGaa                                                                                                                            <  5:139465/148‑1 (MQ=255)           gTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGaaa                                                                                                                           >  3:262805/1‑148 (MQ=255)                 gTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATacagc                                                                                                                     >  4:148214/1‑148 (MQ=255)                   gACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATacagcaac                                                                                                                  <  5:76383/149‑1 (MQ=255)                       cAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GACGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATacagcaacagca                                                                                                              >  5:165728/1‑149 (MQ=255)                               tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                      >  5:86823/1‑149 (MQ=255)                                       tCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCaa                                                                                               <  6:165728/148‑1 (MQ=255)                                            gggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCAc                                                                                                             <  8:235419/129‑1 (MQ=255)                                            gggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCAc                                                                                                             >  7:235419/1‑129 (MQ=255)                                                ttttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACga                                                                                                           >  8:8660/1‑127 (MQ=255)                                                ttttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACga                                                                                                           <  7:8660/127‑1 (MQ=255)                                                 tttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTa                                                                                     >  5:212656/1‑148 (MQ=255)                                                  ttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTa                                                                                     <  4:262805/147‑1 (MQ=255)                                                    ccccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATg                                                                                 >  7:57231/1‑149 (MQ=255)                                                     cccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCg                                                                                       >  7:150022/1‑142 (MQ=255)                                                     cccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCg                                                                                       <  8:150022/142‑1 (MQ=255)                                                         tCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                            <  3:255190/149‑1 (MQ=255)                                                            taGTCTTTCTCGGGCGGGGTCTTGCGTTCGGCTGACTTGTGCTCGATCACCGCATTCATTTCGGC‑GCGGAGCAGCAGCA‑CGGCAGCGGTAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                      <  5:287811/119‑1 (MQ=255)                                                            aaGTCTTTCTCGCCCCGGATCTTGCCTTCGGCTGACATGTGCTCCATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATACAAAAACAGCCACAGCACGATGAtt                                                                                                      >  6:287811/1‑120 (MQ=255)                                                              gTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGt                                                                       <  7:45631/149‑1 (MQ=255)                                                               tCTTTCTCGCCCGGGTTCTTTCCTTCTGCTGACATGTGCTCGATCACCTCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTc                                                                      <  8:72530/149‑1 (MQ=255)                                                                    ctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGaaa                                                                                                                                   <  3:5540/83‑1 (MQ=255)                                                                    ctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGaaa                                                                                                                                   >  4:5540/1‑83 (MQ=255)                                                                    ctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCa                                                                  <  6:86823/148‑1 (MQ=255)                                                                       gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAgg                                                              >  7:290447/1‑149 (MQ=255)                                                                       gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                               <  6:198266/148‑1 (MQ=255)                                                                           gggTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCtt                                                          <  3:148214/149‑1 (MQ=255)                                                                              ttCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCAc                                                       <  6:212656/149‑1 (MQ=255)                                                                                    ccTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagta                                                 <  3:71167/149‑1 (MQ=255)                                                                                      ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCATCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                  <  1:228338/99‑1 (MQ=255)                                                                                      ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                  >  2:228338/1‑98 (MQ=255)                                                                                      ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAgg                                               >  5:220522/1‑149 (MQ=255)                                                                                      ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGAAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAgg                                               >  5:183728/1‑149 (MQ=255)                                                                                          gCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCa                                             >  2:105674/1‑147 (MQ=255)                                                                                                     cTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAgg                                <  8:57231/149‑1 (MQ=255)                                                                                                                  tatattTCGGCTGCCGTGCTGCTGCT‑CGGC‑GCCGAAATGTAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGCACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac                   >  3:97026/4‑149 (MQ=37)                                                                                                                    tattTCGGCTGCCGTGCTGCTGCT‑CGGC‑GCCGAAATGTAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGCTCATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacacc                  <  4:97026/147‑1 (MQ=37)                                                                                                                      tttCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacaccac                >  8:298478/1‑148 (MQ=255)                                                                                                                          ggctgcCGTGCTGCTGCT‑CGGC‑GCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCg           >  5:145675/5‑149 (MQ=255)                                                                                                                             tgcCGTGCTGCTGCT‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAgg                                    >  2:135500/2‑122 (MQ=255)                                                                                                                             tgcCGTGCTGCTGCT‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAgg                                    <  1:135500/121‑1 (MQ=255)                                                                                                                                cGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCGTCACGTTGTAGGCAAAGCCCAGGGAGGCCCCGATCCACACCACACCCGCCAg       >  8:219083/1‑148 (MQ=255)                                                                                                                                    cagcagca‑CGGCAGCCGAAATATAGAATTACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGGAGTAGGCAAAGCCCAGGGAGGGCAGGTTCCCCCCCACAACCGCCGGCACCg  >  3:239869/1‑149 (MQ=255)                                                                                                                                    cagcagca‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCAcc   >  4:170051/1‑148 (MQ=255)                                                                                                                                           |                                                                                                                                               TCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCG  >  NC_002947/5181018‑5181296

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