Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,160,638 T→A 9.4% P139P (CCA→CCT) ‡ PP_4543 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,160,6380TA 9.4% 145.1 / 5.3 53P139P (CCA→CCT) ‡PP_4543hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/27);  new base A (2/3);  total (23/30)
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.65e-01

CATCCACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCG  >  NC_002947/5160491‑5160780
                                                                                                                                                   |                                                                                                                                              
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)
                                                                                                                                                   |                                                                                                                                              
CATCCACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCG  >  NC_002947/5160491‑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: ‑