Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,927 T→C 8.3% intergenic (‑96/+90) 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,9270TC 8.3% 131.6 / 2.7 48intergenic (‑96/+90)PP_3629/opdT‑IIhypothetical protein/tyrosine‑specific outer membrane porin D
Reads supporting (aligned to +/- strand):  ref base T (24/20);  new base C (2/2);  total (26/22)
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.66e-01

TCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACTTCGTCCAGATCCGACTC  >  NC_002947/4125782‑4126073
                                                                                                                                                 |                                                                                                                                                  
tctTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCGCCAAGATCGAAGGTGAGCCGTGAGTCTAGCCCCATAGCCCCGGCGAACCACATCCACCGGACAACCTGATGCAAATGCAGGAATGATAgg                                                                                                                                                 >  1:302402/1‑149 (MQ=255)
     ggCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATc                                                                                                                                            <  1:191853/149‑1 (MQ=255)
      gCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTTATAGGCGATc                                                                                                                                            >  7:314268/1‑148 (MQ=255)
           gTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGCTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCCGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGcc                                                                                                                                       >  5:8854/1‑148 (MQ=255)
              gTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCGTGCCACCTGATGCAAATGCAGGATTGACAGGCGATCAAGCCGAt                                                                                                                                    >  1:120690/1‑148 (MQ=255)
                 gTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATaaa                                                                                                                                 >  3:308419/1‑148 (MQ=255)
                    ccACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATaaaaaaa                                                                                                                             <  4:308419/149‑1 (MQ=255)
                    ccACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATaaaaaaa                                                                                                                             <  2:184940/149‑1 (MQ=255)
                       caGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAgcg                                                                                                                          >  2:279501/1‑149 (MQ=255)
                        aGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAgcgc                                                                                                                         <  8:13288/149‑1 (MQ=255)
                                       ttGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATa                                                                                                                                                   <  5:86108/108‑1 (MQ=255)
                                       ttGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATa                                                                                                                                                   >  6:86108/1‑108 (MQ=255)
                                         ggCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgcgc                                                                                                           <  6:259032/146‑1 (MQ=255)
                                                   gTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCTGAAAATCCAcc                                                                                              >  5:53283/1‑149 (MQ=255)
                                                    tCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCAcc                                                                                              >  4:338689/1‑148 (MQ=255)
                                                         cTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAAc                                                                                         <  1:279501/148‑1 (MQ=255)
                                                          tCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAAcc                                                                                        >  1:273031/1‑148 (MQ=255)
                                                          tCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAAcc                                                                                        >  4:34291/1‑148 (MQ=255)
                                                            cAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGaaa                                                                                                     >  2:128242/1‑133 (MQ=255)
                                                            cAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGaaa                                                                                                     <  1:128242/133‑1 (MQ=255)
                                                            cAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGACAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCTCCTAAACCa                                                                                       >  1:355944/1‑147 (MQ=255)
                                                                aTCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAg                                                                                 >  3:73849/1‑149 (MQ=255)
                                                                       gTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGACAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAg                                                                          <  2:355944/149‑1 (MQ=255)
                                                                         gAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAg                                                                                 <  7:307735/140‑1 (MQ=255)
                                                                         gAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAg                                                                                 >  8:307735/1‑140 (MQ=255)
                                                                          aGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGtt                                                                       >  2:83/1‑149 (MQ=255)
                                                                                  gTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGt                                                               <  2:273031/149‑1 (MQ=255)
                                                                                   tCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTg                                                              >  3:75530/1‑149 (MQ=255)
                                                                                    cTAGCCCCATAGCCCCGGCGTACCTGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATGAAGCCGATAAAAAAAGCGCCCCAAGCTGAGCGCCAGAAAATCCACCTAAACCAAAg                                                                                    <  7:202377/126‑1 (MQ=255)
                                                                                    cTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAACATCAACCTAAACCAAAg                                                                                    >  8:202377/1‑126 (MQ=255)
                                                                                       gCCGCATAGCCCCGGCGTACCGGTAGCACCGGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGa                                                            <  6:8854/147‑1 (MQ=255)
                                                                                          ccATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATca                                                        >  2:292691/1‑148 (MQ=255)
                                                                                           cATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGACAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATcaca                                                      <  2:120690/149‑1 (MQ=255)
                                                                                                ccccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgc                                                                                                             >  8:215740/1‑89 (MQ=255)
                                                                                                ccccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgc                                                                                                             <  7:215740/89‑1 (MQ=255)
                                                                                                ccccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATc                                                 <  8:314268/149‑1 (MQ=255)
                                                                                                ccccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATc                                                 >  3:172558/1‑149 (MQ=255)
                                                                                                 cccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATc                                                 <  3:338689/148‑1 (MQ=255)
                                                                                                  ccGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGa                                               <  1:83/149‑1 (MQ=255)
                                                                                                                   ccTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGa                              <  4:75530/149‑1 (MQ=255)
                                                                                                                            cTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGCATCGACAACGTGTACTCCACGATCATCCGCa                     >  6:94742/1‑149 (MQ=255)
                                                                                                                              gATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACt                   >  4:123835/1‑149 (MQ=255)
                                                                                                                               aTGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGc                      >  2:214663/1‑145 (MQ=255)
                                                                                                                               aTGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGc                      <  1:214663/145‑1 (MQ=255)
                                                                                                                                            aTTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACTTCGTCCAGATCCGa     >  5:15929/1‑149 (MQ=255)
                                                                                                                                             ttGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACTTCGTCCAGATCCGAc    <  3:34291/149‑1 (MQ=255)
                                                                                                                                               gATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACTTCGTCCAGATCCGACTc  >  2:320298/1‑149 (MQ=255)
                                                                                                                                                aTAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCTGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACTTCGTCCAGATCCGACTc  <  6:53283/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGGAGCAACAAGGTTGCCAAGGTGGAATCACAGAATCGACAACGGGTACTCCACGATCATCCGCACTTCGTCCAGATCCGACTC  >  NC_002947/4125782‑4126073

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