Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,256,199 G→C 8.5% intergenic (+259/‑445) asd → / → PP_1990 aspartate‑semialdehyde dehydrogenase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,256,1990GC 8.5% 115.6 / 2.5 45intergenic (+259/‑445)asd/PP_1990aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base C (2/2);  total (24/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 = 9.95e-01

ATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAA  >  NC_002947/2256057‑2256342
                                                                                                                                              |                                                                                                                                                    
attgCGAAATTTCTGCTCTGACAGAGCCGTTGCAGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGccc                                                                                                                                                <  6:120123/146‑1 (MQ=255)
   gCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt                                                                                                                                              >  5:122817/1‑148 (MQ=255)
      aaaTGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCacc                                                                                                                                          <  8:251576/149‑3 (MQ=255)
             tGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCAggc                                                                                                                                   >  3:217761/1‑148 (MQ=255)
               ctctGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg                                                                                                                                    <  6:233177/146‑1 (MQ=255)
                       gagCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg                                                                                                                         <  6:122817/149‑1 (MQ=255)
                          ccgttgccgtcgccgttgCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgggg                                                                                                                      <  2:215867/149‑1 (MQ=255)
                                ccgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCAGGCCCCATCCTGGGGcccagcg                                                                                                               <  4:217761/149‑6 (MQ=14)
                                 cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt                                                                                                                                                    >  8:146654/1‑112 (MQ=255)
                                 cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt                                                                                                                                                    <  7:146654/112‑1 (MQ=255)
                                      ccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGTCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGcgcg                                                                                                        <  2:900/148‑1 (MQ=37)
                                      ccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCTT‑‑‑Gctcccgg                                                                                                        >  3:38991/1‑142 (MQ=255)
                                        gttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGtc                                                                                                     <  4:92443/149‑1 (MQ=35)
                                            cttttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGtcctc                                                                                                  <  4:198143/148‑1 (MQ=34)
                                             ttttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTcc‑‑c                                                                                              >  7:6765/1‑149 (MQ=34)
                                                 gcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑Cgg                                                                                            >  8:26171/1‑147 (MQ=32)
                                                  cttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGa                                                                                          <  7:26171/148‑1 (MQ=32)
                                                   ttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAg                                                                                         >  5:305313/1‑148 (MQ=32)
                                                    tttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGc                                                                                        >  5:314571/1‑148 (MQ=32)
                                                       gcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGc                                                                                    >  5:97079/1‑149 (MQ=32)
                                                         ttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGccc                                                                                  >  3:217945/1‑149 (MQ=32)
                                                            tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCg                                                                                 >  5:138586/1‑147 (MQ=32)
                                                            tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGa                                                                               <  6:305313/149‑1 (MQ=32)
                                                                   gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGTATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg                                                                        <  4:145396/149‑1 (MQ=17)
                                                                   gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg                                                                        <  4:217945/149‑1 (MQ=32)
                                                                   gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg                                                                        >  2:293671/1‑149 (MQ=32)
                                                                               gCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCg                                                            >  1:267569/1‑149 (MQ=32)
                                                                                             gACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCg                                              >  1:268407/1‑149 (MQ=34)
                                                                                                   aggcggCCGTGCGCTGGGTTGGGGGCGGTTGGTGAGGGGCAGGGATCCCCGTAAAGCGGTGGGCCCCTGTTTGGGGCCT‑‑‑GCTGCGGGGTCCTCC‑‑CGGTAGCAAGCCGGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                        <  8:6765/149‑1 (MQ=2)
                                                                                                     gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATgg                                      >  3:160175/1‑149 (MQ=34)
                                                                                                      tggccttgcgcTGGGCTTGCGGTGGGAGGTGACGGGCAGGGTGGCCCGTCAAGCGCTGGGCCCCTGGTTGGGGCCG‑‑‑GCAGGGGGGTCCTCC‑‑CGGGAGAAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATggg                                     <  6:97079/142‑1 (MQ=11)
                                                                                                         ccGTGCGCTTGGTTTGGGGTGGGTGGTTCCGGGCTTGGATGCCCGTCTAGCGTGGGCCCCCTGTAGGGGGCCG‑‑‑GCAGCGGGGTCCACC‑‑CGGCAGCAAGCCGGGTGGGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGc                                  <  2:268407/149‑1 (MQ=2)
                                                                                                           gAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGcc                                 >  7:53447/1‑148 (MQ=34)
                                                                                                             gcgcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGAc                              >  6:289682/1‑149 (MQ=34)
                                                                                                                         gggagggtgGTGACGGGCCTGGATGCCCGTCAAGCGGTGGGCCCCTGGTTGGGGCCTCTCTCTCGGGGGTCCTACATCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGttttt                       <  1:293671/148‑1 (MQ=2)
                                                                                                                            agggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCtt               <  5:289682/149‑1 (MQ=255)
                                                                                                                             agcgcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttt              >  4:11131/7‑149 (MQ=255)
                                                                                                                                gcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGt                                                           <  6:3948/98‑1 (MQ=255)
                                                                                                                                gcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGt                                                           >  5:3948/4‑101 (MQ=255)
                                                                                                                                  gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGc          >  1:9825/1‑148 (MQ=255)
                                                                                                                                       gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt            <  7:110770/141‑1 (MQ=255)
                                                                                                                                       gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt            >  8:110770/1‑141 (MQ=255)
                                                                                                                                       gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGa     >  3:278061/1‑148 (MQ=255)
                                                                                                                                       gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGa     >  1:181579/1‑148 (MQ=255)
                                                                                                                                          cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACaa  >  3:88009/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAA  >  NC_002947/2256057‑2256342

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