Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 2,817,877 C→T 14.4% intergenic (‑17/+182) argQ ← / ← argZ tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi2,817,8770CT14.4% 106.0 / 3.6 51intergenic (‑17/+182)argQ/argZtRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base T (5/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTC  >  NC_000913_3_bme_pgi/2817737‑2818023
                                                                                                                                            |                                                                                                                                                  
tACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGG‑AAATTTATTTTACTGCTGa                                                                                                                                           <  2:485695/149‑1 (MQ=9)
     ttttcttTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATAtt                                                                                                                                       <  1:245882/146‑1 (MQ=255)
       gttttttCGATACCTTCGTTGCTGAAGTTGACGGAGAATATGGTGTATACTGGAGGATTCGAACCTGTTACCGCTGGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATCTTATTCTACTGCTGATATTGa                                                                                                                                     <  2:150802/144‑1 (MQ=11)
       aTCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGa                                                                                                                                     <  1:860498/149‑1 (MQ=255)
          tctTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTAATTTTACTGCTGATTTTGTTAc                                                                                                                                  >  1:80327/1‑149 (MQ=11)
          tctTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGGACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACCTTATTCTCCTTCTATTATTTATAc                                                                                                                                  >  2:158554/1‑149 (MQ=14)
           ctccGATACCTTCATTCCTGAAGGTTACGGAGACTATGGTGCCTCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTATTCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATAcc                                                                                                                                 <  2:49312/146‑1 (MQ=14)
             tCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGc                                                                                                                               <  1:178665/149‑1 (MQ=255)
               tattccTTAATTGCTGAGGTTTACGGAGCATATGGTGCATCCGGGAGGATGCGATCCTCCGACCTCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGTGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGcta                                                                                                                             <  2:46197/145‑1 (MQ=14)
               gATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTAGCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGcta                                                                                                                             <  1:347958/149‑1 (MQ=255)
                atatCTTCATTGATGAAGGTTATTGATAATATTGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGctac                                                                                                                            <  1:216918/149‑1 (MQ=14)
                aTACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGTATGCATCGGGAAAACTTATTTTACTGATGATCTTTATAccaccac                                                                                                                            >  2:280506/1‑144 (MQ=14)
                aTACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTActgc                                                                                                                                              <  2:792097/131‑1 (MQ=255)
                aTACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTActac                                                                                                                                              >  1:792097/1‑129 (MQ=255)
                 tACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCCGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGctact                                                                                                                           >  2:270730/1‑149 (MQ=255)
                  aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGctacta                                                                                                                          <  1:831196/149‑1 (MQ=255)
                    ctttaatGCTGAAGGGTATGGACAATATGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCTCCACTCATAGCCGCGGACTCTGTCGAGCTGATCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTaaa                                                                                                                        <  2:143394/143‑1 (MQ=14)
                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGGACATCCGACCGCGCGGTGCGAAGCCGAGTACTCTATCCAACTGAGCGACGGATGCCTCGGGAAAACTAATTAAACGGCTGACATTGATACCCCAATTAAAg                                                                                                                       >  2:550026/1‑149 (MQ=14)
                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTTTACTGCTATTTTTGATACCCCTACCAAAg                                                                                                                       >  2:660357/1‑149 (MQ=14)
                     ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTaagg                                                                                                                       >  1:186879/1‑147 (MQ=255)
                          tGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCata                                                                                                                  <  1:672462/149‑1 (MQ=255)
                             tGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCa                                                                                                               >  1:566712/1‑149 (MQ=255)
                               aaGGTTACGGAGAATATGGTACATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAt                                                                                                             <  1:663260/149‑1 (MQ=255)
                                aGGTTACGGAGAATATGGTGCATCCGGGTGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAtt                                                                                                            <  2:271671/149‑1 (MQ=255)
                                    tACGGACAAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAaga                                                                                                        <  1:754176/149‑1 (MQ=255)
                                       ggAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGAt                                                                                                     >  1:554664/1‑149 (MQ=255)
                                                     aTCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCCTCGGGCATTCTTATTCTACTGCCGATATTGATCGCGCTACTAATGTCATTTTAATTAAGAGATGGTGCATCCGggag                                                                                       <  2:572335/149‑1 (MQ=11)
                                                      tCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTTTACTCCTGATATTGATACCGCTACAAAAGCCATATCAATTAAGAAATTCTGAATGCGggagg                                                                                      >  2:754086/1‑149 (MQ=11)
                                                        ctggaggaTGAGAACCTCTGAGCGCTCGGTTCGTAGCGGAGTACTCTATCCAGCGGAGCTACGGATGCATCGGGAGAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGTCATATCAATTAAGAGATGGTGCATCCGGGAGGAt                                                                                    <  2:887442/147‑1 (MQ=14)
                                                                 tCGAACCTCCGACCGCTCGGTTCGTAGCCGCGTACTCTATCCAGCTGAGCTACGGCTGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGCAGGATTACTCGCCt                                                                           >  1:285210/1‑149 (MQ=255)
                                                                   gaccttCCGACCGCTCTGTTCGTGGCTGAGTACTCTTTCCAGCTGTGCTACGGATGCATTGGGAGTACTTTTTTTTCTGTTGATATTGATACCGCTACTAAAGCTATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTg                                                                         <  1:323628/144‑1 (MQ=14)
                                                                   gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTATTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGGGGATTACTCGCCTTg                                                                         >  1:25112/1‑149 (MQ=255)
                                                                              cGCTCGGTTCGTAGCCGAGTACTGTATAGAGCTGAGCTGCGGATGCATCGGGGAAAGTTGTTGTTCTGGTGGTATTGATACTGCTACTAAAAACATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                              >  2:161859/1‑149 (MQ=17)
                                                                              cGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                              >  1:607007/1‑149 (MQ=255)
                                                                               gCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAAATGGTGCATCCGGGAGGATTAATCGCCTTGCGGCTCGCGCtt                                                             >  1:460027/1‑149 (MQ=255)
                                                                                  cGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                          <  1:700297/149‑1 (MQ=255)
                                                                                         tAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGG‑AAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTCCATCCGGCAGAATAACTCGCCTTCCGGCTCGCCCTTCGGGCCGttcc                                                  >  1:63587/1‑147 (MQ=11)
                                                                                                   ctctATCCTGCTGGTCTAGGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCCCTACTCAAACCATCTCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCTCTTCGGGCCGTTGCTAAAGCAAc                                         <  2:1623/149‑1 (MQ=17)
                                                                                                    tctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTGTAGTGCTGATATGGATGCAGCGATTAAGATCATGGCAGTTAAGAGATCGTGAATCCGGGAGGATTAGTAGCCTTGACACTTGGCTTTGGAGCCGATGCTTAAGCacag                                        >  2:384493/1‑146 (MQ=9)
                                                                                                          cAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatca                                  >  1:600919/1‑149 (MQ=255)
                                                                                                                 gCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCtcg                           >  1:874304/1‑147 (MQ=255)
                                                                                                                    aCGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGACCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGc                        >  1:43866/1‑149 (MQ=255)
                                                                                                                     cGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCCGGCCGTTTCTTAAGCAACGTTTTCCTCCCTTGTTcc                       >  2:741137/1‑149 (MQ=17)
                                                                                                                      ggATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCt                      >  1:242075/1‑149 (MQ=255)
                                                                                                                        aTGCTTCGGGAAAACTTATTGTACTGCTGTTATTGATACCGCTTCTATTGTTGGTTTGGTTTTGTGTTGGTGCTTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGc                    <  1:248379/149‑1 (MQ=14)
                                                                                                                        aTGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGc                    >  1:612214/1‑149 (MQ=255)
                                                                                                                             tCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTa               <  1:148214/149‑1 (MQ=255)
                                                                                                                              tgggAAAACTTATTCTACTGCTGATATTTATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa              <  1:144736/148‑1 (MQ=255)
                                                                                                                              cGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa              <  1:270730/149‑1 (MQ=255)
                                                                                                                                     aCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCt       >  1:428490/1‑149 (MQ=255)
                                                                                                                                          ttCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCACCGTTATCCTCCCTAGTGCCTGCGGTTATTTCCTCCACctc  >  2:429567/1‑149 (MQ=14)
                                                                                                                                            |                                                                                                                                                  
TACCGCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTC  >  NC_000913_3_bme_pgi/2817737‑2818023

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