Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,914,280 T→C 14.3% intergenic (+116/+37) spuI → / ← argA glutamylpolyamine synthetase/amino acid N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,914,2800TC14.3% 70.4 / 3.9 28intergenic (+116/+37)spuI/argAglutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (14/10);  new base C (2/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

ACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCTCGCGGAACCA  >  NC_002947/5914139‑5914427
                                                                                                                                              |                                                                                                                                                   
ccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagccttggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCgg                                                                                                                                             >  1:658406‑M1/75‑150 (MQ=255)
ccggcggcgcgggcttcgctgtacaactaccagcgcaattcgaagattttcgagaagcccctgtaagcctgggtCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCgg                                                                                                                                             >  1:793599‑M1/75‑150 (MQ=255)
 cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCCAAATCCTGTGACTTACCTGTGGGAGCggg                                                                                                                                            <  1:222436/150‑1 (MQ=255)
             gCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                                <  1:501028/150‑1 (MQ=255)
                        aaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGccc                                                                                                                     >  2:234866/1‑150 (MQ=39)
                                          tCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCCGGGTAATTTCGAAAGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCGtct                                                                                                   >  2:589350/1‑150 (MQ=255)
                                                 cggcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAAt                                                                                            >  2:673710/1‑150 (MQ=255)
                                                  ggcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATc                                                                                           >  2:376986/1‑150 (MQ=255)
                                                   gagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGa                                                                                                                                                 >  2:536050‑M1/24‑95 (MQ=255)
                                                   gagaagcccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGACTGTGACTTACCTGTGGGa                                                                                                                                                 <  1:536050‑M1/72‑1 (MQ=255)
                                                         cccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                                                                                                   <  2:122836‑M1/70‑1 (MQ=255)
                                                         cccctgtaagcctgggcCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                                                                                                                   >  1:122836‑M1/18‑87 (MQ=255)
                                                              tgAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                                                                                                            >  1:431436/1‑105 (MQ=255)
                                                              tgAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                                                                                                            <  2:431436/105‑1 (MQ=255)
                                                                      gtgtCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTgc                                                                                 >  2:565727/1‑140 (MQ=255)
                                                                      gtgtCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTgc                                                                                 <  1:565727/140‑1 (MQ=255)
                                                                                              CCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCgg                                               <  2:789177‑M1/150‑1 (MQ=255)
                                                                                                  TCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                                                                                                            >  2:537729‑M1/1‑69 (MQ=255)
                                                                                                  TCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                                                                                                            <  1:537729‑M1/69‑1 (MQ=255)
                                                                                                    CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCtt                                                                                         <  2:423380‑M1/102‑1 (MQ=255)
                                                                                                    CCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCtt                                                                                         >  1:423380‑M1/1‑102 (MQ=255)
                                                                                                            tAATTTCGAAGGCTGAAGCCTGTGACTTACCTGCGGGAGCGGGCGTGCCcg                                                                                                                                    <  1:552686/51‑1 (MQ=255)
                                                                                                            tAATTTCGAAGGCTGAAGCCTGTGACTTACCTGCGGGAGCGGGCGTGCCcg                                                                                                                                    >  2:552686/1‑51 (MQ=255)
                                                                                                                        cTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCAc                                   <  2:319726/136‑1 (MQ=255)
                                                                                                                        cTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCAc                                   >  1:319726/1‑136 (MQ=255)
                                                                                                                           aaGCCTGTGACTTACCCTGCGGGAGCGGGCGTGCCCGCGAAGAGggc                                                                                                                          >  2:569065/1‑47 (MQ=11)
                                                                                                                           aaGCCTGTGACTTACCCTGCGGGAGCGGGCGTGCCCGCGAAGAGggc                                                                                                                          <  1:569065/47‑1 (MQ=11)
                                                                                                                                      ttACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCTCGCgg       >  2:254745/1‑150 (MQ=255)
                                                                                                                                            tgtgGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCTCGCGGAACca  >  2:471495/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAATCTTCGAATTGCGCTGGTAGTTGTACAGCGAAGCCCGCGCCGCCGGCAGCCGCTCCACGCCACTCGGCGCAAAGCCGCGCTCGCGGAACCA  >  NC_002947/5914139‑5914427

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