Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,181,782 C→T 7.8% G54D (GGT→GAT)  PP_4560 ← ribonuclease BN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,181,7820CT 7.8% 146.6 / 2.1 51G54D (GGT→GAT) PP_4560ribonuclease BN
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base T (2/2);  total (25/26)
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.41e-01

TCACGATACCCACCGACAGCAGCCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGCACGCGGTG  >  NC_002947/5181639‑5181909
                                                                                                                                               |                                                                                                                               
tCACGATACCCACCGACAGCAGCCCACCCTTGGATTGTTGCAGCTGGTTGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCTAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATc                                                                                                                            <  1:176505/149‑1 (MQ=255)
tCACGATACCCACCGACAGCAGCCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGCCGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGAt                                                                                                                             >  5:5222/1‑148 (MQ=255)
         ccACCGACAGCAGCCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAAt                                                                                                                    >  1:16250/1‑148 (MQ=255)
            ccGACAGCAGCCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAgg                                                                                                                <  5:16931/149‑1 (MQ=255)
              gACAGCAGCCAACCCTTGGATTGTTGCAGCTTGTCGATCACCGGGTTCACCTGATCCAGCGCCTGGTGCGGCAACACCAGTTCCGATTTCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGaa                                                                                                              <  5:267138/149‑1 (MQ=255)
                   cagcCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGaacaa                                                                                                           >  8:233513/1‑147 (MQ=255)
                   cagcCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAacaaca                                                                                                         <  7:60177/149‑1 (MQ=255)
                      cccACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAgg                                                                                                       >  4:3181/1‑148 (MQ=255)
                            cTTGGATTGTTGCAGCTGGTCGATCACCTGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGa                                                                                                <  3:151258/149‑1 (MQ=255)
                               ggATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCATTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGTAAGGGGaa                                                                                               <  7:186461/147‑1 (MQ=255)
                                  ttgttgCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCg                                                                                          <  3:3181/149‑1 (MQ=255)
                                    gttgCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCga                                                                                         <  7:233513/148‑1 (MQ=255)
                                       gCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGaa                                                                                                     <  8:69205/133‑1 (MQ=255)
                                       gCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGaa                                                                                                     >  7:69205/1‑133 (MQ=255)
                                       gCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGaa                                                                                      >  4:60629/1‑148 (MQ=255)
                                              gTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCt                                                                              >  7:186092/1‑149 (MQ=255)
                                                     aCCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGc                                                                       >  6:13328/1‑149 (MQ=255)
                                                     aCCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGc                                                                       >  7:219341/1‑149 (MQ=255)
                                                         ggTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAgc                                                                   <  7:203758/149‑1 (MQ=255)
                                                              aCCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAgcgcg                                                                >  4:26272/1‑146 (MQ=255)
                                                               ccTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAg                                                             >  3:72819/1‑149 (MQ=255)
                                                                    tCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAATCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAg                                                                         <  6:167697/132‑1 (MQ=255)
                                                                    tCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAATCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAg                                                                         <  6:163788/132‑1 (MQ=255)
                                                                    tCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAATCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAg                                                                         >  5:167697/1‑132 (MQ=255)
                                                                    tCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAATCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAg                                                                         >  5:163788/1‑132 (MQ=255)
                                                                    tCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCa                                                          >  6:148524/1‑147 (MQ=255)
                                                                      cAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATa                                                        <  1:159827/147‑1 (MQ=255)
                                                                       aGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGt                                                     >  1:66624/1‑149 (MQ=255)
                                                                          gcCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGa                                                  >  6:45082/1‑149 (MQ=255)
                                                                              ggggCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACAtc                                              >  4:59789/1‑149 (MQ=255)
                                                                               gggCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAg                                                             <  7:224602/133‑1 (MQ=255)
                                                                               gggCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAg                                                             >  8:224602/1‑133 (MQ=255)
                                                                                ggcggcAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACAtctc                                            <  3:60629/149‑1 (MQ=255)
                                                                                 gcggcAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCtcg                                           >  3:84203/1‑149 (MQ=255)
                                                                                   ggcAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCtcgtc                                         <  8:219341/149‑1 (MQ=255)
                                                                                   ggcAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGAGCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGCGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCtcgtc                                         <  1:199572/149‑1 (MQ=255)
                                                                                     cAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAATGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCtcgtcgt                                       >  3:119794/1‑149 (MQ=255)
                                                                                      aaCACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCGGGTAGGCCAGCGCTGAGGCATAGGTGGACATCtcgtcgtc                                      >  5:267947/1‑149 (MQ=255)
                                                                                        cacCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGa                                    <  2:148612/149‑1 (MQ=255)
                                                                                                 cGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTcc                            <  4:84203/148‑1 (MQ=255)
                                                                                                       cAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCAAGGAATTCCTCGACga                     <  5:45082/149‑2 (MQ=255)
                                                                                                        aGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGt                    >  1:156405/1‑149 (MQ=255)
                                                                                                              aGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAAc              <  6:5222/149‑1 (MQ=255)
                                                                                                              aGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAAc              <  5:13328/149‑1 (MQ=255)
                                                                                                               gCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGATTAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAAc              >  8:53932/1‑148 (MQ=255)
                                                                                                                 cAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAg           >  5:139336/1‑149 (MQ=255)
                                                                                                                  aGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAATGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGc          <  4:119794/149‑1 (MQ=255)
                                                                                                                  aGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGc          <  5:148524/149‑1 (MQ=255)
                                                                                                                   ggAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGCa         >  8:184702/1‑149 (MQ=255)
                                                                                                                    gagaAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGCAc        <  2:66624/149‑1 (MQ=255)
                                                                                                                          aaGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGCACGCGGTg  <  4:72819/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                               
TCACGATACCCACCGACAGCAGCCCACCCTTGGATTGTTGCAGCTGGTCGATCACCGGGTTCACCTGCTCCAGCGCCTGGGGCGGCAACACCAGTTCCGATTGCAGACGCAGCCAGGAGAAAAAGTCCGGCAAATGTAGGAAACCGATCAGGGCAATCAGGAACAACAGGAAGGGGAATAGCGAGAACAGCGCCTGGTAGGCCAGCGCTGAGGCATAGGTGGACATCTCGTCGTCGAGGAATTCCTTGACGGTGCGAACCAGCACGCGGTG  >  NC_002947/5181639‑5181909

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