Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,818,044 A→G 16.5% intergenic (‑186/+13) argQ ← / ← argZ tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,818,0440AG16.5% 89.0 / 11.0 42intergenic (‑186/+13)argQ/argZtRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base G (3/4);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAG  >  NC_000913/2817899‑2818192
                                                                                                                                                 |                                                                                                                                                    
ctaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAttt                                                                                                                                                   <  1:314134/149‑1 (MQ=255)
        atatCAATTAAGGGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAat                                                                                                                                           <  1:555080/149‑1 (MQ=255)
           tCAATTAAGAGATGGTGCATCCGGGGGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                        <  1:123925/149‑1 (MQ=255)
           tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGGTCGCCCTTAGGGCCGTTGCGGGAGAAACGGTATCATCCCGAGTGGCTGCGGTTAGTTCCGCTACCTTTTCGATACCTTCATTGCTGAAGATTTAGGCGGATATg                                                                                                                                        >  2:627647/1‑149 (MQ=17)
            cAATTACGAGGTGTCGAATCCGGGAGTATTACCCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg                                                                                                                                       <  1:84313/149‑1 (MQ=18)
              aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCTTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAAATTTCGGAGAATATGGTg                                                                                                                                     >  2:395212/1‑149 (MQ=255)
              aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTAGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTg                                                                                                                                     >  1:476434/1‑149 (MQ=255)
                tAAGAGATGGTGCATCCGGGAGGATTACTCGCCTCGCGCCCCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTGGTGCCGGCGGTTAATGCCGCTACCGCCCCGATACCTTCATTGCAGAAGAGTTCGGAGAACATgggcca                                                                                                                                   >  2:462535/1‑145 (MQ=14)
                        ggTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggag                                                                                                                           >  2:664494/1‑149 (MQ=255)
                             aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTc                                                                                                                      <  1:607493/149‑1 (MQ=255)
                                cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTGCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGAATTCGaa                                                                                                                   >  2:311441/1‑149 (MQ=255)
                                       aTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                            >  2:361728/1‑149 (MQ=255)
                                            tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGGGAGAATGGTGCATCCGGGAGGATTCGAACCCCCGACCGCt                                                                                                       >  1:548590/1‑149 (MQ=255)
                                               ccTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTGGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCgg                                                                                                    <  1:166142/149‑1 (MQ=255)
                                                cTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGt                                                                                                   <  2:580777/149‑1 (MQ=255)
                                                   gCGGCTCCCCTTTCGGGCCGTTGCTAAAGCAACGTTCTCACCCATTCTGCCCGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCg                                                                                                <  1:240076/149‑1 (MQ=17)
                                                      gCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                                             >  2:559531/1‑149 (MQ=255)
                                                                     cGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTCCTCTATcc                                                                              >  1:260297/1‑149 (MQ=255)
                                                                       tttcTAAAGCAACGTTATCATCCCTAGTGCTTGCGACATATTCCGTTGCCGCTGCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCCGGAGGATTCGAACCTCCGACTGCTCGGTTCGTAGCCGAGTACTCTATcccg                                                                            <  2:515084/146‑3 (MQ=17)
                                                                       ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCCCGGTTCGTAGCCGAGTACTCTATCCAg                                                                            >  2:321331/1‑149 (MQ=255)
                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                        >  2:263701/1‑149 (MQ=255)
                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGGTCGTAGCCGAGTACTCTATCCACCtac                                                                        >  2:623500/1‑147 (MQ=255)
                                                                                  aCGTTATGTTACCGAGTGCGTCCGGTTAATTCCGCTACCTCTTCCATATCTTCATTGCTGAAGATTTCGGAGAATAGGGTGCACCCGGGAGGCTTCGAACCTCCGTCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACgg                                                                 <  2:184366/149‑1 (MQ=14)
                                                                                  aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGGACCCCCTACCGCCCTTTTCGAACCCGAATACTCTATCCTGCTGAGcgccgg                                                                 >  1:698454/1‑144 (MQ=17)
                                                                                   cGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTGTATCCAGCTGGGCTACgcc                                                                >  2:348109/1‑147 (MQ=255)
                                                                                           tcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                        >  1:590747/1‑149 (MQ=255)
                                                                                                    ccTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttat                                               >  1:728135/1‑149 (MQ=255)
                                                                                                     cTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGGACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttaat                                              >  2:108279/1‑147 (MQ=255)
                                                                                                                     tACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAGCGGGAAATTTATTTTACTGCTGATATTGa                              <  1:446432/149‑1 (MQ=255)
                                                                                                                        ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGAATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAc                           <  1:321331/149‑1 (MQ=255)
                                                                                                                         tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTCATTTTACTGCTGATATTAACAcc                          >  1:310562/1‑149 (MQ=255)
                                                                                                                             cGATACCTTCATTGCTGACGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                      <  1:32243/149‑1 (MQ=255)
                                                                                                                                 aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                           <  1:265447/124‑1 (MQ=2)
                                                                                                                                 aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                           >  2:265447/1‑124 (MQ=2)
                                                                                                                                  ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                 <  1:263701/149‑1 (MQ=255)
                                                                                                                                   cTTCATTGCTGAATGTTTCGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg                <  1:593282/149‑1 (MQ=255)
                                                                                                                                   cTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAGCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg                >  1:115315/1‑149 (MQ=255)
                                                                                                                                     tCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     <  1:164211/126‑1 (MQ=255)
                                                                                                                                     tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     <  1:164299/126‑1 (MQ=255)
                                                                                                                                     tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     >  2:164299/1‑126 (MQ=255)
                                                                                                                                     tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     >  2:164211/1‑126 (MQ=16)
                                                                                                                                          gCTGAAGATGTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCTCGGTTCGTAGCCGAGTACTCGATCCAGCTGAGCTACGGTTGCATCTGGAAATTTGTTTTTCTGCTGATTTTGATACCGCTACTAAAGCCATATc         <  2:380894/149‑1 (MQ=14)
                                                                                                                                                 aTTGCGGAGTATATGGTGCATCCGGGAGGCTTCGGACCTCCGTCCGCTCGGTTTGTAGCCGACTACCCAAACCTGCTGGCCTACGGACGCAGTGGGAAATTTATTTGACTGCCGCAATTGATACCGCTACTAAAGCCATATCAAGTAag  <  2:229080/149‑1 (MQ=9)
                                                                                                                                                 |                                                                                                                                                    
CTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAG  >  NC_000913/2817899‑2818192

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