Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,449,467 C→T 8.5% R74R (CGG→CGA yhaV ← toxin antitoxin system SohB(PrlF)‑YhaV toxin endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,449,4670CT 8.5% 131.7 / 3.5 47R74R (CGG→CGAyhaVtoxin antitoxin system SohB(PrlF)‑YhaV toxin endoribonuclease
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGTTGAGCCTCGACTTGAGCGCTTAGCGCGTCC  >  NC_002947/1449327‑1449611
                                                                                                                                            |                                                                                                                                                 
gCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCg                                                                                                                                          >  1:734869/1‑150 (MQ=255)
 cccGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCgg                                                                                                                                         <  2:786084/150‑1 (MQ=255)
               tCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCgg                                                                                                                               >  1:110204/1‑146 (MQ=255)
               tCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCgg                                                                                                                               <  2:110204/146‑1 (MQ=255)
               tCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCt                                                                                                                           >  1:506992/1‑150 (MQ=255)
                cGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCtt                                                                                                                          <  2:734869/150‑1 (MQ=255)
                 gTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTg                                                                                                                         >  1:107187/1‑150 (MQ=255)
                 gTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTg                                                                                                                         >  1:693171/1‑150 (MQ=255)
                    aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                                                                                                                      >  2:513134/1‑150 (MQ=255)
                    aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                                                                                                                      <  1:513134/150‑1 (MQ=255)
                    aTTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                                                                                                                      >  2:363887/1‑150 (MQ=255)
                                cAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGaa                                                                                                          <  2:693171/150‑1 (MQ=255)
                                   cACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGc                                                                                                       >  2:436073/1‑150 (MQ=255)
                                                           gTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCt                                                                               >  1:276052/1‑150 (MQ=255)
                                                            tAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCtt                                                                              >  2:754077/1‑150 (MQ=255)
                                                             aGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTTGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCttt                                                                             <  2:506992/150‑1 (MQ=255)
                                                              gCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCtttt                                                                            <  1:436073/150‑1 (MQ=255)
                                                                ggAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTaa                                                                         >  2:366776/1‑150 (MQ=255)
                                                                 gaagaaCAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTaaa                                                                        <  1:754077/150‑1 (MQ=255)
                                                                     aacaacCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGACCACCGACGCATCGCCGCCAGCCGTGTAAAggga                                                                    >  2:171480/1‑148 (MQ=255)
                                                                          ccGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGtt                                                                                                                                      <  1:454098/80‑1 (MQ=255)
                                                                          ccGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGtt                                                                                                                                      >  2:454098/1‑80 (MQ=255)
                                                                               aCTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACAACCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCg                                                          <  1:171480/150‑1 (MQ=255)
                                                                                   ctgGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCTGATGGCGGCCAGCCTTGTAAAGGTATTTCTTTTCGtata                                                      >  1:130828/1‑150 (MQ=255)
                                                                                        agaaCTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                                                                                                                      >  1:356174/1‑82 (MQ=255)
                                                                                        agaaCTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAgg                                                                                                                      <  2:356174/82‑1 (MQ=255)
                                                                                            cTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGTCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATACCCCGTTgg                                             <  2:276052/150‑1 (MQ=255)
                                                                                             tttGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTggg                                            >  2:669847/1‑150 (MQ=255)
                                                                                                  tctGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGt                                           >  1:382300/1‑146 (MQ=255)
                                                                                                  tctGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGt                                           <  2:382300/146‑1 (MQ=255)
                                                                                                  tctGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCttt                                       >  1:780240/1‑150 (MQ=255)
                                                                                                    tGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCt                                         <  1:637823/146‑1 (MQ=255)
                                                                                                    tGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCt                                         >  2:637823/1‑146 (MQ=255)
                                                                                                         cAGTGCTTGTGGTCACCGCCCAGTGTTGACACTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGt                                <  2:130828/150‑1 (MQ=255)
                                                                                                           gTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGtt                               <  1:92998/149‑1 (MQ=255)
                                                                                                           gTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGtt                               >  2:92998/1‑149 (MQ=255)
                                                                                                              cTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTa                                                                    >  1:562673/1‑109 (MQ=255)
                                                                                                              cTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTa                                                                    <  2:562673/109‑1 (MQ=255)
                                                                                                                 gtgGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGTTGAGCCTc                        >  2:479354/1‑150 (MQ=255)
                                                                                                                   ggTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGTTGAGCCTCGa                      <  1:366776/150‑1 (MQ=255)
                                                                                                                        ccGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGTTGAGCCTCGACTTGa                 <  1:363887/150‑1 (MQ=255)
                                                                                                                         cGCCAAGTGTTGCCCCTTGTCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTc                                                                >  1:583264/1‑102 (MQ=255)
                                                                                                                         cGCCAAGTGTTGCCCCTTGTCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTc                                                                <  2:186053/102‑1 (MQ=255)
                                                                                                                         cGCCAAGTGTTGCCCCTTGTCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTc                                                                <  2:583264/102‑1 (MQ=255)
                                                                                                                         cGCCAAGTGTTGCCCCTTGTCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTc                                                                >  1:186053/1‑102 (MQ=255)
                                                                                                                                tgtTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGTTGAGCCTCGACTTGAGCGCTTAg         <  2:791373/150‑1 (MQ=255)
                                                                                                                                       ccTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTCCGCGGGATCGCTTTGGTATAATCCGTTGGGACTTTCTCCCCTCGTGCCCTGGCTCTATTGCTTAGCGCGTcc  >  1:413226/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GCCCGCTTGCTCGACTCGTCATTGATCCAGGCCAGCACGATGATCCTGCTGGGGGTGTGGTAGCGGAAGAACAACCGGTACTGCTGGAAGAACTTTGCTCTGAACCAGTGCTTGTGGTCACCGCCAAGTGTTGCCCCTTGCCGGTATTCGGGTTTGGTCGGGTCTTGAGGGATTACATCGAACGCCAACCGCCGGATGGCGGCCAGCCTTTTAAAGGTATTTCTTTTCGTATAACCCGTTGGGTCTTTGTCCCGTTGAGCCTCGACTTGAGCGCTTAGCGCGTCC  >  NC_002947/1449327‑1449611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: