Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,812,818 C→T 10.5% intergenic (+69/‑41) rplT → / → pheS 50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,812,8180CT10.5% 92.7 / 3.8 38intergenic (+69/‑41)rplT/pheS50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (18/16);  new base T (2/2);  total (20/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 = 5.38e-01

TGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAG  >  NC_002947/2812676‑2812963
                                                                                                                                              |                                                                                                                                                 
tgtcTGATCTGGCCAGGCACGCACCAGCCGCCTATGCTGCGATCGTCGCGAACGCTACAGCCCTACTCGCCTAACTACCCGCGGCAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAg                                                                                                                                            <  2:283002/147‑1 (MQ=255)
             aGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCTGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTAtt                                                                                                                               >  2:182717/1‑150 (MQ=255)
                gAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTNAACGTCNTCGATAGGGGAANAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                            <  1:504826/150‑1 (MQ=255)
                 aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCg                                                                                                                           >  1:174183/1‑150 (MQ=255)
                                       cGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGaa                                                                                                     >  2:590535/1‑150 (MQ=255)
                                            gTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                                                                      >  1:42364/1‑128 (MQ=255)
                                            gTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                                                                      <  2:42364/128‑1 (MQ=255)
                                             tCGAGAACGCTAAAGCCGTTCTGGCTTAAGTCCCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCCAAGCTCTTCCCCTATTTTCGTATCTGGCGTGTGTACATGGAAAACCTg                                                                                               >  2:13791/1‑150 (MQ=255)
                                               gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGa                                                                                             >  2:21876/1‑150 (MQ=255)
                                                      cTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTg                                                                                      >  2:266875/1‑150 (MQ=255)
                                                       tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTg                                                                                                                     <  1:356508/118‑1 (MQ=255)
                                                       tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTg                                                                                                                     >  2:356508/1‑118 (MQ=255)
                                                               ttCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGtctc                                                                                 >  2:25222/1‑146 (MQ=255)
                                                               ttCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGtctc                                                                                 <  1:25222/146‑1 (MQ=255)
                                                                              ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg                                                              >  2:488266/1‑150 (MQ=255)
                                                                              ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg                                                              >  2:141716/1‑150 (MQ=255)
                                                                               cgcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTa                                                                       <  2:198347/140‑1 (MQ=255)
                                                                               cgcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTa                                                                       >  1:198347/1‑140 (MQ=255)
                                                                                     aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg                                                       >  2:40305/1‑150 (MQ=255)
                                                                                          cccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                            >  2:335189/1‑76 (MQ=255)
                                                                                          cccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                            <  1:335189/76‑1 (MQ=255)
                                                                                               cggTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATTCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACAt                                             <  1:13791/150‑1 (MQ=255)
                                                                                                  tGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCaa                                          <  1:141716/150‑1 (MQ=255)
                                                                                                         gcTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTATGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTg                                   >  1:47110/1‑150 (MQ=255)
                                                                                                            gCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGaa                                <  1:488266/150‑1 (MQ=255)
                                                                                                               gggCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAg                             >  2:462869/1‑150 (MQ=255)
                                                                                                                  cATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                                                                     >  2:239628/1‑107 (MQ=255)
                                                                                                                  cATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                                                                     <  1:239628/107‑1 (MQ=255)
                                                                                                                  cATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATc                          <  1:590535/150‑1 (MQ=255)
                                                                                                                                gATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTAtct            <  1:40305/150‑1 (MQ=255)
                                                                                                                                 aTAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTAtctc           <  1:462869/150‑1 (MQ=255)
                                                                                                                                     gggAAGAGCCTTCAAAGCTGTTTCACTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCa       <  2:47110/150‑1 (MQ=255)
                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCa       >  2:448474/1‑150 (MQ=255)
                                                                                                                                      ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCt                                                                                                >  2:346408/1‑60 (MQ=21)
                                                                                                                                      ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCt                                                                                                <  1:346408/60‑1 (MQ=21)
                                                                                                                                      ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa                  >  1:430051/1‑138 (MQ=255)
                                                                                                                                      ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa                  <  2:430051/138‑1 (MQ=255)
                                                                                                                                          gagCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTCTACATGGAAATCCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCaagaag  <  2:174183/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAG  >  NC_002947/2812676‑2812963

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: