Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,246:1 +GA 100% intergenic (+46/+96) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2461.G86.7% 127.5 / 15.6 45intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base G (18/21);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 6.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029471,070,2462.A86.7% 119.8 / 15.8 45intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base A (18/21);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 6.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAC  >  NC_002947/1070100‑1070378
                                                                                                                                                      ||                                                                                                                                    
tGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                        >  6:118282/1‑147 (MQ=255)
  aCGGCGCCGATTCGCATCTCTGTTATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Ga                                                                                                                                    <  7:299402/147‑1 (MQ=255)
  aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                        >  5:361793/1‑145 (MQ=255)
  aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                        <  6:361793/145‑1 (MQ=255)
  aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Ga                                                                                                                                    <  4:154174/147‑1 (MQ=255)
    gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCt‑‑gt                                                                                                                                        >  1:240014/1‑143 (MQ=255)
      cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcga                                                                                                                              <  5:235533/149‑7 (MQ=255)
      cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcga                                                                                                                              <  5:156864/149‑7 (MQ=255)
       gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcgag                                                                                                                             <  2:201679/149‑8 (MQ=255)
            ttCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGACAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGcgccg                                                                                                                          >  3:342999/1‑149 (MQ=255)
               gCATCTCGGTATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCcgcgcg                                                                                                                      >  2:120837/1‑149 (MQ=255)
                cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCcgcgcg                                                                                                                      <  7:6300/149‑1 (MQ=255)
                   ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTGTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                    >  1:85633/1‑149 (MQ=255)
                   ctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTATCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcgggc                                                                                                                   >  5:269779/1‑149 (MQ=255)
                      ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                    <  7:52790/145‑1 (MQ=255)
                         aaTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGACGCGCGGGCGGCGCt                                                                                                             >  4:341747/1‑149 (MQ=255)
                          aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                             <  3:289936/148‑1 (MQ=37)
                           tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                           <  6:35883/149‑1 (MQ=38)
                             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGTGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                         <  6:225152/149‑1 (MQ=38)
                             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                          <  2:240014/148‑1 (MQ=25)
                             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                         >  5:225152/1‑149 (MQ=38)
                                aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                       <  5:118282/148‑1 (MQ=255)
                                     ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                 <  3:341747/149‑1 (MQ=255)
                                         tGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGGGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                             <  8:359441/149‑1 (MQ=255)
                                          ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                             <  4:342999/148‑1 (MQ=255)
                                               ggtgtgctgtgGGTGGCTTTCATGCTTAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                       <  2:382338/147‑1 (MQ=38)
                                                   tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTGGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                     <  1:120837/147‑1 (MQ=255)
                                                       gtgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcggg                                                                                                                    >  8:145075/1‑112 (MQ=255)
                                                          ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                                            <  1:241300/149‑1 (MQ=255)
                                                              gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGACGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                        <  4:246354/149‑1 (MQ=38)
                                                                   tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                    >  4:212940/1‑149 (MQ=255)
                                                                   tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                    >  4:140349/1‑149 (MQ=255)
                                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                              >  5:349853/1‑149 (MQ=255)
                                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                             <  2:262359/102‑1 (MQ=18)
                                                                         tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                             >  1:262359/1‑102 (MQ=18)
                                                                            agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                           <  3:140349/149‑1 (MQ=255)
                                                                                 gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                      >  5:171662/1‑149 (MQ=255)
                                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                            <  1:173075/93‑1 (MQ=25)
                                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                            >  2:173075/1‑93 (MQ=25)
                                                                                       cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                >  6:85975/1‑149 (MQ=255)
                                                                                           ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                       <  4:343187/106‑1 (MQ=255)
                                                                                           ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                       >  3:343187/1‑106 (MQ=255)
                                                                                            cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                           <  2:85633/149‑1 (MQ=255)
                                                                                                                    aaGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                   >  1:56303/1‑149 (MQ=255)
                                                                                                                     aGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                   <  6:269779/148‑1 (MQ=255)
                                                                                                                                gTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg       >  6:325707/1‑149 (MQ=255)
                                                                                                                                    ggTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa   <  2:56303/149‑1 (MQ=255)
                                                                                                                                           ggCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc  >  7:43961/1‑143 (MQ=255)
                                                                                                                                           ggCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc  <  8:43961/143‑1 (MQ=255)
                                                                                                                                                      ||                                                                                                                                    
TGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAC  >  NC_002947/1070100‑1070378

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: