Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,890 G→A 9.7% intergenic (‑59/+127) PP_3629 ← / ← opdT‑II hypothetical protein/tyrosine‑specific outer membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,8900GA 9.7% 110.8 / 2.5 41intergenic (‑59/+127)PP_3629/opdT‑IIhypothetical protein/tyrosine‑specific outer membrane porin D
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base A (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

CTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTAC  >  NC_002947/4125746‑4126038
                                                                                                                                                |                                                                                                                                                    
cTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAg                                                                                                                                                  >  8:279504/1‑149 (MQ=255)
       ggAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGa                                                                                                                                           <  3:325207/149‑1 (MQ=255)
               cGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCAc                                                                                                                                               >  8:781/1‑137 (MQ=255)
               cGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCAc                                                                                                                                               <  7:263213/137‑1 (MQ=255)
               cGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCAc                                                                                                                                               <  7:781/137‑1 (MQ=255)
               cGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCAc                                                                                                                                               >  8:263213/1‑137 (MQ=255)
                       gAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATg                                                                                                                           >  5:107671/1‑149 (MQ=255)
                          gTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAg                                                                                                                        >  7:223107/1‑149 (MQ=255)
                          gTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAg                                                                                                                        >  1:191023/1‑149 (MQ=255)
                                      ttCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCgg                                                                                                                                                     >  6:283008/1‑108 (MQ=255)
                                      ttCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCgg                                                                                                                                                     <  5:283008/108‑1 (MQ=255)
                                           cATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCaa                                                                                                       <  8:223107/149‑1 (MQ=255)
                                                   tCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATaa                                                                                               <  8:75074/149‑1 (MQ=255)
                                                              gCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTTATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGccc                                                                                    >  5:31476/1‑149 (MQ=255)
                                                              gCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGccc                                                                                    <  7:279504/149‑1 (MQ=255)
                                                                      tgtTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCaaa                                                                                                                             >  2:280357/1‑100 (MQ=255)
                                                                      tgtTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCaaa                                                                                                                             <  1:280357/100‑1 (MQ=255)
                                                                        ttCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgc                                                                          <  3:250036/149‑1 (MQ=255)
                                                                        ttCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgc                                                                          <  4:134436/149‑1 (MQ=255)
                                                                             ggCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgcgc                                                                        >  4:351491/1‑146 (MQ=255)
                                                                              gCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAg                                                                     >  2:183106/1‑148 (MQ=255)
                                                                                 gTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGaaaa                                                                 <  3:330734/149‑1 (MQ=255)
                                                                                    aTGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATc                                                               >  5:18070/1‑148 (MQ=255)
                                                                                            gCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCAGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAAc                                                      >  4:261238/1‑149 (MQ=255)
                                                                                            gCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCAGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAAc                                                      >  4:256617/1‑149 (MQ=255)
                                                                                             cTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAAc                                                      >  7:203560/1‑148 (MQ=255)
                                                                                                 aaGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGc                                                                               >  4:327664/1‑119 (MQ=255)
                                                                                                   gATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAg                                                 >  1:36357/1‑147 (MQ=255)
                                                                                                   gATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGa                                               <  2:36357/149‑1 (MQ=255)
                                                                                                    aTCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAg                                              <  6:107671/149‑1 (MQ=255)
                                                                                                        aaGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGa                                               >  2:311526/1‑144 (MQ=255)
                                                                                                        aaGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGa                                               <  1:311526/144‑1 (MQ=255)
                                                                                                                  ttGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCGCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGcc                                 >  5:107498/1‑148 (MQ=255)
                                                                                                                     cgtcaccccccATAGCCCCGGCGTACCGGTCCCACCTGACCACCTTCTGAAAAGGCCTCATTGATCGCCCATCAAGCCTATAAAAAAAGCCCCCCAACCGTAGCGCCACAAAAACCACCTAAACCAACCGCCCAACAAGGTTGCCAAgg                             <  6:31476/142‑1 (MQ=255)
                                                                                                                           gCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGa                                                                            >  1:12164/1‑96 (MQ=255)
                                                                                                                           gCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGa                                                                            <  2:12164/96‑1 (MQ=255)
                                                                                                                           gCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAAt                       <  6:18070/149‑1 (MQ=255)
                                                                                                                              ccATAGCCCCGGCGTACCGGTAGCACCCGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGcaacaa                                        <  8:140310/129‑1 (MQ=255)
                                                                                                                              ccATAGCCCCGGCGTACCAGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATcac                    <  3:261238/149‑1 (MQ=255)
                                                                                                                              ccATAGCCCCGGCGTACCAGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATcac                    <  3:256617/149‑1 (MQ=255)
                                                                                                                                    ccccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATc              <  3:96785/149‑1 (MQ=255)
                                                                                                                                         gCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGcaacaa                                        >  8:27050/1‑118 (MQ=255)
                                                                                                                                              ccGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGt    <  3:351491/149‑1 (MQ=255)
                                                                                                                                                ggTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTAc  <  1:183106/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTAC  >  NC_002947/4125746‑4126038

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: ‑