Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,372,499 A→G 14.9% intergenic (‑1040/‑532) PP_5636 ← / → PP_3849 hypothetical protein/hemolysin‑type calcium‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,372,4990AG14.9% 53.2 / 4.2 27intergenic (‑1040/‑532)PP_5636/PP_3849hypothetical protein/hemolysin‑type calcium‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (12/11);  new base G (2/2);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

TCGTAAATGGCAGGTGAGGTCTGATTACCTCAGAGAATTACCGTCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCT  >  NC_002947/4372351‑4372636
                                                                                                                                                    |                                                                                                                                         
tCGTAAATGGCAGGTGAGGTCTGATTACCTCAGAGAATTACCGTCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTAt                                                                                                                                          >  1:686919/1‑150 (MQ=255)
tCGTAAATGGCAGGTGAGGTCTGATTACCTCAGAGAATTACCGTCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTAt                                                                                                                                          <  1:48596/150‑1 (MQ=255)
                            cTCAGAGAATTACCGTCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCa                                                                                                              >  1:597495/1‑150 (MQ=255)
                                     ttACCGTCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTAcc                                                                                                     >  2:738754/1‑150 (MQ=255)
                                           tCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCAccc                                                                                               >  2:174401/1‑150 (MQ=255)
                                               gAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGa                                                                                                                         <  2:440100/120‑1 (MQ=255)
                                               gAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGa                                                                                                                         >  1:440100/1‑120 (MQ=255)
                                                       ccATTATGTTTAAAATCGACATGAGCCCGAGACCTTCTGCGTCGCATATGTTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGtcatca                                                                                   <  2:686919/150‑1 (MQ=255)
                                                                    aaTCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCt                                                                      <  2:597495/150‑1 (MQ=255)
                                                                                cccGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATc                                                          >  2:76875/1‑150 (MQ=255)
                                                                                    cgcCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTc                                                      <  2:329591/150‑1 (MQ=255)
                                                                                         ttCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTc                                                                         <  1:516976/126‑1 (MQ=255)
                                                                                         ttCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTc                                                                         >  2:516976/1‑126 (MQ=255)
                                                                                                    acaAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGt                                      >  1:378664/1‑150 (MQ=255)
                                                                                                    acaAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGt                                      >  1:282474/1‑150 (MQ=255)
                                                                                                      aaaCTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTGTGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCt                                    >  2:266621/1‑150 (MQ=255)
                                                                                                      aaaCTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTGTGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCt                                    >  2:266693/1‑150 (MQ=255)
                                                                                                            tGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAAt                              <  2:282474/150‑1 (MQ=255)
                                                                                                              gCGAGCATCAGTCCAGCGTATACGCTCTAGATTTATGGATGGATAGATCATCATTGACCCACGCTCCAGACACTACCGCACCCAAACTGTCATACCTAACAATTCCTTCAATAGTAAATCTTTCGTCAGGCCGGTCTGGTCTTACAATAg                            >  1:581292/1‑150 (MQ=255)
                                                                                                                  gCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACgg                        <  1:76875/150‑1 (MQ=255)
                                                                                                                                  taCGATCTAGATTAATTTGTGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGAt        <  1:266693/150‑1 (MQ=255)
                                                                                                                                  taCGATCTAGATTAATTTGTGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGAt        <  1:266621/150‑1 (MQ=255)
                                                                                                                                  taCGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGAt        >  2:636307/1‑150 (MQ=255)
                                                                                                                                        cTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCt  <  1:738754/150‑1 (MQ=255)
                                                                                                                                        cTAGATTAATTTATGGATAAATCATCAATGAGACACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCTGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCt  <  1:175133/150‑1 (MQ=255)
                                                                                                                                            aTTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAgggcgggc                                          >  2:435786/1‑106 (MQ=255)
                                                                                                                                            aTTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAgggcgggc                                          <  1:435786/106‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
TCGTAAATGGCAGGTGAGGTCTGATTACCTCAGAGAATTACCGTCCTGAACTGCTCCAGTATGTTGAAAATCGACATGAGCCCGCGCCATTCTGAGTAGCACAAACTTTGGCGAGCATCAGTCCAGCGTATACGATCTAGATTAATTTATGGATAAATCATCAATGAGCCACGCTCCAGACACTACCGCACCCAAACTGTCATCACTAACAATTCCCTCGATAGTAAATCTTTCGTCAGGGCGGGCTGGTCTTACAATAGACGGCCAAGCCACGGACGATTTATCT  >  NC_002947/4372351‑4372636

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