Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 22,773 A→T 13.9% E128V (GAA→GTA)  ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091322,7730AT13.9% 125.9 / 12.3 43E128V (GAA→GTA) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base T (4/2);  total (25/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 = 9.92e-01

TGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGT  >  NC_000913/22626‑22918
                                                                                                                                                   |                                                                                                                                                 
tGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGaa                                                                                                                                                  >  1:310400/1‑149 (MQ=255)
tGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGa                                                                                                                                                   <  3:301972/148‑1 (MQ=255)
 gAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAAt                                                                                                                                                 >  1:241644/1‑149 (MQ=255)
     tCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATAc                                                                                                                                               >  2:338838/1‑147 (MQ=255)
     tCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATAc                                                                                                                                               <  1:338838/147‑1 (MQ=255)
       cAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcgg                                                                                                                                           >  2:391249/1‑149 (MQ=255)
          aGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggcac                                                                                                                                        <  1:391249/149‑3 (MQ=255)
              cTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTTTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCa                                                                                                                                    >  2:374764/1‑149 (MQ=255)
               tttCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAg                                                                                                                                   >  2:200288/1‑149 (MQ=255)
                     gTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGAc                                                                                                                             <  2:403497/149‑1 (MQ=255)
                      ttATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGACTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGAc                                                                                                                             >  4:234588/1‑148 (MQ=255)
                             tCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACg                                                                                                                     >  1:81285/1‑149 (MQ=255)
                                          ctggctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATc                                                                                                        >  2:284670/1‑149 (MQ=255)
                                                gggACTGCCACGGTTTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTg                                                                                                  <  1:374764/149‑1 (MQ=255)
                                                 ggACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTgg                                                                                                 <  2:81285/149‑1 (MQ=255)
                                                   aCTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGc                                                                                               >  1:229898/1‑149 (MQ=255)
                                                                 ccGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACt                                                                                  >  4:236611/1‑148 (MQ=255)
                                                                     tCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTc                                                                              >  4:141996/1‑148 (MQ=255)
                                                                            gAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggc                                                                                                                                          >  3:84920/1‑81 (MQ=255)
                                                                            gAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggc                                                                                                                                          <  4:84920/81‑1 (MQ=255)
                                                                            gAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGt                                                                      >  2:22612/1‑149 (MQ=255)
                                                                                gTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTAcc                                                                   >  4:171085/1‑148 (MQ=255)
                                                                                 tCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCt                                                                  >  3:58905/1‑148 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                               <  2:353962/100‑1 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                               >  1:353962/1‑100 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGAcc                                                              <  2:310400/149‑1 (MQ=255)
                                                                                      cAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCa                                                             >  4:268752/1‑148 (MQ=255)
                                                                                         gAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATgg                                                          <  3:268752/148‑1 (MQ=255)
                                                                                         gAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATgg                                                          >  4:256056/1‑148 (MQ=255)
                                                                                         gAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGa                                                         <  2:374236/149‑1 (MQ=255)
                                                                                         gAATACGGTAAGCCGGATGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATgg                                                          <  3:115419/148‑1 (MQ=255)
                                                                                                  aaGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCaaaa                                                 <  3:136786/148‑1 (MQ=255)
                                                                                                    gCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTg                                              >  1:356996/1‑149 (MQ=255)
                                                                                                                     ccgccgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCcg                              <  3:234588/148‑1 (MQ=255)
                                                                                                                            cgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTg                       <  3:236611/148‑1 (MQ=255)
                                                                                                                                 tCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcatcat                                 <  3:358978/133‑1 (MQ=255)
                                                                                                                                 tCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcatcat                                 >  4:358978/1‑133 (MQ=255)
                                                                                                                                     cgcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatc              >  4:402532/1‑148 (MQ=255)
                                                                                                                                      gcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat            <  1:29429/149‑1 (MQ=255)
                                                                                                                                      gcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat            >  2:190107/1‑149 (MQ=255)
                                                                                                                                      gcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatca             <  3:403216/148‑1 (MQ=255)
                                                                                                                                       ccAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat            <  3:402532/148‑1 (MQ=255)
                                                                                                                                            cgCCGCGTATTCGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCa       >  4:138622/2‑148 (MQ=255)
                                                                                                                                             gCCGCGTATTCGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCaa      >  3:359247/1‑148 (MQ=255)
                                                                                                                                                 cgAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCATAATCATCGGCAACGGt  >  4:241928/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGT  >  NC_000913/22626‑22918

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