Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,143,074 G→C 19.4% P68P (CCG→CCC arcD‑II → arginine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,143,0740GC19.4% 69.7 / 14.1 31P68P (CCG→CCCarcD‑IIarginine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base G (10/15);  new base C (3/3);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CGGCGGGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCA  >  NC_002947/1142939‑1143204
                                                                                                                                       |                                                                                                                                  
cggcggGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCg                                                                                                                                    <  2:302964/136‑1 (MQ=255)
   cggGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGac                                                                                                                                 <  1:366557/136‑2 (MQ=255)
               gCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCgggg                                                                                                                     <  2:141273/136‑1 (MQ=255)
                                ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCaa                                                                                                    >  1:568094/1‑136 (MQ=255)
                                ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCaa                                                                                                    >  1:292509/1‑136 (MQ=255)
                                ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGCGGTATATGCCTACGCCaa                                                                                                    >  1:561278/1‑136 (MQ=255)
                                ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  <  2:338656/106‑1 (MQ=255)
                                ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  >  1:338656/1‑106 (MQ=255)
                                ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGACCTGGATGGTGGGGTGTATGCGTATGCCaa                                                                                                    >  1:541816/1‑136 (MQ=14)
                                  gcgcgTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  >  1:565213/1‑104 (MQ=255)
                                  gcgcgTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  <  2:565213/104‑1 (MQ=255)
                                          tGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCttt                                                                                          <  1:27157/136‑1 (MQ=255)
                                          tGAGGTGGGGGCAGTGCTGATCGGCGGGGGGAGCACCGCGGCCGGCATGTTGGCCCGGGCCTTCGTGTTCCCGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCttt                                                                                          <  1:29213/136‑1 (MQ=21)
                                                 ggCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGa                                                                                                                                  >  1:25448/1‑89 (MQ=255)
                                                 ggCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGa                                                                                                                                  <  2:25448/89‑1 (MQ=255)
                                                 ggCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGt                                                                                   >  2:501546/1‑136 (MQ=255)
                                                                ggCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCg                                                                                                                                    <  2:121052/72‑1 (MQ=255)
                                                                ggCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCg                                                                                                                                    >  1:121052/1‑72 (MQ=255)
                                                                  cTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCg                                                                  <  1:88140/136‑1 (MQ=255)
                                                                   tGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCgg                                                                 <  1:104830/136‑1 (MQ=255)
                                                                     ggggATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTg                                                                                         >  1:533274/1‑110 (MQ=255)
                                                                     ggggATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTg                                                                                         <  2:533274/110‑1 (MQ=255)
                                                                      gggATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCt                                                              <  1:488825/136‑1 (MQ=255)
                                                                      gggATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCt                                                              <  1:214007/136‑1 (MQ=255)
                                                                                           ttGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCct                                         >  2:62654/1‑136 (MQ=255)
                                                                                                      cTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATggg                                                                            <  1:522356/90‑1 (MQ=255)
                                                                                                      cTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATggg                                                                            >  2:522356/1‑90 (MQ=255)
                                                                                                       ttCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACg                             >  1:12605/1‑136 (MQ=255)
                                                                                                                           cAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGctgctg         >  1:153822/1‑136 (MQ=255)
                                                                                                                            aaCCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGt        <  1:91603/136‑1 (MQ=255)
                                                                                                                             aCCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGtt       <  1:162805/136‑1 (MQ=255)
                                                                                                                                  aaGCCGGAACTGGACTCGGTGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCAGCCTGGGGCTACTGGATCAGTGCCTGGCTGGGCAACGTCGGCTACTTCGTGCTGCTGTTCAGCa  <  2:561278/136‑1 (MQ=255)
                                                                                                                                  aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  1:62654/136‑1 (MQ=255)
                                                                                                                                  aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGAGCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  2:568094/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CGGCGGGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCA  >  NC_002947/1142939‑1143204

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