Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 413,107 C→G 11.0% A686A (GCC→GCG glnE → glutamate‑ammonia‑ligase adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947413,1070CG11.0% 95.8 / 2.4 34A686A (GCC→GCGglnEglutamate‑ammonia‑ligase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base G (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC  >  NC_002947/412959‑413256
                                                                                                                                                    |                                                                                                                                                     
gcCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCAACGTAAATCAGCGGTAACCTGCAGTTGATGAAAGTCAGCGTCTACCCGACCGCGCCGTCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGGCCGTGGccc                                                                                                                                                      >  1:408842/1‑150 (MQ=255)
     ttCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCAc                                                                                                                                                 <  1:515126/150‑1 (MQ=255)
        aaGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGc                                                                                                                                              <  2:31104/150‑1 (MQ=255)
         aGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtc                                                                                                                                             <  2:762332/150‑3 (MQ=255)
          gCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtct                                                                                                                                            <  1:335529/150‑4 (MQ=255)
           cTGGCCCACAGCCTGCGCGTCGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCaa                                                                                                                                           <  1:838096/150‑1 (MQ=255)
                       cTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGc                                                                                                                               <  2:565772/150‑1 (MQ=255)
                                    ccTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCtgt                                                                                                                  >  2:505350/1‑150 (MQ=255)
                                                  ggCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTGtc                                                                                                    >  2:175186/1‑150 (MQ=255)
                                                           ccGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATc                                                                                           <  1:5607/150‑1 (MQ=255)
                                                            cGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCg                                                                                          <  1:6072/150‑1 (MQ=255)
                                                                 atgaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTAc                                                                                     >  2:575083/1‑150 (MQ=255)
                                                                   gaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTTGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACgg                                                                                   >  2:504033/1‑150 (MQ=255)
                                                                                  cctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACgg                                                                                   >  2:126039/1‑135 (MQ=255)
                                                                                  cctgacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACgg                                                                                   <  1:126039/135‑1 (MQ=255)
                                                                                          ggctggctGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTTGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg                                                            <  1:266547/150‑1 (MQ=255)
                                                                                                     gCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCg                                                 >  1:339395/1‑150 (MQ=255)
                                                                                                               aCCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACGCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctggacc                                       >  1:497186/1‑150 (MQ=255)
                                                                                                                ccAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctggacct                                      >  2:769910/1‑150 (MQ=255)
                                                                                                                  aGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGacctggacctgg                                    <  1:81496/150‑1 (MQ=255)
                                                                                                                   ggTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGt                                   >  1:839031/1‑150 (MQ=255)
                                                                                                                   ggTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGt                                   <  1:645193/150‑1 (MQ=255)
                                                                                                                        tggccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGTGCGTCCTGGAGCGCGGCCACGGTTCGGACCTGGACCTGGTCTTCa                              <  1:575083/150‑1 (MQ=255)
                                                                                                                          gccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATc                            >  2:410839/1‑150 (MQ=255)
                                                                                                                            cctggcctggCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATggg                                                                          >  2:663811/1‑102 (MQ=255)
                                                                                                                            cctggcctggCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATggg                                                                          <  1:663811/102‑1 (MQ=255)
                                                                                                                             gaggcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCAt                         >  2:549464/3‑150 (MQ=255)
                                                                                                                             ctggcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCCGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCAt                         <  1:175186/150‑1 (MQ=255)
                                                                                                                                gcctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGAc                      <  2:408842/150‑1 (MQ=255)
                                                                                                                                       gcCAGACCGTGGCGCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATc               >  2:527090/1‑150 (MQ=255)
                                                                                                                                       gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATc               >  1:635532/1‑150 (MQ=255)
                                                                                                                                       gcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATc               <  1:505350/150‑1 (MQ=255)
                                                                                                                                         ttggccGTGGCGGGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTCCGGGGAGCTGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCg             >  1:148766/5‑150 (MQ=255)
                                                                                                                                          aGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGc            >  1:102567/1‑150 (MQ=255)
                                                                                                                                              cGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGCCGGCGATCCGCAGGc        >  2:146861/1‑150 (MQ=255)
                                                                                                                                                   cccGCCACGGCGAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGCTGGGCGGCCTGGAGCTCGGCGACGGTACGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGaaa   <  2:379383/150‑1 (MQ=255)
                                                                                                                                                    ccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAc  <  1:758378/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCATTTCAAGCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC  >  NC_002947/412959‑413256

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: