Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,670,221 C→T 75.0% intergenic (+37/‑6) trmD → / → rplS tRNA (guanosine(37)‑N1)‑methyltransferase TrmD/50S ribosomal protein L19

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,670,2210CT75.0% 88.6 / 29.2 56intergenic (+37/‑6)trmD/rplStRNA (guanosine(37)‑N1)‑methyltransferase TrmD/50S ribosomal protein L19
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (20/22);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCTTGGTAGGACCTTCGAACGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTGCAGGCGTTCGAAGGC  >  NC_002947/1670076‑1670364
                                                                                                                                                 |                                                                                                                                                  
ccTTGGTAGGACCTTCGAACGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGtaca                                                                                                                                                 <  3:9389/149‑1 (MQ=255)
  ttGGTAGGACCTTCGAACGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAg                                                                                                                                                >  5:83972/1‑148 (MQ=255)
             ttCGAACGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACcaac                                                                                                                                      >  4:238766/1‑147 (MQ=255)
               cGAACGACGCTCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACcaacaa                                                                                                                                   <  1:301262/147‑1 (MQ=255)
               cGAACGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACcaacaa                                                                                                                                   <  2:169885/147‑1 (MQ=255)
                 aaCGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGAACGACAAGat                                                                                                                                >  1:113538/1‑148 (MQ=255)
                  acgacgCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGatc                                                                                                                              >  6:176852/1‑148 (MQ=255)
                     acgCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGatcatc                                                                                                                           >  5:203169/1‑148 (MQ=255)
                     acgCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATcc                                                                                                                          <  6:203169/149‑1 (MQ=255)
                      cgcgCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCca                                                                                                                         <  8:203049/149‑1 (MQ=255)
                       gcgcCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCcag                                                                                                                        >  2:110733/1‑149 (MQ=255)
                                 tctGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAg                                                                                                              >  1:315710/1‑149 (MQ=255)
                                     gAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACcaac                                                                                                                                      >  4:24999/1‑123 (MQ=255)
                                     gAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACcaac                                                                                                                                      <  3:24999/123‑1 (MQ=255)
                                      agaGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAg                                                                                                         >  2:183086/3‑149 (MQ=255)
                                       aaGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAg                                                                                                         <  3:238766/148‑1 (MQ=255)
                                            gctcgctTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGAt                                                                                                    <  6:83972/148‑1 (MQ=255)
                                            gctcgctTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGCGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGAt                                                                                                    >  8:17816/1‑148 (MQ=255)
                                                gctTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAg                                                                                                 <  3:283483/147‑1 (MQ=255)
                                                    tCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCCATCCAGCAGCTCGAAGCCGAGCAGATGAGCaa                                                                                              <  1:207036/147‑1 (MQ=255)
                                                        gagaagagaaGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGcc                                                                                                            <  3:257537/128‑1 (MQ=255)
                                                        gagaagagaaGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGcc                                                                                                            >  4:257537/1‑128 (MQ=255)
                                                        gagaagagaaGAAGCTGCGCGAGGAATATCTCAGCGAGCGGTAAGATAGTTAAAAGTATCGATTGAGCATAGCGTATACATCTTAGTAGTACATCATGCACAACAAGATCATCCAGCAGCTCGACGACGACCAGATGAGCAACGAAAta                                                                                      >  5:140038/1‑148 (MQ=255)
                                                         agaagagaagAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATc                                                                                      <  7:17816/148‑1 (MQ=255)
                                                             gagaagaAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGa                                                                                  >  7:188758/1‑148 (MQ=255)
                                                             gagaagaAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGAc                                                                                 >  8:162368/1‑149 (MQ=255)
                                                                      ctgctCGAGGAATATCTCCGCGAGCGGGACGATAGTAAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCAc                                                                         >  1:262508/1‑148 (MQ=255)
                                                                           cGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCAGAGCAGATGAGAAAGGAAATCCAGACCTTCGCACCAGGCg                                                                   >  8:288474/1‑149 (MQ=255)
                                                                            gAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGa                                                                  >  7:213083/1‑149 (MQ=255)
                                                                              ggAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGa                                                                                           >  7:303661/1‑123 (MQ=255)
                                                                              ggAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGa                                                                                           <  8:303661/123‑1 (MQ=255)
                                                                                      tcCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTc                                                        >  2:81952/1‑149 (MQ=255)
                                                                                        cgctAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTc                                                        <  7:257675/147‑1 (MQ=255)
                                                                                          cgAGGGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATGCATCTTAGGAGCACAGCATGACGTACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCAGCAGGCGACACCGTTGTCGTTc                                                    <  2:262508/149‑1 (MQ=255)
                                                                                               gggACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGtt                                               <  7:288474/149‑1 (MQ=255)
                                                                                                    gATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTaa                                             <  6:122723/146‑1 (MQ=255)
                                                                                                    gATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTaa                                             >  5:122723/1‑146 (MQ=255)
                                                                                                         ttAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTgaag                                      >  8:125528/1‑148 (MQ=255)
                                                                                                            aaCGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATTCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTgaagga                                    <  1:81952/148‑1 (MQ=255)
                                                                                                            aaCGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTgaaggaag                                  <  7:133187/149‑1 (MQ=255)
                                                                                                              cGTATCTATGGTGGATCGCGTATCAATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGGGAAACCGTAGTCGTTCAGGTTAAAGTGAAGGAAGGt                                <  1:183086/149‑1 (MQ=255)
                                                                                                              cGTATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGt                                >  3:97695/1‑149 (MQ=255)
                                                                                                                tATCGATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGa                              <  8:188758/149‑1 (MQ=255)
                                                                                                                    gATGGTGGATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGCGCGt                          >  8:195200/1‑149 (MQ=255)
                                                                                                                     aTGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGc                            >  5:254867/1‑146 (MQ=255)
                                                                                                                     aTGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGc                            <  6:254867/146‑1 (MQ=255)
                                                                                                                      tggtggATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTc                        <  7:125528/149‑1 (MQ=255)
                                                                                                                         tggATCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTccc                      <  7:195200/148‑1 (MQ=255)
                                                                                                                           gATCGCGTATCCATCTTAGGAGTACAGCATGACCAACTAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGt                    >  1:225312/1‑149 (MQ=255)
                                                                                                                            aTCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCt                  <  5:176852/149‑1 (MQ=255)
                                                                                                                             tCGCGTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTg                 <  2:315710/149‑1 (MQ=255)
                                                                                                                                 gTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCg                           >  7:267120/1‑135 (MQ=255)
                                                                                                                                 gTATCCATCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCg                           <  8:267120/135‑1 (MQ=255)
                                                                                                                                        tCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTGCa               <  5:193670/140‑1 (MQ=255)
                                                                                                                                        tCTTAGGAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTGCa               >  6:193670/1‑140 (MQ=255)
                                                                                                                                            aGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTGCAGGCGTTCGAAgg   <  4:232190/148‑1 (MQ=255)
                                                                                                                                             ggAGTACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTGCAGGCGTTCGAAGGc  <  8:213083/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCTTGGTAGGACCTTCGAACGACGCGCCGATCTTCTGGAAAGTCGCTCGCTTTCTGGAGAAGAGAAGAAGCTGCTCGAGGAATATCTCCGCGAGCGGGACGATAGTTAAACGTATCGATGGTGGATCGCGTATCCATCTTAGGAGCACAGCATGACCAACAAGATCATCCAGCAGCTCGAAGCCGAGCAGATGAGCAAGGAAATCCCGACCTTCGCACCAGGCGACACCGTTGTCGTTCAGGTTAAAGTGAAGGAAGGTGAGCGTTCCCGTCTGCAGGCGTTCGAAGGC  >  NC_002947/1670076‑1670364

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