Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,806 G→A 14.9% intergenic (+72/+70) PP_2589 → / ← PP_2590 aldehyde dehydrogenase family protein/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,958,8060GA14.9% 107.4 / 9.2 48intergenic (+72/+70)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (24/16);  new base A (3/4);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAG  >  NC_002947/2958672‑2958954
                                                                                                                                      |                                                                                                                                                    
tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                                           <  8:72241/146‑1 (MQ=255)
tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCt                                                                                                                                        >  5:262742/1‑149 (MQ=255)
tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                          <  4:159999/147‑1 (MQ=38)
tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTga                                                                                                                                            <  1:187001/145‑2 (MQ=38)
    tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCAAGCCCGCTCCCACAGGATCAGCACAACGCTGAGCGTACTGCGGGACCGGTgg                                                                                                                                    >  8:132662/1‑149 (MQ=255)
            ttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgc                                                                                                                            <  1:352254/149‑1 (MQ=255)
               gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgcg                                                                                                                         >  7:63175/1‑149 (MQ=255)
                aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                    >  8:214004/1‑137 (MQ=255)
                aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                    <  7:214004/137‑1 (MQ=255)
                       ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcgcg                                                                                                                   >  3:273588/1‑147 (MQ=255)
                         cGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                                 >  6:232388/1‑147 (MQ=255)
                            aaCTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGc                                                                                                            >  4:62293/1‑149 (MQ=255)
                                cAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTa                                                                                                        >  3:166298/1‑149 (MQ=255)
                                    accaccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGATCCAGCACAACCCTCAGACTTGTGCTga                                                                                                                                            <  3:149345/109‑2 (MQ=18)
                                    accaccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGATCCAGCACAACCCTCAGACTTGTGCTga                                                                                                                                            >  4:149345/1‑108 (MQ=18)
                                         cTGGTCCCAGTTGCGCGGGTGATGTATTGATCTTTGCGCGCGTCGTCTGCTTCGGGGGCATGCCCGCTCCCAGAGGTTCAGCACAAGTCTGAGGGTTGTGGTGGACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGcgcc                                                                                               <  4:274115/149‑1 (MQ=38)
                                              aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                          >  3:144129/1‑149 (MQ=255)
                                                      cgcgGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGGCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGAccc                                                                                   <  5:168403/148‑1 (MQ=255)
                                                          gcttttttATTGATCATTTACTGCGTTGTGTGCTTCGCGGGATTGCCGGCTCCCACTGGTTCAGCACATGTATGAGGGTTGTGCTGGACCTGTGGTTGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                              <  6:262742/142‑1 (MQ=255)
                                                              aTGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                                           <  2:40691/148‑1 (MQ=255)
                                                                    tGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                                     <  4:273588/148‑1 (MQ=255)
                                                                     gATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCAAAGGTTCAGCACAATTGTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGc                                                                   >  6:37206/1‑149 (MQ=255)
                                                                        cATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGGTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATATGTCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCAAg                                                                >  2:155728/1‑149 (MQ=255)
                                                                        cATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAg                                                                >  5:194225/1‑149 (MQ=255)
                                                                          ttGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGcc                                                                                                           >  1:76212/1‑104 (MQ=255)
                                                                          ttGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGcc                                                                                                           <  2:76212/104‑1 (MQ=255)
                                                                              ccGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                                          <  8:174313/149‑1 (MQ=255)
                                                                                gcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTGAGCACAACTCTGAGGGGTGTGGTGGCGCTGTGGGAGCGGGCGCGCCCGGGAATGGGGCGGTTAAGGCGCCGCCGCGGGGCCCGGGTAGGGGTCCACGCGAGCCCAAAGc                                                        >  2:190425/1‑149 (MQ=16)
                                                                                 cgtcgtCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCt                                                       >  1:280990/1‑149 (MQ=255)
                                                                                     gtcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                                     >  5:11537/1‑147 (MQ=255)
                                                                                      tcTGCTTCGAGGGCATGCCCGCTCCCACAGGTTCACCACAAGACTGAGGGTTGAGCTGGAGCTGTGGGAGCGGGCGCGCCCGCGAACAGGACGTATCTGCAGCCCCGGCGGGCCCCGGTTACAAGTCGACGCAAGCAGCAACCTGATAg                                                  >  8:109885/1‑149 (MQ=25)
                                                                                                        ccGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAg                                                  >  1:268831/1‑131 (MQ=255)
                                                                                                              ccACCGGTGCAGCACACTGCTGCGGGTTGGGCTGGACGTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCa                            <  3:216060/147‑1 (MQ=255)
                                                                                                              ccACAGGTTCAGCACATGGCTGCGGGTTGTGCTGGACCTGTGGTAGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCa                            <  4:144129/147‑1 (MQ=255)
                                                                                                              ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                          >  1:246654/1‑149 (MQ=255)
                                                                                                              ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                          <  3:62293/149‑1 (MQ=255)
                                                                                                                 caGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                       >  2:279010/1‑149 (MQ=255)
                                                                                                                  aGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCt                      >  5:194895/1‑149 (MQ=255)
                                                                                                                    gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                    >  8:130713/4‑149 (MQ=37)
                                                                                                                    gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGCGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                          >  1:214379/4‑143 (MQ=25)
                                                                                                                    gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGCGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                          <  2:214379/140‑1 (MQ=25)
                                                                                                                    gTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATa                                                   <  7:218875/118‑1 (MQ=255)
                                                                                                                    gTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATa                                                   >  8:218875/1‑118 (MQ=255)
                                                                                                                          cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                    <  5:120961/143‑1 (MQ=255)
                                                                                                                          cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                    >  6:120961/1‑143 (MQ=255)
                                                                                                                          cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt              >  6:100427/1‑149 (MQ=255)
                                                                                                                               gTCTGAGGGTTGTGCTGGACGTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGgcgc                                       <  5:166695/119‑1 (MQ=255)
                                                                                                                               gTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGgcgc                                       >  6:166695/1‑119 (MQ=255)
                                                                                                                                      ggTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCag  <  6:194895/149‑1 (MQ=255)
                                                                                                                                      ggTTGGGCTGGACCGGTGGTAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCag  <  6:194225/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                    
TCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAG  >  NC_002947/2958672‑2958954

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