Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,281,219 C→T 10.5% intergenic (+39/+4) PP_5742 → / ← PP_1118 membrane protein/recombinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,281,2190CT10.5% 92.4 / 3.7 39intergenic (+39/+4)PP_5742/PP_1118membrane protein/recombinase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base T (2/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

CTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTC  >  NC_002947/1281072‑1281345
                                                                                                                                                   |                                                                                                                              
cTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCAc                                                                                                                              <  2:270753/150‑1 (MQ=255)
 tgtgATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACAAAACTGAATTAAGAGAGCCTGCCTTAAGGCGGGTCTATTATGAGTCACg                                                                                                                             >  2:354045/1‑150 (MQ=255)
    gATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTa                                                                                                                          <  2:433582/150‑1 (MQ=255)
       aaTTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCCTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTAttt                                                                                                                       <  2:719476/150‑1 (MQ=255)
               tgttgCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAatat                                                                                                               <  1:586852/150‑1 (MQ=255)
               tgttgCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAatat                                                                                                               <  1:399766/150‑1 (MQ=255)
                   gCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCt                                                                                                           <  1:646132/150‑1 (MQ=255)
                          aCTATGTGTTTTTGACAGAGATGGGGCGGGAGTGTGTCTTGCTAACACTGAAGCCGCTTATATTGGAATCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATCAGTCCCGTTATTTGATAATATGTCTGGCTGTc                                                                                                    <  2:53190/150‑1 (MQ=255)
                                  tttttGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGc                                                                                            >  1:347478/1‑150 (MQ=255)
                                              aTGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAgg                                                                                <  2:632568/150‑1 (MQ=255)
                                              aTGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAgg                                                                                >  1:632568/1‑150 (MQ=255)
                                              aTGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCAAATTATTTGATAATATGTCTGGCTTTCCACATGTTTAGCATATGAgg                                                                                >  2:579107/1‑150 (MQ=255)
                                               tGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGc                                                                               >  2:324146/1‑150 (MQ=255)
                                                         gTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGt                                                                     <  2:223771/150‑1 (MQ=255)
                                                            tGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGGGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTccc                                                                  >  2:634373/1‑150 (MQ=255)
                                                                ttGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCAtgt                                                                                    <  1:154586/128‑1 (MQ=255)
                                                                ttGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCAtgt                                                                                    >  2:154586/1‑128 (MQ=255)
                                                                    aTCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCa                                                          >  1:658548/1‑150 (MQ=255)
                                                                        acaGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGttt                                                      >  2:655253/1‑150 (MQ=255)
                                                                               ccGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTTACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCaaa                                               >  1:808870/1‑150 (MQ=255)
                                                                               ccGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTTACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCaaa                                               >  1:386290/1‑150 (MQ=255)
                                                                                cGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTg                                                                      <  2:51997/126‑1 (MQ=255)
                                                                                cGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTg                                                                      >  1:51997/1‑126 (MQ=255)
                                                                                                  gACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTc                            <  2:347478/150‑1 (MQ=255)
                                                                                                    ccAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAg                          <  1:324146/150‑1 (MQ=255)
                                                                                                     cAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAgg                         >  2:788458/1‑150 (MQ=255)
                                                                                                       aaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGc                                                                                            <  2:334870/81‑1 (MQ=255)
                                                                                                       aaCTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGc                                                                                            >  1:334870/1‑81 (MQ=255)
                                                                                                              ttAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTTACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                <  2:386290/150‑1 (MQ=255)
                                                                                                              ttAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTTACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                <  2:808870/150‑1 (MQ=255)
                                                                                                              ttAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                <  1:655253/150‑1 (MQ=255)
                                                                                                              ttAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTgg                <  2:658548/150‑1 (MQ=255)
                                                                                                                    gagCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAAc          >  1:337535/1‑150 (MQ=255)
                                                                                                                    gagCCCGCCATAAGGCGGGTCTATTATGAGTAACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTg                 >  1:761469/1‑143 (MQ=255)
                                                                                                                    gagCCCGCCATAAGGCGGGTCTATTATGAGTAACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTg                 <  2:761469/143‑1 (MQ=255)
                                                                                                                        ccGCCATAAGGCGGGTCTATTATGAGTCACGTTATATGATAATAACTTTCGCATTTCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTGCAGGTCGCCTTGGGTCAACTGaa      >  1:363320/1‑150 (MQ=255)
                                                                                                                            cATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTc  >  2:562241/1‑150 (MQ=255)
                                                                                                                                           ttatGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                    >  2:562458/1‑85 (MQ=255)
                                                                                                                                           ttatGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTg                                                    <  1:562458/85‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                              
CTGTGATAATTCTGTTGTTGCAGAAGACTATGTGTTTTTGACAGAGATGGGGCGGGAGTATGTCTTGCATCCACAGAAGCCGCTTCTATGGGAAGCTAGACCAAACTGAATTAAGAGAGCCCGCCATAAGGCGGGTCTATTATGAGTCACGTTATTTGATAATATGTCTGGCTTTCCACATGGCTAGCATGTGAGGCGGTAGACTGTCCCAAGATTCAGTTTTGGCAAAGGCGGCTTCCACCTCGTTCAGGTCGCCTTGGGTCAACTGAAGGTC  >  NC_002947/1281072‑1281345

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