Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,256,198 G→C 8.9% intergenic (+258/‑446) 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,1980GC 8.9% 116.6 / 4.7 45intergenic (+258/‑446)asd/PP_1990aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base G (22/19);  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 = 8.97e-01

CTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAA  >  NC_002947/2256050‑2256342
                                                                                                                                                    |                                                                                                                                                     
cTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATg                                                                                                                                                       <  6:887/149‑1 (MQ=255)
cTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTCCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATg                                                                                                                                                       <  7:193004/149‑1 (MQ=255)
       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                                                                                                                                                    <  7:146654/112‑1 (MQ=255)
                                        cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt                                                                                                                                                    >  8:146654/1‑112 (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                                                                        >  2:293671/1‑149 (MQ=32)
                                                                          gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg                                                                        <  4:217945/149‑1 (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                                                           >  5:3948/4‑101 (MQ=255)
                                                                                                                                       gcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGt                                                           <  6:3948/98‑1 (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)
                                                                                                                                                    |                                                                                                                                                     
CTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAA  >  NC_002947/2256050‑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: