Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,774,420 A→C 6.7% S535R (AGT→CGT) ‡ PP_1581 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,774,4200AC 6.7% 189.1 / 4.8 60S535R (AGT→CGT) ‡PP_1581BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (23/33);  new base C (2/2);  total (25/35)
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.50e-01

GAGCATATTTCCCAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAA  >  NC_002947/1774279‑1774556
                                                                                                                                             |                                                                                                                                        
gAGCATATTTCCCAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGa                                                                                                                                   >  1:41575/1‑149 (MQ=255)
 aGCATATTTCCCAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGa                                                                                                                                   >  4:134935/1‑148 (MQ=255)
      atTTCCCAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACatgat                                                                                                                             >  5:195836/1‑149 (MQ=255)
        ttacaaGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGAAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGTAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTAGCGACatgatg                                                                                                                            >  8:63240/6‑148 (MQ=255)
         tCCCAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACAtgatga                                                                                                                           <  2:21258/148‑1 (MQ=255)
          cccAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGaa                                                                                                                          <  8:253320/148‑1 (MQ=255)
           ccAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGa                                                                                                                        >  5:191303/1‑149 (MQ=255)
           ccAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGa                                                                                                                        >  3:298612/1‑149 (MQ=255)
                ttCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACat                                                                                                                                <  1:48317/136‑1 (MQ=255)
                ttCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACat                                                                                                                                >  2:48317/1‑136 (MQ=255)
                    gCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGc                                                                                                               >  5:68501/1‑149 (MQ=255)
                       caacaaCGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCtt                                                                                                            <  2:41575/149‑1 (MQ=255)
                        accaacGGATGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGCCGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTg                                                                                                           <  6:310246/147‑1 (MQ=255)
                        aacaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCt                                                                                                                     >  4:83098/1‑139 (MQ=255)
                        aacaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCt                                                                                                                     <  3:83098/139‑1 (MQ=255)
                        aacaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTg                                                                                                           >  6:164816/1‑149 (MQ=255)
                         acaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCTAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTgg                                                                                                          >  3:282258/1‑149 (MQ=255)
                         acaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTgg                                                                                                          <  3:134935/149‑1 (MQ=255)
                         acaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTgg                                                                                                          <  4:102612/149‑1 (MQ=255)
                         acaacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTgg                                                                                                          <  4:229614/149‑1 (MQ=255)
                          caacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTggg                                                                                                         <  6:68501/149‑1 (MQ=255)
                           aacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGCACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGCCCTGGACGCCTTGGGt                                                                                                        >  7:110590/1‑149 (MQ=255)
                           aacGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGt                                                                                                        <  4:298612/149‑1 (MQ=255)
                              ggAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGaa                                                                                                                          <  5:74406/128‑1 (MQ=255)
                              ggAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGACGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGaa                                                                                                                          >  6:74406/1‑128 (MQ=255)
                                 aGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTa                                                                                                  <  6:191303/149‑1 (MQ=255)
                                     aGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTACTGACCGACatcag                                                                                                                             <  2:80147/118‑4 (MQ=255)
                                     aGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTACTGACCGACatcag                                                                                                                             >  1:80147/1‑115 (MQ=255)
                                      gCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTgg                                                                                             <  1:317306/149‑1 (MQ=255)
                                          aGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGtt                                                                                         >  7:29375/1‑149 (MQ=255)
                                          aGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGtt                                                                                         >  6:120854/1‑149 (MQ=255)
                                           gCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTg                                                                                        <  5:164816/149‑1 (MQ=255)
                                           gCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTg                                                                                        <  1:148575/149‑1 (MQ=255)
                                                   ttcGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTACTGACCGACatca                                                                                                                              >  6:218510/1‑101 (MQ=255)
                                                   ttcGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTACTGACCGACatca                                                                                                                              <  5:218510/103‑3 (MQ=255)
                                                     cGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCt                                                                                               <  8:91634/132‑1 (MQ=255)
                                                     cGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCt                                                                                               >  7:91634/1‑132 (MQ=255)
                                                     cGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCcttct                                                                               <  8:29375/148‑1 (MQ=255)
                                                     cGAAACTGAATTTGTGCAGCGGGCCGAGGTTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCcttc                                                                                <  4:59361/147‑1 (MQ=255)
                                                           tGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGa                                                                                                      <  1:186077/119‑1 (MQ=255)
                                                           tGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGa                                                                                                      >  2:186077/1‑119 (MQ=255)
                                                              aTTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTcgc                                                                     <  6:277301/149‑1 (MQ=255)
                                                                   tgCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATc                                                                <  1:128374/149‑1 (MQ=255)
                                                                         gggCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCtt                                                          <  5:120854/149‑1 (MQ=255)
                                                                             cGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTAGCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGcc                                                      >  6:83594/1‑149 (MQ=255)
                                                                               aTGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCtt                                                          >  2:74666/1‑143 (MQ=255)
                                                                               aTGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCtt                                                          <  1:74666/143‑1 (MQ=255)
                                                                                     cTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTcc                                               <  1:279759/148‑1 (MQ=255)
                                                                                                   tGCCGCTGTAACCCAGTAGATTCTGGTGTACTGATTTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTTAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGtt                                <  8:110590/149‑1 (MQ=255)
                                                                                                     ccGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGtt                                <  6:190481/147‑1 (MQ=255)
                                                                                                     ccGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGtt                                >  5:190481/1‑147 (MQ=255)
                                                                                                               ccAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGaa                    <  6:195836/149‑1 (MQ=255)
                                                                                                                  ggAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGcc                                                      <  7:2008/112‑1 (MQ=255)
                                                                                                                  ggAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGcc                                                      >  8:2008/1‑112 (MQ=255)
                                                                                                                   gagaTTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGaacaac                <  3:72905/149‑1 (MQ=255)
                                                                                                                            gggTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATc          <  4:282258/146‑1 (MQ=255)
                                                                                                                             ggTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGAt                             <  4:234944/126‑1 (MQ=255)
                                                                                                                                  tGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGaa  >  7:148181/1‑148 (MQ=255)
                                                                                                                                          gTCAGTACCGACATGATGAAGACCTGGACGCCTTgg                                                                                                          >  5:191927/1‑36 (MQ=255)
                                                                                                                                          gTCAGTACCGACATGATGAAGACCTGGACGCCTTgg                                                                                                          <  6:191927/36‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                        
GAGCATATTTCCCAGCTTCCGCCCAACAACGGAAGTCAGCACAGCGGTTTCTTCGAAACTGAATTTGTGCAGCGGGCCGATGGTTCTGTTTTGGCGACGTGCCGCTGGAACCCAGGAGATTCTGGGGTACTGATGTCGGTCAGTACCGACATGATGAAGACCTGGACGCCTTGGGTGAACTACCTGGCGTTGGCTCCTTCTTCGTCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAA  >  NC_002947/1774279‑1774556

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