Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,224,546 T→G 16.1% intergenic (‑614/+207) PP_3703 ← / ← PP_3704 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,224,5460TG16.1% 63.6 / 9.3 31intergenic (‑614/+207)PP_3703/PP_3704hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base G (2/3);  total (13/18)
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.40e-01

GTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGA  >  NC_002947/4224397‑4224692
                                                                                                                                                       |                                                                                                                                                  
gtcgCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACGACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGTAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTCAGTGGC‑‑CCAGCAGTGCCCCt                                                                                                                                                    >  2:578264/1‑150 (MQ=255)
gtcgCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCt                                                                                                                                                    <  2:719692/150‑1 (MQ=255)
  ctcccTATCAATCCAAACGTAAATGTCCCTAAGCCTGCTTGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCtgt                                                                                                                                                  <  1:209413/148‑1 (MQ=255)
     ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGACAACCCAGCGCCAGAGACACCTAGGAGGC‑‑CCAGCAGTGCACCtgtgag                                                                                                                                               >  2:209413/1‑148 (MQ=255)
     ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCa                                                                                                                                               <  1:110271/150‑1 (MQ=255)
           aaTCTAAACGTAAATTTCCTTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATATGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCActgctg                                                                                                                                         <  1:371348/150‑1 (MQ=255)
             tCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTg                                                                                                                                            <  2:581353/145‑1 (MQ=255)
             tCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTg                                                                                                                                            >  1:581353/1‑145 (MQ=255)
                aaaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCACTGCTgg                                                                                                                                        <  2:364068/146‑1 (MQ=255)
                aaaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCACTGATgg                                                                                                                                        >  1:364068/1‑146 (MQ=255)
                   gactaaTGGCTCTAGGCCTGACCCCGGGGAACAACCTCCGGCCTGCAACAAACAGGGCGAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCAc                                                                                                                                 <  2:351301/146‑1 (MQ=255)
                              cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCt                                                                                                                      >  1:357274/1‑150 (MQ=255)
                                aaGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTa                                                                                                                    <  2:301045/150‑1 (MQ=255)
                                 aGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                   <  1:95762/150‑1 (MQ=255)
                                  gCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCAACCAGCAATGCACAGGTGCACTGCt                                                                                                                                          >  1:162537/1‑128 (MQ=38)
                                  gCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCAACCAGCAATGCACAGGTGCACTGCt                                                                                                                                          <  2:162537/128‑1 (MQ=38)
                                   ccTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGTCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGa                                                                                                                 <  1:652106/150‑1 (MQ=255)
                                            ggggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATCCGTACTACCAGGAGGTATGCGGGAAAACCCCGCGAAAAAGTAACCTAGGTGGG‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGAc                                                                                                        <  1:776391/150‑1 (MQ=255)
                                            ggggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGAc                                                                                                        <  1:106211/150‑1 (MQ=255)
                                                       tCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTa                                                                                             >  1:442606/1‑150 (MQ=255)
                                                           tcgcGTAACAAAGACGGCCAGATAAGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATcgc                                                                                         <  2:395528/147‑1 (MQ=255)
                                                                           gCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGTCTTTCGCGCCCTTCTGTAGCGga                                                                         <  2:272283/150‑1 (MQ=255)
                                                                           gCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                                         <  2:357274/150‑1 (MQ=255)
                                                                            ccAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGtt                                                                                                                                      <  2:144924/88‑1 (MQ=255)
                                                                            ccAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGtt                                                                                                                                      >  1:144924/1‑88 (MQ=255)
                                                                                    gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACa                                                                >  1:445950/1‑150 (MQ=255)
                                                                                      gATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGa                                                              >  2:554421/1‑150 (MQ=255)
                                                                                       aTAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGaa                                                             <  2:445950/150‑1 (MQ=255)
                                                                                                                               cTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaagga                     >  2:53266/1‑150 (MQ=255)
                                                                                                                                      gC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCt              >  1:183853/1‑150 (MQ=255)
                                                                                                                                          ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCt                                                                                                              >  2:37243/1‑52 (MQ=255)
                                                                                                                                          ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCt                                                                                                              <  1:37243/52‑1 (MQ=255)
                                                                                                                                                  gCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                                         <  1:37903/81‑1 (MQ=255)
                                                                                                                                                  gCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                                         >  2:37903/1‑81 (MQ=255)
                                                                                                                                                    aCCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGa  >  2:42983/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
GTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGA  >  NC_002947/4224397‑4224692

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