Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,831,538 T→C 14.4% intergenic (‑8/+131) PP_2484 ← / ← PP_2485 ArsR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,831,5380TC14.4% 64.5 / 3.4 30intergenic (‑8/+131)PP_2484/PP_2485ArsR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/13);  new base C (2/2);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

ATACTGACGCACACACCTACCGTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTATTGTCCTGTCAGGTAGAACTATCGG  >  NC_002947/2831391‑2831685
                                                                                                                                                   |                                                                                                                                                   
atACTTACGCACTCACCTACCGTTTCTGGGTCGACCTCATGCTCAGGAAAGACTGTCCTGGTATCCTTGAGCAAATTTAGGATGGCAAGGCGCATCTGTTTCGCCAATGCTTTCAGGGCTTAGTTGAAGTCAATAGTCATGAAttttt                                                                                                                                                     <  3:233042/148‑1 (MQ=255)
    tGACGCACACACCTACCGTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTCTGATa                                                                                                                                                >  2:65946/1‑149 (MQ=255)
    tGACGCACACACCTACCGTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTCTGATa                                                                                                                                                >  2:69408/1‑149 (MQ=255)
               ccTACCGTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTa                                                                                                                                     <  5:45915/149‑1 (MQ=255)
                     gTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGa                                                                                                                                  >  3:2496/1‑146 (MQ=255)
                     gTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGa                                                                                                                                  <  4:2496/146‑1 (MQ=255)
                        tCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCt                                                                                                                             <  2:182550/148‑1 (MQ=255)
                          tGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTTGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATAc                                                                                                                                               <  5:177945/128‑1 (MQ=255)
                          tGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATAc                                                                                                                                               >  6:177945/1‑128 (MQ=255)
                              tCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATc                                                                                                                      >  6:220210/1‑149 (MQ=255)
                                  cctcctGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTcc                                                                                                                  >  5:198168/1‑149 (MQ=255)
                                         cTCAGGAAAGACTGTCCTTGGATCCTTGAGCCAATTCAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAAATTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTc                                                                                                            <  4:61227/148‑1 (MQ=255)
                                            aGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATc                                                                                                                      >  3:113843/1‑135 (MQ=255)
                                            aGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATc                                                                                                                      <  4:113843/135‑1 (MQ=255)
                                            aGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGAGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTAGGATATCTAGATATCTTCATGGTTcgcg                                                                                                         >  8:174696/1‑148 (MQ=255)
                                                       tgcTGGTATCCTTGAGGCTATTGAGGATGGCAAGGCGGATCGTGTTCGCCAATGCTTTCTGGGCTTCGTTGACGTCAATAGTCTTTAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAAtt                                                                                                <  8:70374/144‑1 (MQ=255)
                                                       tCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTg                                                                                               <  3:119855/147‑1 (MQ=255)
                                                         cTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCAAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGTAAg                                                                                           <  5:220210/149‑1 (MQ=255)
                                                                        cAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGAtc                                                                             >  7:73764/1‑148 (MQ=255)
                                                                                                 ggTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGAc                                                   >  7:78121/1‑149 (MQ=255)
                                                                                                       ccAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTCTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGActggctg                                             <  1:69408/149‑1 (MQ=255)
                                                                                                       ccAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTCTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGActggctg                                             <  1:65946/149‑1 (MQ=255)
                                                                                                                           ttGACGTCAATAGTCATGAATTTTTTTATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTctgga                                      <  7:7184/136‑1 (MQ=255)
                                                                                                                           ttGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTctgga                                      >  8:7184/1‑136 (MQ=255)
                                                                                                                           ttGACGTCAATAGTCATGAATTTTTTGATACCTGTTGTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTAtt                         <  6:198168/149‑1 (MQ=255)
                                                                                                                                     tAGTCATGAATTTTATGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTATTGTCCTGTCac               <  7:174696/149‑2 (MQ=255)
                                                                                                                                          aTGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTa                           <  7:173342/132‑1 (MQ=255)
                                                                                                                                          aTGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTa                           >  8:173342/1‑132 (MQ=255)
                                                                                                                                           tGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTATTGTCCTGTCAGGTAGaa         >  8:135266/1‑149 (MQ=255)
                                                                                                                                              aTTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTct                                         >  6:201116/1‑114 (MQ=255)
                                                                                                                                              aTTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTct                                         <  5:201116/114‑1 (MQ=255)
                                                                                                                                                tttttGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTATTGTCCTGTCAGGTAGAACTATc    >  7:94617/1‑149 (MQ=255)
                                                                                                                                                  tttGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTATTGTCCTGTCAGGTAGAACTATCgg  >  4:190065/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATACTGACGCACACACCTACCGTTTCTGGGTCGACCTCCTGCTCAGGAAAGACTGTCCTGGGATCCTTGAGCCAATTGAGGATGGCAAGGCGCATCGGGTTCGCCAATGCTTTCAGGGCTTCGTTGACGTCAATAGTCATGAATTTTTTGATACCTGTTTTTTACGATATCTAGATATCTTCCTGGTTCGCGCATCAAATTGGAAGGCTTCGCCGTGATCTCTGCGCGCCTGCCCGTCGCAATGACTGGCTGGTCTGGACTGGAATCCTATTGTCCTGTCAGGTAGAACTATCGG  >  NC_002947/2831391‑2831685

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