Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,220 G→T 26.8% intergenic (‑128/‑140) PP_5175 ← / → PP_5176 HlyD family secretion protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,903,2200GT26.8% 36.6 / 10.7 26intergenic (‑128/‑140)PP_5175/PP_5176HlyD family secretion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base T (5/2);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

ATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAA  >  NC_002947/5903091‑5903360
                                                                                                                                 |                                                                                                                                             
aTGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGTGCGAGCCGGCGAAGAGGCCTCca                                                                                                                           >  2:486277/1‑150 (MQ=255)
          cttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTa                                                                                                                 <  2:353136/150‑1 (MQ=255)
           ttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTAt                                                                                                                >  1:735707/1‑150 (MQ=255)
                   gCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCAATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCt                                                                                                        >  1:204223/1‑150 (MQ=255)
                     acaAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGc                                                                                                      <  1:701160/150‑1 (MQ=255)
                                                    tAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATc                                                                      >  2:88162/1‑150 (MQ=255)
                                                       ttCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGACGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATcccc                                                                   >  2:107706/1‑150 (MQ=255)
                                                                        aTCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCtt                                                  <  1:88162/150‑1 (MQ=255)
                                                                             tGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATAATCATACACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTAGGTTTcaca                                             >  2:621752/1‑150 (MQ=255)
                                                                                 gggggCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCAGCGCAGATGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGa                                         <  1:407539/150‑1 (MQ=255)
                                                                                      ccATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAggg                                          >  1:652505/1‑144 (MQ=255)
                                                                                      ccATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAggg                                          <  2:652505/144‑1 (MQ=255)
                                                                                          gCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAAGTGGTGCCATATCCAGCCACGACACTCAACACACTTTCGCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACAt                                >  1:178755/1‑150 (MQ=255)
                                                                                                          ccTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCaaaa                <  2:735707/150‑1 (MQ=255)
                                                                                                              cGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGcc            <  1:709978/150‑1 (MQ=255)
                                                                                                               gggggCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCg           >  1:464223/1‑150 (MQ=255)
                                                                                                                   gCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCg                                                                            <  1:81032/81‑1 (MQ=255)
                                                                                                                   gCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCg                                                                            >  2:81032/1‑81 (MQ=255)
                                                                                                                    ccTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGc      <  2:178755/150‑1 (MQ=255)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATc                                                                                                <  2:334630/57‑1 (MQ=25)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATc                                                                                                >  1:334630/1‑57 (MQ=25)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCaa                                                            >  2:562287/1‑92 (MQ=255)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCaa                                                            <  1:562287/92‑1 (MQ=255)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGaa  >  2:339789/1‑150 (MQ=255)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGaa  >  1:467736/1‑150 (MQ=255)
                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATATATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGaa  >  1:217532/1‑150 (MQ=255)
                                                                                                                                 |                                                                                                                                             
ATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAA  >  NC_002947/5903091‑5903360

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: