Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,425,798 A→G 9.2% T331A (ACC→GCC)  U712_07320 → Spore photoproduct lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,425,7980AG 9.2% 116.1 / 2.3 45T331A (ACC→GCC) U712_07320Spore photoproduct lyase
Reads supporting (aligned to +/- strand):  ref base A (29/11);  new base G (2/2);  total (32/13)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGAAGAAAAGCGCCGTTATAAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAATTGC  >  CP006881/1425690‑1425936
                                                                                                            |                                                                                                                                          
gATGAAGAAAAGCGCCGTTATAAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAatat                                                                                                               <  8:50129/138‑1 (MQ=255)
      gaaAAGCGCCGTTATATATGGGGCCGTTTCGGGATCGGATATTATTTTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCAGTTGAATCCTATATTGATTCCTTTTTCCCTAAGGCAAATATTGAATATTTCACt                                                                                                         >  5:341900/1‑138 (MQ=255)
      gaaAAGCGCCGTTATAAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACt                                                                                                         >  5:1204/1‑138 (MQ=255)
        aaaTCGCCGTTAATAATGGGGCCGTTTCGGGATCGGTAAATAAATTTATCATAAAGATGCAGCCCATGGAATTCGAGCCGCACTTGTATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTaa                                                                                                      <  6:86827/139‑1 (MQ=255)
        aaaGCGCCGTTATAAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTaa                                                                                                      <  6:90364/139‑1 (MQ=255)
                 ttATAAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCtgt                                                                                             >  4:379875/1‑139 (MQ=255)
                   ataAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCtgtt                                                                                            <  6:1204/138‑1 (MQ=255)
                           ggCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACggg                                                                                                 >  8:245951/1‑125 (MQ=255)
                           ggCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACggg                                                                                                 <  7:245951/125‑1 (MQ=255)
                           ggCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCaa                                                                                    >  1:181924/1‑138 (MQ=255)
                           ggCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCaa                                                                                    >  5:8368/1‑138 (MQ=255)
                            gCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAACCGGGCTGTTGTGATCAAcc                                                                                  >  8:78100/1‑138 (MQ=255)
                            gCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCaa                                                                                    >  1:129547/1‑137 (MQ=255)
                            gCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCaa                                                                                    >  8:15404/1‑137 (MQ=255)
                                     ggATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCttt                                                                         >  6:99323/1‑139 (MQ=255)
                                                          aGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTtgtg                                                                                         >  3:93639/1‑102 (MQ=255)
                                                          aGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTtgtg                                                                                         <  4:93639/102‑1 (MQ=255)
                                                                      agCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAAc                                        >  1:261933/1‑139 (MQ=255)
                                                                        cATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAACCGc                                      >  7:147196/1‑139 (MQ=255)
                                                                        cATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGc                                      >  6:45276/1‑139 (MQ=255)
                                                                          tGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCg                                    >  1:114589/1‑139 (MQ=255)
                                                                          tGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCg                                    >  8:320895/1‑139 (MQ=255)
                                                                          tGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCg                                    >  8:219014/1‑139 (MQ=255)
                                                                          tGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCg                                    >  8:132429/1‑139 (MQ=255)
                                                                          tGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCg                                    >  3:86630/1‑139 (MQ=255)
                                                                          tGCACTTAGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTAAAACGGGATGTTGTGATCAACAGCTCGCTATTTTATAAAAAACATTTTTCAATTCCATCAAAACGCCg                                    >  2:246353/1‑139 (MQ=255)
                                                                                       gCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCaaaa                                         >  8:185645/1‑121 (MQ=255)
                                                                                       gCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCaaaa                                         <  7:185645/121‑1 (MQ=255)
                                                                                       gCACTTGAATCCTATATAGATACCTTTTTCCCTAACGCAAAAATTCAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTATTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGAtt                       >  1:370287/1‑139 (MQ=255)
                                                                                        cACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTT‑AACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCt                     <  2:224861/139‑1 (MQ=255)
                                                                                        cACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTT‑AACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCt                     <  2:219609/139‑1 (MQ=255)
                                                                                                  cTATATTGATGCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTAtcctc             >  8:181431/1‑138 (MQ=255)
                                                                                                  cTATATTGATGCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTAAATGATTCTTTAtcctc             >  8:175745/1‑138 (MQ=255)
                                                                                                   tataTTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTAtcctcc            >  5:93560/1‑138 (MQ=255)
                                                                                                     tatTGATGCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCa           <  7:181431/137‑1 (MQ=255)
                                                                                                     tatTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCaaa         >  7:38734/1‑139 (MQ=255)
                                                                                                     taatGATGCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACCTTTTTCCATTCCATCAAACCGCCGTAATTCAATGATTCTTTATCCTCCa           <  7:175745/134‑1 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAACCAGTTTTCCATTCCATCAAAACGCCGTCTTTCATTGATTCTTTATCCTCGAACTTAtt    >  2:265344/1‑139 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCaaaa                                         <  1:97744/102‑1 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCaaaa                                         >  2:97744/1‑102 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAAtt    >  2:211059/1‑139 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAAtt    >  5:91182/1‑139 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAAtt    >  5:8619/1‑139 (MQ=255)
                                                                                                          aTACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCAATCAAAACGCCTTAATTCAATGATTCTTTATCATCCAAATAAtt    >  7:176613/1‑139 (MQ=255)
                                                                                                           tACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAATTg   <  5:353553/139‑1 (MQ=255)
                                                                                                            aCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAATTGc  >  6:185926/1‑139 (MQ=255)
                                                                                                            aCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAATTGc  >  4:251603/1‑139 (MQ=255)
                                                                                                            aCCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCACCAAATAATTGc  >  8:3928/1‑139 (MQ=255)
                                                                                                            |                                                                                                                                          
GATGAAGAAAAGCGCCGTTATAAATGGGGCCGTTACGGGATCGGAAAATATATTTATCAGAAAGATGAAGAGCATGCACTTCGAGAGGCACTTGAATCCTATATTGATACCTTTTTCCCTAACGCAAAAATTGAATATTTCACTTAAACGGGCTGTTGTGATCAACAGCTCGCTTTTTTATAAAAAACATTTTTCCATTCCATCAAAACGCCGTAATTCAATGATTCTTTATCCTCCAAATAATTGC  >  CP006881/1425690‑1425936

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