Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,058,264:1 +T 13.5% intergenic (‑190/+30) PP_0914 ← / ← sodB GGDEF domain‑containing protein/superoxide dismutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,058,2641.T13.5% 114.6 / 16.1 37intergenic (‑190/+30)PP_0914/sodBGGDEF domain‑containing protein/superoxide dismutase
Reads supporting (aligned to +/- strand):  ref base . (15/17);  new base T (3/2);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATG  >  NC_002947/1058120‑1058399
                                                                                                                                                 |                                                                                                                                       
aaGTGTATCGGTTGCCTGACCACCCAACTTGAGTGCGCGGAATCACCGTGCCGGATGTTGGATGGTATCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCgg                                                                                                                                    <  2:571139/149‑1 (MQ=255)
       tCNGTTGCCTGACCCCGCAACTTGAATNCGCTGTATCAATNTTCCAGNGNTTTGATGTTAGCGTGGTAGGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtctc                                                                                                                             >  1:496237/1‑150 (MQ=255)
       tCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTtctc                                                                                                                             >  2:377915/1‑150 (MQ=255)
          gTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTACCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGtt                                                                                                                          <  1:151631/150‑1 (MQ=255)
             gCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTATGATGTTAGCGTAGAATGCGGCTTAAGCAATCAAAACTGTGCCAGCCTGTTTGCGGCAAAACAAAAAAAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCaa                                                                                                                       <  2:180367/150‑1 (MQ=255)
              ccTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCCATGTGCCGGGtt                                                                                                                                 >  1:325946/1‑139 (MQ=255)
              ccTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCCATGTGCCGGGtt                                                                                                                                 <  2:325946/139‑1 (MQ=255)
               cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAg                                                                                                                      <  1:541881/149‑1 (MQ=255)
               cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGGTTCTCGTTCAAg                                                                                                                      >  2:541881/1‑149 (MQ=255)
                       gCAACTTGAACGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGTCCGAATGa                                                                                                             <  1:358455/150‑1 (MQ=255)
                           cTTGAATGCGCGGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCa                                                                                                         <  2:496237/150‑1 (MQ=255)
                            ttGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAg                                                                                                        <  2:430808/150‑1 (MQ=255)
                                       gTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGt                                                                                             <  1:271685/150‑1 (MQ=255)
                                          tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                          <  1:377915/150‑1 (MQ=255)
                                                  ccAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAGGCCAAACCGTGACCGCATGTTTACGGAAACCCACAAACGCGAGAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGa                                                                                  >  1:96870/1‑150 (MQ=255)
                                                    agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGttt                                                                                                                                >  1:594978/1‑102 (MQ=255)
                                                    agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGttt                                                                                                                                <  2:594978/102‑1 (MQ=255)
                                                    agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCgttg                                                                                >  2:91005/1‑147 (MQ=255)
                                                     gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTg                                                                                         >  1:596545/1‑140 (MQ=255)
                                                     gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTg                                                                                         <  2:596545/140‑1 (MQ=255)
                                                         tttGATGTTAGCGTGGTATGCGGCCTAAGCAAGCAAAAACGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc                                                                           >  1:36957/1‑150 (MQ=255)
                                                             aTGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCa                                                                       >  1:368152/1‑150 (MQ=255)
                                                                           tGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGtt                                                         >  1:720274/1‑150 (MQ=255)
                                                                            gCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTg                                                        >  2:94341/1‑150 (MQ=255)
                                                                                     gcaagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGTCACAAAAGAACCCGGCACATTGCCGGGTTTTTCGTTAAAGACCGCAGGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGACCAGGTGATCAGGTTc                                               >  1:578175/1‑150 (MQ=255)
                                                                                      caagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc                                              <  2:368152/150‑1 (MQ=255)
                                                                                       aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa                                             >  1:40128/1‑150 (MQ=255)
                                                                                                                 ggCAACACAGGCACAAAATAACACGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTGGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCa                   <  2:96870/150‑1 (MQ=255)
                                                                                                                   cAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAgg                 <  2:40128/150‑1 (MQ=255)
                                                                                                                             acaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                          <  1:602239/67‑1 (MQ=255)
                                                                                                                             acaAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt                                                                                          >  2:602239/1‑67 (MQ=255)
                                                                                                                               cgagaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg     >  2:598232/5‑150 (MQ=255)
                                                                                                                                gagaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGa    >  1:390959/4‑150 (MQ=255)
                                                                                                                                 aaaaaCCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt   <  2:36957/150‑1 (MQ=255)
                                                                                                                                  gaaaCCCGGCCATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATg  <  1:598232/149‑1 (MQ=255)
                                                                                                                                            acaTGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCa                                                            >  2:621732/1‑82 (MQ=255)
                                                                                                                                            acaTGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCa                                                            <  1:621732/82‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
CAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATGGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATG  >  NC_002947/1058120‑1058399

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