Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,075,430 A→G 12.9% intergenic (+146/+174) PP_4468 → / ← PP_t58 Cro/CI family transcriptional regulator/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,075,4300AG12.9% 68.5 / 3.9 31intergenic (+146/+174)PP_4468/PP_t58Cro/CI family transcriptional regulator/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base G (2/2);  total (15/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 = 2.52e-01

ATCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGCTTCAAA  >  NC_002947/5075284‑5075577
                                                                                                                                                  |                                                                                                                                                   
aTCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGcc                                                                                                                                                  <  2:592686/150‑1 (MQ=255)
    ggCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAg                                                                                                                                              >  1:80708/1‑150 (MQ=255)
          aaGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGtt                                                                                                                                        <  1:116999/150‑1 (MQ=255)
             gaggACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAag                                                                                                                                     <  2:184356/148‑1 (MQ=255)
             gCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCATGCACATTCCAGGAAAGTGCCGGGCGAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAag                                                                                                                                     <  1:383008/150‑1 (MQ=255)
                gACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAgg                                                                                                                                  >  2:565682/1‑150 (MQ=255)
                       aaaCGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTagaa                                                                                                                           >  1:50913/1‑150 (MQ=255)
                            tAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTagaagagaa                                                                                                                      <  2:359292/150‑1 (MQ=255)
                               aTCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGcc                                                                                                                                                  <  2:579827/119‑1 (MQ=255)
                               aTCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGcc                                                                                                                                                  >  1:579827/1‑119 (MQ=255)
                                   aTTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAAt                                                                                                               >  2:281391/1‑150 (MQ=255)
                                        aCTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATg                                                                                                                                          <  1:524163/118‑1 (MQ=255)
                                        aCTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATg                                                                                                                                          >  2:524163/1‑118 (MQ=255)
                                                 aTCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCAt                                                                                                 <  1:439918/150‑1 (MQ=255)
                                                   cAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAg                                                                                               <  1:565682/150‑1 (MQ=255)
                                                         aaCGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAg                                                                                         >  1:246805/1‑150 (MQ=255)
                                                           cGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGaa                                                                                          >  2:273972/1‑147 (MQ=255)
                                                           cGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGaa                                                                                          <  1:273972/147‑1 (MQ=255)
                                                                                                     ttGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCg                                             <  2:80708/150‑1 (MQ=255)
                                                                                                         cAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGata                                         >  1:491949/1‑150 (MQ=255)
                                                                                                              tAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAg                                                                                               >  1:28204/1‑91 (MQ=255)
                                                                                                              tAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAg                                                                                               <  2:28204/91‑1 (MQ=255)
                                                                                                                    tGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGGAGACTAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACACAACTCAAGCAAAGCGGATATGAAGAAcaca                              >  1:111987/1‑150 (MQ=255)
                                                                                                                              gTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGgcag                    >  2:436586/1‑150 (MQ=255)
                                                                                                                                 acatgcttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCaa                 >  2:528035/7‑150 (MQ=255)
                                                                                                                                       ttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCgaga                                                                                 >  2:570134/1‑80 (MQ=255)
                                                                                                                                       ttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCgaga                                                                                 <  1:570134/80‑1 (MQ=255)
                                                                                                                                       ttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTgg           <  1:528035/150‑1 (MQ=255)
                                                                                                                                                tCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAAc                                                                         >  2:73135/1‑79 (MQ=255)
                                                                                                                                                tCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGCTTCaaa  <  2:246805/150‑1 (MQ=255)
                                                                                                                                                gcAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAAc                                                                         <  1:73135/78‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGGTGCTTCAAA  >  NC_002947/5075284‑5075577

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: