Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,897,761 A→G 7.7% I50I (ATT→ATC)  PP_3440 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,897,7610AG 7.7% 145.6 / 2.4 52I50I (ATT→ATC) PP_3440hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (2/2);  total (27/25)
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.72e-01

GCAACGTGTGGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGATGTCGC  >  NC_002947/3897622‑3897903
                                                                                                                                           |                                                                                                                                              
gCAACGTGTGGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCaa                                                                                                                                       >  4:320059/1‑149 (MQ=255)
     gtgtGGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCtt                                                                                                                                  <  3:320059/149‑1 (MQ=255)
      tgtgGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTtg                                                                                                                                 <  7:112938/149‑1 (MQ=255)
      tgtgGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTtg                                                                                                                                 >  4:20446/1‑149 (MQ=255)
      tgtgGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTtg                                                                                                                                 >  8:200236/1‑149 (MQ=255)
      tgtgGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTtg                                                                                                                                 <  1:367069/149‑1 (MQ=255)
                tCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGc                                                                                                                       >  7:309010/1‑149 (MQ=255)
                  ccGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTacc                                                                                                                                         >  2:234813/1‑129 (MQ=255)
                   cGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGctct                                                                                                                    >  5:261744/1‑149 (MQ=255)
                      cATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTga                                                                                                                 <  3:351919/149‑1 (MQ=255)
                              gcgcCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCTGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTga                                                                                                                 <  2:267716/141‑1 (MQ=255)
                              gcgcCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTga                                                                                                                 >  1:267716/1‑141 (MQ=255)
                                 ccGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGcc                                                                                                      >  8:167643/1‑149 (MQ=255)
                                          cGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCa                                                                                             <  7:200236/149‑1 (MQ=255)
                                                  tAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCg                                                                                     <  2:157607/149‑1 (MQ=255)
                                                    aaaCCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGtt                                                                                   >  3:14115/1‑149 (MQ=255)
                                                    aaaCCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGtt                                                                                   >  5:160051/1‑149 (MQ=255)
                                                         gCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTgcgc                                                                              >  4:109291/1‑149 (MQ=255)
                                                              ccTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCa                                                                          <  4:92535/148‑1 (MQ=255)
                                                               cTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATg                                                                        >  3:227179/1‑149 (MQ=255)
                                                                 tCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCGATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTa                                                                      >  1:301569/1‑149 (MQ=255)
                                                                 tCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCGATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTa                                                                      >  7:243497/1‑149 (MQ=255)
                                                                   ccAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGcc                                                                                                      <  5:206283/115‑1 (MQ=255)
                                                                   ccAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGcc                                                                                                      >  6:206283/1‑115 (MQ=255)
                                                                   ccAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTAt                                                                     <  4:62737/148‑1 (MQ=255)
                                                                    cAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCg                                                                   <  7:167643/149‑1 (MQ=255)
                                                                       cccGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTa                                                                <  6:121224/149‑1 (MQ=255)
                                                                        ccGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAg                                                               <  4:227179/149‑1 (MQ=255)
                                                                        ccGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAg                                                               <  3:109291/149‑1 (MQ=255)
                                                                          gTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAggc                                                             <  1:352300/149‑1 (MQ=255)
                                                                             aCCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCa                                                                                         >  8:356970/1‑118 (MQ=255)
                                                                             aCCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCa                                                                                         <  7:356970/118‑1 (MQ=255)
                                                                             aCCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGGATCGGTAggcggc                                                          >  8:321691/1‑149 (MQ=255)
                                                                              cccGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGcc                                                         <  4:14115/149‑1 (MQ=255)
                                                                                           gCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGgtgt                                            <  3:20446/149‑1 (MQ=255)
                                                                                           gCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGgt                                              <  1:402912/147‑1 (MQ=255)
                                                                                             caggcaggAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGgtgtgt                                          >  5:322333/1‑149 (MQ=255)
                                                                                              aggcaggAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGtt                                         >  7:100458/1‑149 (MQ=255)
                                                                                                      aaaTGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCt                                 <  8:100458/149‑1 (MQ=255)
                                                                                                            ccGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGc                           >  7:337910/1‑149 (MQ=255)
                                                                                                                 ggTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATg                       >  3:199136/1‑148 (MQ=255)
                                                                                                                    gAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTCGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACt                   >  2:83294/1‑149 (MQ=255)
                                                                                                                       aCCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGc                               >  6:82797/1‑134 (MQ=255)
                                                                                                                       aCCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGc                               <  5:82797/134‑1 (MQ=255)
                                                                                                                               ttCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGc            >  1:64452/1‑145 (MQ=255)
                                                                                                                               ttCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGc            >  2:346661/1‑145 (MQ=255)
                                                                                                                                 cGCCCGGGGCGATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGAt       <  2:301569/148‑1 (MQ=255)
                                                                                                                                 cGCCCGGGGCGATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGAt       <  8:243497/148‑1 (MQ=255)
                                                                                                                                  gCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCCATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGAATCGGTAGGCGGCCAAGGTCAAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGATGt     >  7:348402/1‑149 (MQ=255)
                                                                                                                                   cccGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGATGTc    >  6:281730/1‑149 (MQ=255)
                                                                                                                                    ccGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGATGTCg   <  6:261744/149‑1 (MQ=255)
                                                                                                                                     cGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGATGTCGc  <  8:337910/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GCAACGTGTGGATACCTCCCGCCATGGTCTGCGCCGGGTGCACGTAGGGGTAAAACCGCTTGCCTTCCCAGCCCGTCACCCGATCCAGCTGGCCAGGCAGGAAAATGGCCGGCGGTGAAACCGTGGCTTCGCCCGGGGCAATGTCCCAACTCTTGGTGGCGACGCTCTTGAGATCGATCGCCGTTTTTCCACCCAGCCGTTTGCGCGCCATGTATCGGTAGGCGGCCAAGGTCCAGGTGTGTTGCCTGGCTGCGGGCGATGCACTGAACTGCGAGATGTCGC  >  NC_002947/3897622‑3897903

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