Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,160,633 A→G 9.7% L141P (CTG→CCG)  PP_4543 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,160,6330AG 9.7% 135.3 / 5.6 52L141P (CTG→CCG) PP_4543hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base G (2/3);  total (23/29)
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

TTCCATCCACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCG  >  NC_002947/5160488‑5160780
                                                                                                                                                 |                                                                                                                                                   
ttCCATCCACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGcc                                                                                                                                                  >  1:347412/1‑149 (MQ=255)
   catccaCCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTg                                                                                                                                               >  5:132908/1‑149 (MQ=255)
    atccaCCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCGAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTgg                                                                                                                                              >  5:6361/1‑149 (MQ=255)
     tccaCCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTggg                                                                                                                                             <  1:245333/149‑1 (MQ=255)
     tccaCCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTggg                                                                                                                                             <  4:318622/149‑1 (MQ=255)
        aCCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCg                                                                                                                                          >  5:183068/1‑149 (MQ=255)
          ccAGGCAGGACCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTGGCCCAGATACTCCGGCAGGGCCTGCCGCATCTACAAACCGGCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCt                                                                                                                                        >  5:267655/1‑149 (MQ=255)
                aggcCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCgg                                                                                                                                  <  1:352214/149‑1 (MQ=255)
                     cGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTg                                                                                                                              >  4:222195/1‑148 (MQ=255)
                       ccAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAAc                                                                                                                           <  7:339782/149‑1 (MQ=255)
                         aaCTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCTGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   <  8:250848/139‑3 (MQ=255)
                         aaCTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   >  7:250848/1‑137 (MQ=255)
                         aaCTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACgg                                                                                                                         <  4:170616/149‑1 (MQ=255)
                           cTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGca                                                                                                                       <  5:294497/149‑1 (MQ=255)
                              cgcgTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGcaca                                                                                                                     <  4:226002/148‑1 (MQ=255)
                                   ccAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGcc                                                                                                                >  3:310791/1‑148 (MQ=255)
                                     aGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAg                                                                                                             >  1:318450/1‑149 (MQ=255)
                                          cGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACgg                                                                                                                         >  1:136488/1‑132 (MQ=255)
                                          cGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACgg                                                                                                                         <  2:136488/132‑1 (MQ=255)
                                            ggTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGtttt                                                                                                      <  4:310791/149‑1 (MQ=255)
                                              tCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    <  6:6361/149‑1 (MQ=255)
                                              tCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    <  2:347412/149‑1 (MQ=255)
                                               cGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTggg                                                                                                   <  5:89247/149‑1 (MQ=255)
                                                 gTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   <  6:103390/115‑3 (MQ=255)
                                                 gTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   >  5:103390/1‑113 (MQ=255)
                                                    aaGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGgtagt                                                                                               <  6:132908/148‑1 (MQ=255)
                                                          aTGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGcac                                                                                                                      >  6:313949/1‑119 (MQ=255)
                                                           tgtgCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATa                                                                                       <  1:226824/149‑1 (MQ=255)
                                                              gcgcTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCg                                                                                    <  5:90584/149‑1 (MQ=255)
                                                                gcgccGGAAAGCCCACCCTTGCCCAGATCCTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCCGCCGTGGGCCGCTGACCGGGCTGTACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGc                                                                                  <  6:267655/146‑1 (MQ=255)
                                                                      gAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAAGCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTc                                                                            <  2:8255/149‑1 (MQ=255)
                                                                      gAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCt                                                                             >  2:169238/1‑148 (MQ=255)
                                                                      gAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCt                                                                             >  1:49953/1‑148 (MQ=255)
                                                                          gCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        <  1:169238/149‑1 (MQ=255)
                                                                          gCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        <  3:222195/149‑1 (MQ=255)
                                                                            ccACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        >  4:51947/1‑147 (MQ=255)
                                                                            ccACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        <  3:51947/147‑1 (MQ=255)
                                                                                      ccAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCaa                                                            >  7:337886/1‑149 (MQ=255)
                                                                                                      ggCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACAt                                            <  6:183068/149‑1 (MQ=255)
                                                                                                        ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCg                                                                  <  3:226842/125‑1 (MQ=255)
                                                                                                        ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCg                                                                  >  4:226842/1‑125 (MQ=255)
                                                                                                        ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGaa                                                 <  6:153090/142‑1 (MQ=255)
                                                                                                        ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGaa                                                 >  5:153090/1‑142 (MQ=255)
                                                                                                          tGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    >  7:251809/1‑89 (MQ=255)
                                                                                                          tGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    <  8:251809/89‑1 (MQ=255)
                                                                                                            cagcaTCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGc                                      <  2:127404/149‑1 (MQ=255)
                                                                                                                 tCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAAc                                 <  2:318450/149‑1 (MQ=255)
                                                                                                                            tCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGGCACtt                      >  1:80316/1‑149 (MQ=255)
                                                                                                                                        gCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCggg          >  7:287395/1‑149 (MQ=255)
                                                                                                                                         cggTCAGCGGCCCACGGCTGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGc                                                                <  4:366432/92‑1 (MQ=38)
                                                                                                                                         ccGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCggg          >  8:34504/1‑148 (MQ=255)
                                                                                                                                            tCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCCGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCaa       >  3:239717/1‑148 (MQ=255)
                                                                                                                                              agcagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACa     <  2:49953/148‑1 (MQ=255)
                                                                                                                                                cagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAgcg  >  4:173312/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTCCATCCACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCG  >  NC_002947/5160488‑5160780

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