Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,046,859 C→A 9.1% T8N (ACC→AAC)  PP_3568 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,046,8590CA 9.1% 131.6 / 4.4 44T8N (ACC→AAC) PP_3568hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base A (2/2);  total (23/21)
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.49e-01

TTGTCCATGTTCAGCCCGTCTTCCAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGCGAGTGGATCATG  >  NC_002947/4046712‑4047007
                                                                                                                                                     |                                                                                                                                                    
ttGTCCATGTTCAGCCCGTCTTCCAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAcc                                                                                                                                                     >  7:1125/1‑149 (MQ=255)
      aTGTTCAGCCCGTCTTCCAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAccgccctg                                                                                                                                               <  7:54539/149‑1 (MQ=255)
            aGCCCGTCTTCCAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAccgccctgccgcc                                                                                                                                          <  6:223304/148‑1 (MQ=255)
                  tCTTCCAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGAGCGCAAAGGAGATTAGTCGTCCATAGTCGAATCAACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACcgccctgccgccctgc                                                                                                                                      >  1:103343/1‑148 (MQ=255)
                       cAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGc                                                                                                                              <  6:22893/149‑1 (MQ=255)
                                          aCAAACAGCCGGGCGGGTCAAGGAGCGCGCGCCAAGGCGATT‑G‑CGTACATAGTGGAATCCACCAATTTCCAAGGAGGGAAACCATCACCCAGAACAACCTGGGGACCGCCCTGCCGCCCAGCACCCCCGCCAGCGCCGAACGTGTcgac                                                                                                           >  7:50497/1‑149 (MQ=255)
                                                 gCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGc                                                                                                    <  4:44733/149‑1 (MQ=255)
                                                    gggcgggTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCa                                                                                                 <  7:200508/149‑1 (MQ=255)
                                                          gTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAgccgcc                                                                                           >  5:137399/1‑149 (MQ=255)
                                                               gCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAccgccctgccgccctg                                                                                                                                       <  8:40770/100‑1 (MQ=255)
                                                               gCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAccgccctgccgccctg                                                                                                                                       >  7:40770/1‑100 (MQ=255)
                                                                  gcgcgcgcCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATc                                                                                   >  5:283081/1‑149 (MQ=255)
                                                                          cAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCGCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGa                                                                           <  2:103343/149‑1 (MQ=255)
                                                                          cAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCATGGCATCCGCGTg                                                                             >  1:244189/1‑147 (MQ=255)
                                                                           aaGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAg                                                                          >  3:266669/1‑149 (MQ=255)
                                                                           aaGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAg                                                                          >  7:143940/1‑149 (MQ=255)
                                                                            aGGCGATT‑G‑CGTCCAAAGTGGAATCCACCCATTTCCCAGGAGGCTTACCATGACCCTGAACACCCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGc                                                                         <  8:249526/149‑1 (MQ=255)
                                                                                gATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTAt                                                                     >  1:199684/1‑149 (MQ=255)
                                                                                     g‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGt                                                                  <  8:1125/148‑1 (MQ=255)
                                                                                        gTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCTTGGCATCCGCGTGGAGCGTATCGtctc                                                               >  3:195279/1‑149 (MQ=255)
                                                                                        gTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGtctc                                                               >  6:71874/1‑149 (MQ=255)
                                                                                         tCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGt                                                                              >  1:98285/1‑133 (MQ=255)
                                                                                         tCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGt                                                                              <  2:98285/133‑1 (MQ=255)
                                                                                           cATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCa                                                             <  2:199684/148‑1 (MQ=255)
                                                                                                  gAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAg                                                     <  6:283081/149‑1 (MQ=255)
                                                                                                    aTCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGc                                                   >  4:276050/1‑149 (MQ=255)
                                                                                                    aTCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGc                                                   <  5:71874/149‑1 (MQ=255)
                                                                                                     tCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAACGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGcc                                                  <  2:299747/149‑1 (MQ=255)
                                                                                                     tCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAACGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGcc                                                  <  2:295588/149‑1 (MQ=255)
                                                                                                     tCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAACGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGc                                                   >  1:299747/1‑148 (MQ=255)
                                                                                                     tCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGAACGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGc                                                   >  1:295588/1‑148 (MQ=255)
                                                                                                        aCCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGc                                               >  6:299275/1‑149 (MQ=255)
                                                                                                           cATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGcccc                                            >  7:264549/1‑148 (MQ=255)
                                                                                                            aTTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCAc                                           >  2:62673/1‑149 (MQ=255)
                                                                                                                 cAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCTTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGc                                      <  4:195279/149‑1 (MQ=255)
                                                                                                                 cAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGc                                      <  8:50497/149‑1 (MQ=255)
                                                                                                                   aggaggCTTACGATGACCCTGAACAACCTGCTGACCGCCCTGCCGCGCTGCGACCCCGCCAGCTCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGa                                                                           <  3:294365/110‑1 (MQ=255)
                                                                                                                   aggaggCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCCAACTGCGACCCCGCCAGCGCCCAACGTGTCGACGAACAGCTCAGAAGCCCTGGCATCCGCGTgaa                                                                           >  4:294365/1‑108 (MQ=255)
                                                                                                                       ggCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTgg                                <  1:62673/149‑1 (MQ=255)
                                                                                                                        gCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGt                               <  6:137399/149‑1 (MQ=255)
                                                                                                                          ttACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTAc                             >  8:305554/1‑149 (MQ=255)
                                                                                                                                   cccTGAACAACCTGCTGACCGCCCTGCCGACCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGTCTTCTGGTACGACCAGGcc                    <  2:244189/149‑1 (MQ=255)
                                                                                                                                    ccTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCcgc                                                                                >  2:133514/1‑88 (MQ=255)
                                                                                                                                    ccTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCcgc                                                                                <  1:133514/88‑1 (MQ=255)
                                                                                                                                    ccTGAACAACCTGCTGACCGCCCTGACGCCCTGCGACCCCGCCAGCGCAGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCATCTGGTACGACCAGGcc                    >  1:48205/1‑148 (MQ=255)
                                                                                                                                        aacaacCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGaa                 <  3:276050/147‑1 (MQ=255)
                                                                                                                                          caacCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGc              >  6:54293/1‑148 (MQ=255)
                                                                                                                                             ccTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGCGAGt          <  7:305554/149‑1 (MQ=255)
                                                                                                                                                   gACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGCGAGTGGATCa    >  3:95633/1‑149 (MQ=255)
                                                                                                                                                     ccgccctgccgccctgCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGCGAGTGGATCATg  <  5:54293/149‑1 (MQ=255)
                                                                                                                                                     ccgccctgccgccctgCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGCGAGTGGATCATg  <  8:143940/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGTCCATGTTCAGCCCGTCTTCCAAAACTCGGAAAGTTCTTACAAACAGCCGGGCGGGTCAAGCCGCGCGCGCCAAGGCGATT‑G‑CGTCCATAGTGGAATCCACCCATTTCCAAGGAGGCTTACCATGACCCTGAACAACCTGCTGACCGCCCTGCCGCCCTGCGACCCCGCCAGCGCCGAACGTGTCGACGAACTGCTCAGCCGCCCTGGCATCCGCGTGGAGCGTATCGTCTCCAGCGGGCAGGCCAGCCCACCCGGCTTCTGGTACGACCAGGCCGAAGGCGAGTGGATCATG  >  NC_002947/4046712‑4047007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑