Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,485,504 A→G 10.3% K32K (AAA→AAG)  spuC‑I → polyamine:pyruvate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,485,5040AG10.3% 100.7 / 3.7 40K32K (AAA→AAG) spuC‑Ipolyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base G (2/2);  total (21/19)
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.77e-01

CTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGCCGAGGTCGCCAGCCAGCAGATGAAG  >  NC_002947/2485357‑2485651
                                                                                                                                                     |                                                                                                                                                   
cTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATTGAGTTAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAg                                                                                                                                                    <  5:123156/149‑1 (MQ=255)
  gTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGc                                                                                                                                                  >  5:349551/1‑149 (MQ=255)
     aaaaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCa                                                                                                                                               >  1:58250/1‑149 (MQ=255)
      aaaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCAc                                                                                                                                              <  1:120419/149‑1 (MQ=255)
               tCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCAtcatca                                                                                                                                     >  5:98902/1‑149 (MQ=255)
                 gCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCAcc                                                                                                                                   >  1:65818/1‑149 (MQ=255)
                         gCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCaa                                                                                                                           <  6:98902/149‑1 (MQ=255)
                          cTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGTGAAAGGCCCACGCATCATCACCTCGGCCAAg                                                                                                                          <  1:212120/149‑1 (MQ=255)
                           ggcgcAACTTGATGGAGTGACGGCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAgg                                                                                                                         <  8:289055/148‑1 (MQ=255)
                             cgcAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTTGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAgg                                                                                                                         <  6:349551/147‑1 (MQ=255)
                              gcAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAGGGCCCACGcatc                                                                                                                                         >  8:37124/1‑130 (MQ=255)
                              gcAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAGGGCCCACGcatc                                                                                                                                         >  8:31881/1‑130 (MQ=255)
                              gcAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAGGGCCCACGcatc                                                                                                                                         <  7:37124/130‑1 (MQ=255)
                              gcAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAGGGCCCACGcatc                                                                                                                                         <  7:31881/130‑1 (MQ=255)
                               cAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCAGCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGgtgt                                                                                                                     <  2:177698/149‑1 (MQ=255)
                                              aaGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAg                                                                                                                          >  8:27182/1‑129 (MQ=255)
                                              aaGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAg                                                                                                                          <  7:27182/129‑1 (MQ=255)
                                                   aaTGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAgg                                                                                                                         >  8:142757/1‑125 (MQ=255)
                                                   aaTGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAgg                                                                                                                         <  7:142757/125‑1 (MQ=255)
                                                                    cGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGa                                                                                 <  2:65818/148‑1 (MQ=255)
                                                                     gCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGatg                                                                               >  8:22890/1‑149 (MQ=255)
                                                                           ccTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCCGTTCTGCGATGTCAGGCAACT‑GCCCGAGAACGGCCCACGCATCTTCACCTCGGCCAAGGGTGTGTACCTGGGGGCAAGCGAGGGCAACAAGATACTCGatggcatg                                                                          >  1:56741/1‑149 (MQ=255)
                                                                           ccTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGatggcatgg                                                                         >  2:235459/1‑149 (MQ=255)
                                                                                       aaCCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTgtggt                                                             <  2:56741/149‑1 (MQ=255)
                                                                                          cATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCg                                                          >  3:277526/1‑149 (MQ=255)
                                                                                           aTGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTgg                                                                                                           <  5:200652/99‑1 (MQ=255)
                                                                                           aTGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTgg                                                                                                           >  6:200652/1‑99 (MQ=255)
                                                                                              aGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATGATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTgg                                                       <  2:58250/148‑1 (MQ=255)
                                                                                                       cATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTa                                             >  7:99329/1‑149 (MQ=255)
                                                                                                        aTCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTAt                                            <  4:277526/149‑1 (MQ=255)
                                                                                                         tCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATg                                           >  8:67762/1‑149 (MQ=255)
                                                                                                         tCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGAAGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATg                                           <  1:235459/149‑1 (MQ=255)
                                                                                                            ccTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTc                                        >  2:169757/1‑149 (MQ=255)
                                                                                                             cTGG‑CCCCGTTCAGCGATGTCAGGCAACTGGGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTc                                        >  1:26391/1‑149 (MQ=255)
                                                                                                             cTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTcg                                       >  4:287003/1‑149 (MQ=255)
                                                                                                                     cGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATTTTGTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAagag                                <  1:114971/149‑1 (MQ=255)
                                                                                                                          aGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGc                           >  7:102378/1‑149 (MQ=255)
                                                                                                                                gTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGCCGAGGt                     >  4:295998/1‑149 (MQ=255)
                                                                                                                                   aGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGCCGAGGTCgc                  >  2:10236/1‑149 (MQ=255)
                                                                                                                                                 agaAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGCCGAGGTCGCCAGCCAGCAGATg     <  8:102378/149‑1 (MQ=255)
                                                                                                                                                    aaGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGCCGAGGTCGCCAGCCAGCAGATGAag  <  7:67762/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGG‑CCCCGTTCAGCGATGTCAGGCAACT‑GGCCGAGAAAGGCCCACGCATCATCACCTCGGCCAAGGGTGTGTACCTGTGGGACAGCGAGGGCAACAAGATTCTCGATGGCATGGCCGGCCTGTGGTGCGTGGCGGTCGGCTATGGTCGCGAAGAGCTGGCCGAGGTCGCCAGCCAGCAGATGAAG  >  NC_002947/2485357‑2485651

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