Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,820,642 Δ1 bp 35.3% intergenic (+400/+112) PP_2473 → / ← PP_2474 hypothetical protein/glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,820,6420A.35.3% 31.0 / 18.4 17intergenic (+400/+112)PP_2473/PP_2474hypothetical protein/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base A (4/7);  new base . (3/3);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

TTATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAG  >  NC_002947/2820493‑2820732
                                                                                                                                                      |                                                                                          
ttATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATACCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGa                                                                                            >  2:379310/1‑150 (MQ=255)
        gggATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGt                                                                                    >  1:1778/1‑150 (MQ=255)
                cAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGtt                                                                            <  2:1778/150‑1 (MQ=255)
                               aaGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGa                                                             <  2:359437/150‑1 (MQ=255)
                                         tcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCg                                                                               <  1:426430/122‑1 (MQ=255)
                                         tcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCg                                                                               >  2:426430/1‑122 (MQ=255)
                                               gTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCgg                                                                 <  1:239474/130‑1 (MQ=255)
                                               gTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCgg                                                                 >  2:239474/1‑130 (MQ=255)
                                                                 tCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGtt                                                               <  2:48049/114‑1 (MQ=255)
                                                                 tCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGtt                                                               >  1:48049/1‑114 (MQ=255)
                                                                 tCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGctc                           <  1:379310/150‑1 (MQ=255)
                                                                                         gggCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGCGCAAAGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa   <  2:314223/150‑1 (MQ=255)
                                                                                                                                        aaaaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg  >  2:131959/1‑104 (MQ=255)
                                                                                                                                        aaaaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg  <  1:131959/104‑1 (MQ=255)
                                                                                                                                          aaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg  <  1:123158/102‑1 (MQ=255)
                                                                                                                                          aaCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg  >  2:123158/1‑102 (MQ=255)
                                                                                                                                           aCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTa                                               <  1:429136/56‑1 (MQ=25)
                                                                                                                                           aCCCGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTa                                               >  2:429136/1‑56 (MQ=25)
                                                                                                                                                      |                                                                                          
TTATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAG  >  NC_002947/2820493‑2820732

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: