Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,075,425 C→T 14.3% intergenic (+141/+179) 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,4250CT14.3% 63.2 / 3.9 28intergenic (+141/+179)PP_4468/PP_t58Cro/CI family transcriptional regulator/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (2/2);  total (14/14)
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.04e-01

ATCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGG  >  NC_002947/5075284‑5075568
                                                                                                                                             |                                                                                                                                               
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                                                                                 <  1:273972/147‑1 (MQ=255)
                                                           cGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGaa                                                                                 >  2:273972/1‑147 (MQ=255)
                                                                                                     ttGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCg                                    <  2:80708/150‑1 (MQ=255)
                                                                                                         cAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGata                                >  1:491949/1‑150 (MQ=255)
                                                                                                              tAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAg                                                                                      <  2:28204/91‑1 (MQ=255)
                                                                                                              tAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAg                                                                                      >  1:28204/1‑91 (MQ=255)
                                                                                                                    tGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGGAGACTAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACACAACTCAAGCAAAGCGGATATGAAGAAcaca                     >  1:111987/1‑150 (MQ=255)
                                                                                                                              gTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGgcag           >  2:436586/1‑150 (MQ=255)
                                                                                                                                 acatgcttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCaa        >  2:528035/7‑150 (MQ=255)
                                                                                                                                       ttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCgaga                                                                        <  1:570134/80‑1 (MQ=255)
                                                                                                                                       ttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCgaga                                                                        >  2:570134/1‑80 (MQ=255)
                                                                                                                                       ttCGGCTGATCGGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTgg  <  1:528035/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
ATCAGGCACAAAGGCGGACAAGCAAACGTAGATCAATTCCACTCTGCCCATCAATCGAACGACATCCATTTCCCGCAGCATCTCAGAAATTTTCCAGGAAATTGCCAGGCTAGGAATGAGCTGCTGGTAAGCGAGTTCGGCCGATCAGCCGAAGCATGTTAAGAGGCGTAGAAGAGAAGCAGAATAGCCCGACACCCATAGGAGAAGCGCCGAGAGCAGGAACTTCCACCCCACTCAACTCAAGCAAAGCGGATATGAAGAACACAGAGAAGGCAGCAAACCTGG  >  NC_002947/5075284‑5075568

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: