Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 22,773 A→T 10.3% 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,7730AT10.3% 157.6 / 6.9 48E128V (GAA→GTA) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base T (2/3);  total (24/24)
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

AAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGG  >  NC_000913/22635‑22917
                                                                                                                                          |                                                                                                                                                
aaGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggc                                                                                                                                         <  1:434635/148‑1 (MQ=255)
              tATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATAcgcg                                                                                                                                           >  4:615855/1‑132 (MQ=255)
              tATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATAcgcg                                                                                                                                           <  3:615855/132‑1 (MQ=255)
                 gACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATa                                                                                                                                               <  4:128636/125‑1 (MQ=255)
                 gACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATa                                                                                                                                               <  1:440338/125‑1 (MQ=255)
                 gACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATa                                                                                                                                               >  3:128636/1‑125 (MQ=255)
                 gACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATa                                                                                                                                               >  2:440338/1‑125 (MQ=255)
                  aCTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGaa                                                                                                                                                 >  3:450532/1‑122 (MQ=255)
                  aCTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGaa                                                                                                                                                 <  4:450532/122‑1 (MQ=255)
                   cTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAccc                                                                                                                                    <  2:60467/134‑1 (MQ=255)
                   cTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAccc                                                                                                                                    >  1:60467/1‑134 (MQ=255)
                    tCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTcgcg                                                                                                                                           >  4:178426/1‑126 (MQ=255)
                    tCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTcgcg                                                                                                                                           <  3:178426/126‑1 (MQ=255)
                    tCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggc                                                                                                                                         >  1:771615/1‑128 (MQ=255)
                    tCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggc                                                                                                                                         <  2:771615/128‑1 (MQ=255)
                           aTGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                              >  2:464153/1‑148 (MQ=255)
                                  tggctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATc                                                                                                       >  4:510678/1‑148 (MQ=255)
                                      tggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTc                                                                                                   <  1:464153/148‑1 (MQ=255)
                                        ggACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATACACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGa                                                                                                                             <  3:492256/120‑1 (MQ=255)
                                        ggACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATACACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGa                                                                                                                             >  4:492256/1‑120 (MQ=255)
                                               cACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCaa                                                                                                                 >  1:1123307/1‑125 (MQ=255)
                                               cACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCaa                                                                                                                 <  2:1123307/125‑1 (MQ=255)
                                                aCGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCt                                                                                         >  1:722574/1‑148 (MQ=255)
                                                aCGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCt                                                                                         <  3:510678/148‑1 (MQ=255)
                                                        ccGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTc                                                                                                   <  3:810196/130‑1 (MQ=255)
                                                        ccGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTc                                                                                                   >  4:810196/1‑130 (MQ=255)
                                                         cGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTg                                                                                <  2:722574/148‑1 (MQ=255)
                                                           aTCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGa                                                                                                                                                  >  3:106302/1‑80 (MQ=255)
                                                           aTCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGa                                                                                                                                                  <  4:106302/80‑1 (MQ=255)
                                                                gCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTggg                                                                                      >  3:36215/1‑135 (MQ=255)
                                                                gCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTggg                                                                                      <  4:36215/135‑1 (MQ=255)
                                                                    aaaGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                              <  3:704709/107‑1 (MQ=255)
                                                                    aaaGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                              >  4:704709/1‑107 (MQ=255)
                                                                        tCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCa                                                                                                                                   >  3:437613/1‑82 (MQ=255)
                                                                        tCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCa                                                                                                                                   <  4:437613/82‑1 (MQ=255)
                                                                               aGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATg                                                          >  3:9237/1‑148 (MQ=255)
                                                                                       gTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGAc                                                                                                                            >  2:716537/1‑74 (MQ=255)
                                                                                       gTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGAc                                                                                                                            <  1:716537/74‑1 (MQ=255)
                                                                                             cGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGa                                            <  4:9237/148‑1 (MQ=255)
                                                                                                        ttCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGAc                                                                                  <  1:241010/99‑1 (MQ=255)
                                                                                                        ttCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGAc                                                                                  >  2:241010/1‑99 (MQ=255)
                                                                                                               ccgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCAt                                                           >  2:81445/1‑115 (MQ=255)
                                                                                                               ccgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTAGTCGCACCCGTACCTGACCAt                                                           <  1:81445/115‑1 (MQ=255)
                                                                                                                          cgcgCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTa                                                                    <  4:324464/95‑1 (MQ=255)
                                                                                                                          cgcgCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTa                                                                    >  3:324464/1‑95 (MQ=255)
                                                                                                                           gcgcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCaaa                                                 <  4:1102102/113‑1 (MQ=255)
                                                                                                                           gcgcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCaaa                                                 >  3:1102102/1‑113 (MQ=255)
                                                                                                                            cgcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACt                                                      >  4:1131156/1‑107 (MQ=255)
                                                                                                                            cgcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACt                                                      <  3:1131156/107‑1 (MQ=255)
                                                                                                                                      cgcgAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTg                                                                                        <  4:308458/63‑1 (MQ=255)
                                                                                                                                      cgcgAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTg                                                                                        >  3:308458/1‑63 (MQ=255)
                                                                                                                                       gcgAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACgg  >  2:675196/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                
AAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGG  >  NC_000913/22635‑22917

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