Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,622,914 A→G 9.1% intergenic (+1507/+3057) PP_5654 → / ← PP_4091 phosphotriesterase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,622,9140AG 9.1% 112.1 / 2.3 46intergenic (+1507/+3057)PP_5654/PP_4091phosphotriesterase/transposase
Reads supporting (aligned to +/- strand):  ref base A (24/16);  new base G (2/2);  total (26/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 = 9.82e-01

AGCGAAACCAACAGCGCTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCACTGCT  >  NC_002947/4622767‑4623062
                                                                                                                                                    |                                                                                                                                                    
aGCGAAACCAACAGCGCTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGAt                                                                                                                                                     >  3:136596/1‑149 (MQ=255)
                cTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGt                                                                                                                                     <  5:262234/149‑1 (MQ=255)
                   gCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGAt                                                                                                                                  <  2:299659/149‑1 (MQ=255)
                             gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc                                                                                                                                       <  6:175273/134‑1 (MQ=255)
                             gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc                                                                                                                                       <  6:179587/134‑1 (MQ=255)
                             gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc                                                                                                                                       >  5:179587/1‑134 (MQ=255)
                             gTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGGTCAGCGATTTCGAAc                                                                                                                                       >  5:175273/1‑134 (MQ=255)
                                 cGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCAcgacga                                                                                                                    <  3:74282/149‑1 (MQ=255)
                                         cccGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCATATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTgc                                                                                                            <  1:51012/149‑1 (MQ=255)
                                                 gCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCACGCATTGCAAACGTGATCCGGTGCACGCCCCGGTCGTGCGCCCCCag                                                                                                    >  1:131864/1‑148 (MQ=14)
                                                 gCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCAACTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGCTGCGGTGCACGACGAGGTCCTGCGCAacc                                                                                                      >  8:71653/1‑147 (MQ=255)
                                                   aaTCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAgg                                                                                                  >  4:92825/1‑149 (MQ=255)
                                                        gCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAAc                                                                                             >  5:130253/1‑149 (MQ=255)
                                                         cAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACt                                                                                            >  7:7370/1‑149 (MQ=255)
                                                         cAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACt                                                                                            >  7:2717/1‑149 (MQ=255)
                                                          aGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTg                                                                                           <  4:136596/149‑1 (MQ=255)
                                                               cttgTCCTGGCCGAACAAAGTGATCGCAACCCACTGTTCTTCGGCAGGGGGCACCTGCAACATTTCAAACTTAAAGGGCGGGGTTTCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCAc                                                                                      <  2:131864/146‑1 (MQ=2)
                                                               cTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCAc                                                                                      >  6:126720/1‑149 (MQ=39)
                                                                     cTGCCCGACCAACGGGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAc                                                                                <  8:209762/149‑1 (MQ=21)
                                                                        cccGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAcgcg                                                                             >  4:269568/1‑149 (MQ=37)
                                                                          cGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACgcgc                                                                            >  3:169396/1‑148 (MQ=37)
                                                                           gACCAACGTGAACTCAACCACCTGTACTTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAc                                                                                >  1:20376/1‑144 (MQ=255)
                                                                           gACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAc                                                                                <  2:20376/143‑1 (MQ=37)
                                                                               aaCGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGa                                                                       >  3:98710/1‑148 (MQ=37)
                                                                                  gTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGccac                                                                    >  8:284852/1‑148 (MQ=25)
                                                                                        tCAACCACCTGTACTACGGGAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCATGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCgga                                                             >  7:4402/1‑148 (MQ=21)
                                                                                        tCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGt                                                             >  7:184837/1‑149 (MQ=34)
                                                                                          aaCCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCg                                                           >  4:90284/1‑149 (MQ=34)
                                                                                                 tGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGcc                                                    >  4:271767/1‑149 (MQ=37)
                                                                                                  gTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCg                                                   >  3:215433/1‑149 (MQ=37)
                                                                                                       tacGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCAGGCGCTACGACACCAGCTGTCGCCTGGCCGGaaa                                               <  3:185740/149‑1 (MQ=255)
                                                                                                         cGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGGTTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGaaa                                               <  5:119730/147‑1 (MQ=255)
                                                                                                         cGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGaaa                                               <  7:46643/147‑1 (MQ=39)
                                                                                                               cGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTAtt                                        >  5:328270/1‑148 (MQ=39)
                                                                                                                      cTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACcg                                <  5:126720/149‑1 (MQ=39)
                                                                                                                      cTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACcg                                >  4:145586/1‑149 (MQ=39)
                                                                                                                       tGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTAACTAcc                                 >  3:230864/1‑147 (MQ=255)
                                                                                                                         ccccAGAGCAACCTCAACGGCCGGGTGTTCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGa                             <  4:169396/147‑1 (MQ=255)
                                                                                                                            cAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGc                          <  6:328270/149‑1 (MQ=39)
                                                                                                                                       gAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAg               >  1:148359/1‑149 (MQ=39)
                                                                                                                                             ccGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGtgctgc         <  8:184837/149‑1 (MQ=39)
                                                                                                                                             ccGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGATCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGtgctgc         <  6:130253/149‑1 (MQ=255)
                                                                                                                                               gggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTAtt                                        >  7:137081/1‑116 (MQ=255)
                                                                                                                                               gggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTAtt                                        <  8:137081/116‑1 (MQ=255)
                                                                                                                                                ggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCAc      >  8:183210/1‑149 (MQ=40)
                                                                                                                                                ggTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCAc      <  3:145586/149‑1 (MQ=40)
                                                                                                                                                    aTCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTTCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCActgct  >  1:309472/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGAAACCAACAGCGCTGGCTTGCTGCAGTACCGCTATGACCCGCTGGGCAATCTGCAGACCCTGGTCCTGCCCGACCAACGTGAACTCAACCACCTGTACTACGGCAGCGGGCACCTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAACGTGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCCTGGCCGGAAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCACTGCT  >  NC_002947/4622767‑4623062

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: