Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,181,152 A→T 11.6% L264Q (CTG→CAG)  PP_4560 ← ribonuclease BN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,181,1520AT11.6% 113.7 / 6.0 43L264Q (CTG→CAG) PP_4560ribonuclease 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 = 9.07e-01

TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCG  >  NC_002947/5181006‑5181296
                                                                                                                                                   |                                                                                                                                                 
tATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGagcagcag                                                                                                                                                  >  4:255190/1‑148 (MQ=255)
tATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTCCTCGATCACCGCATTCATTTCGGC‑GCCTagcagcag                                                                                                                                                  >  7:72530/1‑148 (MQ=255)
            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                                                                                                             >  7:235419/1‑129 (MQ=255)
                                                       gggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCAc                                                                                                             <  8:235419/129‑1 (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                                                                                       <  8:150022/142‑1 (MQ=255)
                                                                cccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCg                                                                                       >  7:150022/1‑142 (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                                                                                                                                   >  4:5540/1‑83 (MQ=255)
                                                                               ctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGaaa                                                                                                                                   <  3:5540/83‑1 (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)
                                                                                                                                                   |                                                                                                                                                 
TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCA‑CGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCG  >  NC_002947/5181006‑5181296

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