Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,355,046 Δ1 bp 9.8% coding (246/681 nt) modB → molybdate ABC transporter permease ModB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,355,0460G. 9.8% 155.8 / 11.0 51coding (246/681 nt)modBmolybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base G (18/28);  new base . (2/3);  total (20/31)
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.58e-01

GTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAGC  >  NC_002947/4354902‑4355192
                                                                                                                                                   |                                                                                                                                                   
gtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGGACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAg                                                                                                                                                 >  4:150108/1‑149 (MQ=255)
   gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGcc                                                                                                                                             >  2:24762/1‑149 (MQ=255)
     tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTg                                                                                                                                           >  7:5435/1‑149 (MQ=255)
       gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTggg                                                                                                                                         <  6:175286/149‑1 (MQ=255)
       gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                          <  4:32813/148‑1 (MQ=255)
       gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTg                                                                                                                                           <  2:17721/147‑1 (MQ=255)
        cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGGGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggt                                                                                                                                        <  7:19703/149‑2 (MQ=255)
         gcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCa                                                                                                                                       <  3:90413/149‑1 (MQ=255)
         gcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCa                                                                                                                                       <  7:186746/149‑1 (MQ=255)
         gcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCa                                                                                                                                       <  1:318720/149‑1 (MQ=255)
              cgcgCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                  <  6:200770/149‑1 (MQ=255)
              cgcgCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                  <  8:5435/149‑1 (MQ=255)
                     ttGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGtt                                                                                                                             >  7:51304/1‑147 (MQ=255)
                             gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTcaca                                                                                                                   <  2:276892/149‑1 (MQ=255)
                                        ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTca                                                                                                                     <  8:178048/136‑1 (MQ=255)
                                        ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTca                                                                                                                     >  7:178048/1‑136 (MQ=255)
                                          tgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTg                                                                                                      <  4:131173/149‑1 (MQ=255)
                                                     cACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                           <  8:51304/149‑1 (MQ=255)
                                                      aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                            <  8:159288/131‑1 (MQ=255)
                                                      aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                            >  7:159288/1‑131 (MQ=255)
                                                       cTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTa                                                                                         <  3:183676/149‑1 (MQ=255)
                                                          ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTCCTCCGCCCTCACCCCCCAATTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                           <  1:115956/147‑1 (MQ=255)
                                                            gCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCAAGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCa                                                                                    >  4:337096/1‑149 (MQ=255)
                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTc                                                                            >  7:207260/1‑149 (MQ=255)
                                                                    tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTc                                                                            >  6:226940/1‑149 (MQ=255)
                                                                      ccACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgt                                                                          >  4:61135/1‑149 (MQ=255)
                                                                      ccACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgt                                                                          <  8:207260/149‑1 (MQ=255)
                                                                                   tCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTc                                                                                                                      <  4:200142/92‑1 (MQ=255)
                                                                                   tCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTc                                                                                                                      >  3:200142/1‑92 (MQ=255)
                                                                                          ttACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAAt                                                       <  1:24762/148‑1 (MQ=255)
                                                                                            aCCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGc                                                     <  5:226940/148‑1 (MQ=255)
                                                                                               tGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                       >  5:317490/1‑95 (MQ=255)
                                                                                               tGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                       <  6:317490/95‑1 (MQ=255)
                                                                                                   gATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgc                                              <  3:61135/148‑1 (MQ=255)
                                                                                                    aTCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggtg                                                                      <  8:271919/123‑1 (MQ=255)
                                                                                                    aTCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggtg                                                                      >  7:271919/1‑123 (MQ=255)
                                                                                                    aTCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcgc                                            <  3:150108/149‑1 (MQ=255)
                                                                                                     tCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGcgcg                                           <  4:25051/149‑1 (MQ=255)
                                                                                                                ggCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGTGTGGCTCGGCGGTGTACTCCATGCCGTTCGTGGGGCCCCCCCTGCAGAATGCCTTTGGCGCGATCGGc                                     >  5:263925/1‑147 (MQ=255)
                                                                                                                            tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGACCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgc                                              <  1:306802/126‑1 (MQ=255)
                                                                                                                            tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGACCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgc                                              >  2:306802/1‑126 (MQ=255)
                                                                                                                              gctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTGCAGCTTCACAGGCCTGGTGATTGGCTCGACGGGGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                      <  6:83004/148‑1 (MQ=255)
                                                                                                                                   gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                  >  4:294212/1‑35 (MQ=255)
                                                                                                                                   gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                  <  3:294212/35‑1 (MQ=255)
                                                                                                                                        gCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAc           >  3:224885/1‑149 (MQ=255)
                                                                                                                                        gCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAc           >  2:79167/1‑149 (MQ=255)
                                                                                                                                          tgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGc                                     <  5:215894/119‑1 (MQ=255)
                                                                                                                                          tgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGc                                     >  6:215894/3‑121 (MQ=255)
                                                                                                                                            cccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGt                    <  7:271372/136‑1 (MQ=255)
                                                                                                                                            cccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGt                    >  8:271372/1‑136 (MQ=255)
                                                                                                                                               aGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa    <  4:224885/149‑1 (MQ=255)
                                                                                                                                                   gcTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAc                                                                  <  6:61286/83‑1 (MQ=255)
                                                                                                                                                   gcTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAc                                                                  >  5:61286/1‑83 (MQ=255)
                                                                                                                                                   gcTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGagc  >  4:307029/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCT‑‑TGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAGC  >  NC_002947/4354902‑4355192

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