Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,484,077 C→T 7.4% R315C (CGT→TGT)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,484,0770CT 7.4% 157.6 / 2.5 54R315C (CGT→TGT) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/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 = 8.99e-01

CACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAACCCGTATCTGGTA  >  NC_002947/2483936‑2484223
                                                                                                                                             |                                                                                                                                                  
cACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATgg                                                                                                                                             >  3:339642/1‑149 (MQ=255)
cACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATg                                                                                                                                              <  3:61885/148‑1 (MQ=255)
 aCGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGgc                                                                                                                                            >  1:368823/1‑149 (MQ=255)
      aGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCa                                                                                                                                        <  2:52388/148‑1 (MQ=255)
             acgacgCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAAc                                                                                                                                >  5:276420/1‑149 (MQ=255)
               gacgCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAAcgc                                                                                                                              >  5:122248/1‑149 (MQ=255)
                  gccgccGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACgcggc                                                                                                                           >  1:333105/1‑149 (MQ=255)
                  gccgccGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACgcggc                                                                                                                           <  4:339642/149‑1 (MQ=255)
                     gccgGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCAGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAg                                                                                                                        <  8:171757/149‑1 (MQ=255)
                                cATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCt                                                                                                             <  2:368823/149‑1 (MQ=255)
                                       gCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGTCGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGttt                                                                                                      >  4:356719/1‑149 (MQ=255)
                                       gCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGttt                                                                                                      >  2:79042/1‑149 (MQ=255)
                                       gCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGttt                                                                                                      <  4:90262/149‑1 (MQ=255)
                                           gAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGTCGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGaa                                                                                                  <  3:356719/149‑1 (MQ=255)
                                           gAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGaa                                                                                                  <  1:79042/149‑1 (MQ=255)
                                                    aGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGAc                                                                                         <  6:122248/149‑1 (MQ=255)
                                                     gcgcAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGAcg                                                                                        >  6:49755/1‑149 (MQ=255)
                                                      cgcAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGAcgc                                                                                       >  6:40581/1‑149 (MQ=255)
                                                       gcAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGAcgcg                                                                                      >  3:362765/1‑149 (MQ=255)
                                                        cAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGt                                                                                     <  6:276420/149‑1 (MQ=255)
                                                        cAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGt                                                                                     <  2:350001/149‑1 (MQ=255)
                                                        cAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGt                                                                                     <  2:333105/149‑1 (MQ=255)
                                                         aaGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGATGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTg                                                                                    >  5:259520/1‑149 (MQ=255)
                                                         aaGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGATGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTg                                                                                    >  1:325745/1‑149 (MQ=255)
                                                             tGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTgcatc                                                                                >  2:113853/1‑149 (MQ=255)
                                                              gcgcCATGCGGTGGCAGGCTGCCTGGAAGTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTgcatcc                                                                               >  1:212808/1‑149 (MQ=255)
                                                                      cGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTgg                                                                                                           >  4:359775/1‑113 (MQ=255)
                                                                      cGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTgg                                                                                                           <  3:359775/113‑1 (MQ=255)
                                                                         tGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAg                                                                    >  3:166437/1‑149 (MQ=255)
                                                                             aGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTcc                                                                 >  4:89055/1‑148 (MQ=255)
                                                                                       gAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCgcgc                                                        <  7:29789/147‑1 (MQ=255)
                                                                                        aaCTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTgcatccgcatcc                                                                         >  8:207002/1‑129 (MQ=255)
                                                                                        aaCTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTgcatccgcatcc                                                                         <  7:207002/129‑1 (MQ=255)
                                                                                        aaCTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAAc                                                     <  5:40581/149‑1 (MQ=255)
                                                                                             ctgcCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTgcg                                                <  3:89055/149‑1 (MQ=255)
                                                                                             ctgcCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTg                                                  <  5:49755/147‑1 (MQ=255)
                                                                                                     cTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGATGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGc                                        <  2:325745/149‑1 (MQ=255)
                                                                                                     cTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGATGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGc                                        <  6:259520/149‑1 (MQ=255)
                                                                                                      tGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGaa                                                                                                  <  4:73508/90‑1 (MQ=255)
                                                                                                      tGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGaa                                                                                                  >  3:73508/1‑90 (MQ=255)
                                                                                                                       ccccAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGcc                       >  7:38751/1‑148 (MQ=255)
                                                                                                                           aaCCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGgcgc                        >  7:226517/1‑143 (MQ=255)
                                                                                                                           aaCCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGgcgc                        <  8:226517/143‑1 (MQ=255)
                                                                                                                            aCCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGcc                 <  2:212808/149‑1 (MQ=255)
                                                                                                                                 aaTGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATccc                                                                        <  6:102864/89‑1 (MQ=255)
                                                                                                                                 aaTGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATccc                                                                        >  5:102864/1‑89 (MQ=255)
                                                                                                                                 aaTGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAAccc            >  3:298532/1‑149 (MQ=255)
                                                                                                                                     ccTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGa                                          >  7:120173/1‑115 (MQ=255)
                                                                                                                                     ccTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGa                                          <  8:120173/115‑1 (MQ=255)
                                                                                                                                      cTACTGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAACCCGTATc       >  3:36160/1‑149 (MQ=255)
                                                                                                                                      cTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAACCCGt          >  5:387732/1‑146 (MQ=255)
                                                                                                                                      cTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAACCCGt          <  6:387732/146‑1 (MQ=255)
                                                                                                                                      cTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGCGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACg                   >  3:122028/1‑137 (MQ=255)
                                                                                                                                           gACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTCACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAACCAGTATCTGGTa  >  2:284993/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACGTTATGGCGCCATGGTCAACGCGGCGAGCAAGGCCAGCTGGGGTTTCGAAGACCGTGACGCGTGCATCCGCATCCCCGAGTCCGACCCGCGCAACCTGCGCATCGAGCACCGCCTCGCTGGCGCCGACGCCAACCCGTATCTGGTA  >  NC_002947/2483936‑2484223

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: