Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,891,173 G→C 9.5% K171N (AAG→AAC)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,891,1730GC 9.5% 121.9 / 4.8 42K171N (AAG→AAC) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base C (2/2);  total (21/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 = 6.06e-01

CCTGACCACGCCTTGCTGTATCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGTGGAATGGTTCAAGG  >  NC_002947/2891030‑2891321
                                                                                                                                               |                                                                                                                                                    
ccTGACCACGCCTTGCTGTATCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCg                                                                                                                                                  <  7:105674/148‑1 (MQ=255)
   gACCACGCCTTGCTGTATCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGa                                                                                                                                              >  1:95234/1‑149 (MQ=255)
    aCCACGCCTTGCTGTATCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGa                                                                                                                                              <  4:265161/148‑1 (MQ=255)
            ttGCTGTATCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCt                                                                                                                                     >  5:51730/1‑149 (MQ=255)
                    tCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCAc                                                                                                                             >  2:366347/1‑149 (MQ=255)
                             ggTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGCGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAAcg                                                                                                                    >  5:33887/1‑149 (MQ=255)
                              gTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAAcgc                                                                                                                   >  6:170377/1‑149 (MQ=255)
                               tGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAAcgcg                                                                                                                  <  6:51730/149‑1 (MQ=255)
                                gcgcCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGCCTTGCAGCTGCCGCCACGGCACAACgcgg                                                                                                                 >  4:131543/1‑149 (MQ=255)
                                  gcCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAAcgc                                                                                                                   >  4:19271/1‑145 (MQ=255)
                                  gcCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAAcgc                                                                                                                   <  3:19271/145‑1 (MQ=255)
                                          aGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTctacta                                                                                                       >  1:243017/1‑149 (MQ=255)
                                           gCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCTAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCtactac                                                                                                      <  3:117428/149‑1 (MQ=255)
                                              aaGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGAt                                                                                                   <  1:183852/149‑1 (MQ=255)
                                              aaGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGAt                                                                                                   >  7:351239/1‑149 (MQ=255)
                                                      aGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCCGCGTTCTACTACGATTTGCTCAg                                                                                           >  4:381133/1‑149 (MQ=255)
                                                        ccTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAACTTCGCCgc                                                                                                                                              >  1:292086/1‑95 (MQ=255)
                                                        ccTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAACTTCGCCgc                                                                                                                                              <  2:292086/96‑2 (MQ=255)
                                                                tGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTgc                                                                                 <  8:351239/149‑1 (MQ=255)
                                                                         cGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTg                                                                        <  8:116217/149‑1 (MQ=255)
                                                                          ggACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCTCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACTATTTGCTCAGCCCGCTGTGCGCACGGGTg                                                                        >  8:81536/1‑148 (MQ=255)
                                                                                ccTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACgtgt                                                                 >  6:372459/1‑149 (MQ=255)
                                                                                 cTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGtgtg                                                                >  5:264948/1‑149 (MQ=255)
                                                                                    gACGAACCGGCCCACCGGCAGTGGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACCACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGgc                                                              <  1:77824/148‑1 (MQ=255)
                                                                                          ccGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCAc                                                                             >  5:68036/1‑127 (MQ=255)
                                                                                          ccGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCAc                                                                             <  6:68036/127‑1 (MQ=255)
                                                                                               ccACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTatc                                                  <  2:143623/149‑1 (MQ=255)
                                                                                               ccACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTat                                                   >  1:93343/1‑148 (MQ=255)
                                                                                               ccACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTat                                                   >  3:191466/1‑148 (MQ=255)
                                                                                                       aGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTctac                                                                                                         <  5:315560/86‑1 (MQ=255)
                                                                                                       aGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTctac                                                                                                         >  6:315560/1‑86 (MQ=255)
                                                                                                           gcgcGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTtg                                       <  5:170377/148‑1 (MQ=255)
                                                                                                            cgcgAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGTGCACCACCTATCATCACCCGTtggt                                     <  2:243017/149‑1 (MQ=255)
                                                                                                             gcgAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTtggtg                                    <  2:93343/149‑1 (MQ=255)
                                                                                                               gagaTCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTggg                                  >  5:80425/1‑149 (MQ=255)
                                                                                                                          ccAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAg                       <  6:33887/149‑1 (MQ=255)
                                                                                                                            aaGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCg                     >  1:66835/1‑149 (MQ=255)
                                                                                                                                  ccTGGGCGGCGAAGTTCTCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGTg               <  7:81536/149‑1 (MQ=255)
                                                                                                                                   caagtcggcgAACTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGt                >  1:383207/6‑147 (MQ=255)
                                                                                                                                   caagtcggcgAACTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGt                <  2:383207/142‑1 (MQ=255)
                                                                                                                                   cTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGTgg              <  6:264948/149‑1 (MQ=255)
                                                                                                                                      ggcggcGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGTGGAAt           <  4:191466/149‑1 (MQ=255)
                                                                                                                                               gTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGTGGAATGGTTCAAgg  <  3:381133/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCTGACCACGCCTTGCTGTATCCGCACCTGGTGCGCCAGCTCAGCGAAGGCGCCAGCCTGGCCGTGCAGACCCCGGACAACCTTGACGAACCGGCCCACCGGCAGTTGCGCGAGATCGCAAGCCAAGGCCCCTGGGCGGCGAAGTTCGCCGACTTGCAGCTGCCGCCACGGCACAACGCGGCGTTCTACTACGATTTGCTCAGCCCGCTGTGCGCACGGGTGGACGTGTGGCGCACCACCTATCATCACCCGTTGGTGGGGGGCGCCGAAGCGGTAGTGGAATGGTTCAAGG  >  NC_002947/2891030‑2891321

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