Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 185,117 A→G 8.9% intergenic (‑48/+6) yaeI ← / ← dapD phosphodiesterase with model substrate bis‑pNPP/2,3,4,5‑tetrahydropyridine‑2‑carboxylate N‑succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913185,1170AG 8.9% 130.7 / 4.5 45intergenic (‑48/+6)yaeI/dapDphosphodiesterase with model substrate bis‑pNPP/2,3,4,5‑tetrahydropyridine‑2‑carboxylate N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (23/18);  new base G (2/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

TTCAAACCAGCCTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCG  >  NC_000913/184971‑185265
                                                                                                                                                  |                                                                                                                                                    
ttCAAACCAGCCTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATa                                                                                                                                                    >  4:439377/1‑149 (MQ=255)
ttCAAACCAGCCTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATa                                                                                                                                                    >  4:618057/1‑149 (MQ=255)
ttCAAACCAGCCTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATa                                                                                                                                                    >  1:590811/1‑149 (MQ=255)
      ccAGCCTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTa                                                                                                                                              <  4:144462/149‑1 (MQ=255)
           cTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGa                                                                                                                                         >  4:174301/1‑149 (MQ=255)
           cTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGa                                                                                                                                         <  3:149149/149‑1 (MQ=255)
                    aGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACgca                                                                                                                                >  1:118316/1‑149 (MQ=255)
                             tataACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGt                                                                                                                       <  3:618057/149‑1 (MQ=255)
                                    gAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATg                                                                                                                  <  4:246664/147‑1 (MQ=255)
                                    gAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATg                                                                                                                  >  3:246664/1‑147 (MQ=255)
                                         cTGGGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCCTCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTc                                                                                                                                           <  2:478995/117‑1 (MQ=255)
                                         cTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTc                                                                                                                                           >  1:478995/1‑117 (MQ=255)
                                          tGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACttt                                                                                                          >  2:157196/1‑149 (MQ=255)
                                                 gTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGa                                                                                                   >  4:639582/1‑149 (MQ=255)
                                                  tGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAg                                                                                                  >  1:87273/1‑149 (MQ=255)
                                                         cgtcgCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTcgc                                                                                           <  3:439377/149‑1 (MQ=255)
                                                          gtcgCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTcgcg                                                                                          >  1:306325/1‑149 (MQ=255)
                                                           tcgCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGt                                                                                         >  2:154733/1‑149 (MQ=255)
                                                                     aGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCtt                                                                               <  1:157196/149‑1 (MQ=255)
                                                                        cTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAAc                                                                            >  2:410149/1‑149 (MQ=255)
                                                                                       gcgTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTa                                                             <  2:333180/149‑1 (MQ=255)
                                                                                         gTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTaga                                                           >  1:501833/1‑149 (MQ=255)
                                                                                            aaaTCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCGTACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGc                                                        >  3:23843/1‑149 (MQ=255)
                                                                                            aaaTCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCGTACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGc                                                        >  4:282281/1‑149 (MQ=255)
                                                                                               tcatcGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGt                                                     <  2:87273/149‑1 (MQ=255)
                                                                                               tcatcGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGt                                                     <  4:252098/149‑1 (MQ=255)
                                                                                                 atcGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTAt                                                   >  3:523669/1‑149 (MQ=255)
                                                                                                   cGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCGTACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTAttt                                                 <  4:23843/149‑1 (MQ=255)
                                                                                                   cGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCGTACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTAttt                                                 <  3:282281/149‑1 (MQ=255)
                                                                                                    gCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTg                                                <  3:174301/149‑1 (MQ=255)
                                                                                                     cTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGc                                               <  2:675158/149‑1 (MQ=255)
                                                                                                               cAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGa                                     <  2:501833/149‑1 (MQ=255)
                                                                                                                aaaaaaaTCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTa                                                             <  1:123997/124‑1 (MQ=255)
                                                                                                                aaaaaaaTCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTa                                                             >  2:123997/1‑124 (MQ=255)
                                                                                                                aaaaaaaTCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGAc                                    >  2:401378/1‑149 (MQ=255)
                                                                                                                          ggCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTGGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGGAGATTAc                          >  4:655265/1‑149 (MQ=255)
                                                                                                                          ggCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTAc                          >  4:477344/1‑149 (MQ=255)
                                                                                                                           gCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTAcc                         <  2:590811/149‑1 (MQ=255)
                                                                                                                               cacgacgctgcCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGaa                     >  1:326270/7‑149 (MQ=255)
                                                                                                                                gtcgtGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGaaa                    <  4:523669/149‑1 (MQ=255)
                                                                                                                                gtcgtGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGaaa                    <  1:410149/149‑1 (MQ=255)
                                                                                                                                    tgctgcCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAaacaac                >  4:119/1‑149 (MQ=255)
                                                                                                                                     gctgcCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAAcc               <  3:639582/149‑1 (MQ=255)
                                                                                                                                         cccGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAg           >  1:598925/1‑149 (MQ=255)
                                                                                                                                                gCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCccgc    >  3:690920/1‑149 (MQ=255)
                                                                                                                                                  aTACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCccgccg  >  2:686562/1‑149 (MQ=255)
                                                                                                                                                  aTACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTCTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCccgccg  >  4:359232/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTCAAACCAGCCTGGCTCACAGTAATGCATATAACCGAAACCTGAGCTCGTGGCTATCGTCGCGGCAGTAGCCTGCAAAAATCGGCGGCGTGAAATCATCGCTCATCCCTGCAAAAAAATCGGGCAGCGTCGTGCTGCCCGTGTGCATACTTTTAGTCGATGGTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCG  >  NC_000913/184971‑185265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑