Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,181,152 A→T 10.2% 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,1520AT10.2% 130.6 / 6.1 49L264Q (CTG→CAG) PP_4560ribonuclease 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 = 7.97e-01

TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCC  >  NC_002947/5181006‑5181300
                                                                                                                                                   |                                                                                                                                                    
tATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGagcagcag                                                                                                                                                     >  8:57022/1‑148 (MQ=255)
      gCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACgg                                                                                                                                                >  7:71195/1‑147 (MQ=255)
        ttGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGcgc                                                                                                                                             <  4:241789/149‑3 (MQ=255)
              ggggTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTCCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                       <  2:269255/148‑1 (MQ=255)
              gggTTTGGGGGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                       <  8:167459/148‑1 (MQ=255)
                 tttGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGCGCCGAAatgaa                                                                                                                                   <  6:183091/149‑4 (MQ=38)
                 tttGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAAtata                                                                                                                                   <  5:45333/149‑1 (MQ=255)
                    gggTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGa                                                                                                                                 <  1:160793/148‑1 (MQ=255)
                     ggTGCTCGTGACCCAGTACAGTCAACGTATCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGaa                                                                                                                                <  5:249118/148‑1 (MQ=255)
                       tGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAAt                                                                                                                              >  5:186709/1‑148 (MQ=255)
                                cccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                                    >  5:44787/1‑149 (MQ=255)
                                 ccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacag                                                                                                                    >  4:167812/1‑148 (MQ=255)
                                           cAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                          <  7:57022/148‑1 (MQ=255)
                                             aGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCa                                                                                                       >  8:219316/1‑149 (MQ=255)
                                              ggTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTGGATCACCGCATTCATTTCGGCGCAGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                      >  5:114727/1‑149 (MQ=255)
                                                                cccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATg                                                                                     <  6:186709/148‑1 (MQ=255)
                                                                 ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                          <  5:15320/126‑1 (MQ=255)
                                                                 ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGAtt                                                                                                          >  6:15320/1‑126 (MQ=255)
                                                                 ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGc                                                                                   >  1:169086/1‑149 (MQ=255)
                                                                      gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGCGCCGaaa                                                                                                                                       <  3:7191/92‑1 (MQ=25)
                                                                      gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGCGCCGaaa                                                                                                                                       >  4:7191/1‑92 (MQ=25)
                                                                          tCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTc                                                                          <  7:219316/149‑1 (MQ=255)
                                                                            tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATTCTGCCGTACATGGCGTTGTAGTCg                                                                         >  3:38003/1‑148 (MQ=255)
                                                                              tctcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCgg                                                                        <  2:301670/147‑1 (MQ=255)
                                                                                tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACAt                                                                                      <  2:185643/131‑1 (MQ=255)
                                                                                tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACAt                                                                                      >  1:185643/1‑131 (MQ=255)
                                                                                  gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGGTGCCGTAGATGGCGTGGTAGTCGGCTAAgg                                                                  >  4:26224/1‑149 (MQ=255)
                                                                                  gCCCGGGTTCTTGCCTTCGGCTGACAGGTGCTCGCTCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCGGGCGAAATATAGAAATACTGCAACAGACGGATGATTGCACCAATGCAGCGGTAAAGGGCGCTGGAGGCGGCAAAgg                                                                  >  6:155928/1‑149 (MQ=255)
                                                                                            ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                                        >  4:22109/1‑149 (MQ=255)
                                                                                              gCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtag                                                       <  3:167812/148‑1 (MQ=255)
                                                                                                     gCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtagtag                                                    <  7:279485/144‑1 (MQ=255)
                                                                                                     gCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtagtag                                                    >  8:279485/1‑144 (MQ=255)
                                                                                                     gCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCaaa                                               >  6:228168/1‑149 (MQ=255)
                                                                                                       tGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagt                                                      >  4:40119/1‑140 (MQ=255)
                                                                                                       tGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagt                                                      <  3:40119/140‑1 (MQ=255)
                                                                                                         aCATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGccc                                           <  3:119019/149‑1 (MQ=255)
                                                                                                                       acctgattcattTCGGCGCCGAGCAGCAGCACGTCAGCCGAAATATAGAAATAGAGCAAGAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTAGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCATGGAGGCCACGa                              <  4:38003/143‑1 (MQ=255)
                                                                                                                       aCCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGa                              <  3:22109/148‑1 (MQ=255)
                                                                                                                        ccGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATc                            >  5:178990/1‑149 (MQ=255)
                                                                                                                             ttcattTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac                       >  8:53715/1‑149 (MQ=255)
                                                                                                                                attTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCt                                                               >  2:190646/1‑106 (MQ=255)
                                                                                                                                attTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCt                                                               <  1:190646/106‑1 (MQ=255)
                                                                                                                                attTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacacca                     >  2:291722/1‑148 (MQ=255)
                                                                                                                                   tCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACCGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGGCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACTATCCACACCACaa                  <  5:155928/148‑1 (MQ=255)
                                                                                                                                    cGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGATAAGCCCAGGGAGGCCACGATCCACACCACAAc                 <  8:71195/148‑1 (MQ=255)
                                                                                                                                      gctgcCGTGCTGCTGCTCGGCGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc              >  8:170349/4‑149 (MQ=255)
                                                                                                                                      gcgcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc              >  6:274472/1‑149 (MQ=255)
                                                                                                                                      gcgcCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGc              >  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)
                                                                                                                                                   |                                                                                                                                                    
TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCGCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCC  >  NC_002947/5181006‑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: