Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 555,549 T→C R313R (CGT→CGC cysS → cysteinyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913555,5490TC95.6% 144.6 / ‑4.6 45R313R (CGT→CGCcysScysteinyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (23/20);  total (24/21)
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

GAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGTGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTATTGAAGC  >  NC_000913/555453‑555644
                                                                                                |                                                                                               
gAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGcgcgcgcg                                                                                              >  1:578512/1‑100 (MQ=255)
   aTACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCgcggcg                                                                                           >  1:27731/1‑100 (MQ=255)
     actacGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCt                                                                                         >  1:629897/1‑100 (MQ=255)
      ctacGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTg                                                                                        <  1:1156810/100‑1 (MQ=255)
       tacGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTgg                                                                                       <  2:1057599/100‑1 (MQ=255)
           acgCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCg                                                                                   <  2:574566/100‑1 (MQ=255)
                 aaaCCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTAAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTa                                                                             <  2:300824/100‑1 (MQ=255)
                       tGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACActgc                                                                       <  1:522183/100‑1 (MQ=255)
                                tCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTTTTGTTGGCGGGTGCGGCGCCCGAGGAGATCTACACTGCGCTGCGCgg                                                              <  2:749690/100‑1 (MQ=255)
                                 ccTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGc                                                             >  1:997994/1‑100 (MQ=255)
                                 ccTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGc                                                             >  2:176764/1‑100 (MQ=255)
                                 ccTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGc                                                             >  1:514542/1‑100 (MQ=255)
                                 ccTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGc                                                             <  1:433299/100‑1 (MQ=255)
                                  cTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGca                                                            >  1:43662/1‑100 (MQ=255)
                                  cTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGca                                                            >  1:432003/1‑100 (MQ=255)
                                  cTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGca                                                            >  1:249975/1‑100 (MQ=255)
                                    gATGGCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGcaca                                                          <  1:517509/100‑1 (MQ=255)
                                        tCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATa                                                      <  2:1129647/100‑1 (MQ=255)
                                         cGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATaa                                                     >  1:862384/1‑100 (MQ=255)
                                            gCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAAc                                                  <  2:608343/100‑1 (MQ=255)
                                             ccACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAAcc                                                 >  1:533688/1‑100 (MQ=255)
                                             ccACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAAcc                                                 >  1:26507/1‑100 (MQ=255)
                                              cACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCg                                                >  2:1294768/1‑100 (MQ=255)
                                                 tATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTg                                             <  2:1251295/100‑1 (MQ=255)
                                                  aTCGCAGCCAGTTGGACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTgc                                            >  1:34007/1‑100 (MQ=255)
                                                    cGCAGCCAGTTGAACTACAGAGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTgcgc                                          >  2:55660/1‑100 (MQ=255)
                                                     gCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGcc                                         >  2:265309/1‑100 (MQ=255)
                                                          cAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACATTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCg                                    <  2:801401/100‑1 (MQ=255)
                                                          cAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCg                                    >  2:361350/1‑100 (MQ=255)
                                                          cAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCg                                    >  2:146902/1‑100 (MQ=255)
                                                          cAGTTGAACGACAGCGAAGAGATCCTGAAGCGGGCGGGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGcc                                     <  2:85666/99‑1 (MQ=255)
                                                             ttGAACTACAGCGAAGAGAACCTGAAGAAGGCGCGCGCGGCGCTGAAGCGTCTCTACACAGCACTCCGCAGCACAGATAAAACCGTTGCGcccac                                      >  2:12655/1‑92 (MQ=255)
                                                                aaCTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCg                              <  1:1003891/100‑1 (MQ=255)
                                                                 aCTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGa                             <  1:710905/100‑1 (MQ=255)
                                                                  cTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGaa                            >  2:748406/1‑100 (MQ=255)
                                                                  cTACAGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGaa                            >  1:1164444/1‑100 (MQ=255)
                                                                      aGCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCg                         >  1:153406/1‑99 (MQ=255)
                                                                       gCGAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGtt                       >  1:1006441/1‑100 (MQ=255)
                                                                         gAAGAGAACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTg                     <  1:1059205/100‑1 (MQ=255)
                                                                          aaGAGAACCTGAAGCAGGCGCGTGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGa                    >  2:638208/1‑100 (MQ=255)
                                                                            gagaACCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAg                  <  1:265312/100‑1 (MQ=255)
                                                                                aCCTGAAGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCGGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAgcgcg              <  2:487574/100‑1 (MQ=255)
                                                                                      aGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTAt        <  1:830363/100‑1 (MQ=255)
                                                                                      aGCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTAt        >  1:1068022/1‑100 (MQ=255)
                                                                                       gCAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTAtt       <  2:1041709/100‑1 (MQ=255)
                                                                                        cAGGCGCGCGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTATTg      <  2:508867/100‑1 (MQ=255)
                                                                                           gcgcgcgcGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTATTGAAg   <  2:249972/100‑1 (MQ=255)
                                                                                            cgcgcgcgGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTATTGAAGc  >  1:1167398/1‑100 (MQ=255)
                                                                                            cgcgcgcgGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTATTGAAGc  <  2:997985/100‑1 (MQ=255)
                                                                                                |                                                                                               
GAAATACTACGACGCGGAAACCGTGCGTTACTTCCTGATGTCGGGCCACTATCGCAGCCAGTTGAACTACAGCGAAGAGAACCTGAAGCAGGCGCGTGCGGCGCTGGAGCGTCTCTACACTGCGCTGCGCGGCACAGATAAAACCGTTGCGCCTGCCGGTGGCGAAGCGTTTGAAGCGCGCTTTATTGAAGC  >  NC_000913/555453‑555644

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