Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 420,946 T→A 16.1% intergenic (+36/‑40) brnQ → / → proY branched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val/proline‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913420,9460TA16.1% 84.0 / 10.5 31intergenic (+36/‑40)brnQ/proYbranched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val/proline‑specific permease
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base A (2/3);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

TACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTAC  >  NC_000913/420798‑421078
                                                                                                                                                    |                                                                                                                                    
tACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGAt                                                                                                                                      <  3:83901/149‑1 (MQ=255)
    gCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGccc                                                                                                                                  <  4:304660/149‑1 (MQ=255)
        gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTg                                                                                                                              <  1:464796/149‑1 (MQ=255)
          cGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTGGt                                                                                                                            <  3:156341/149‑1 (MQ=255)
           gAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTGGtt                                                                                                                           <  1:323285/149‑1 (MQ=255)
           gAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTGGtt                                                                                                                           >  2:323285/1‑149 (MQ=255)
                   tCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtg                                                                                                                   >  4:289383/1‑149 (MQ=255)
                   tCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATAGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtg                                                                                                                   >  1:544196/1‑149 (MQ=255)
                                     aaCAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAAt                                                                                                 <  2:463929/149‑1 (MQ=255)
                                          tggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATg                                                                                            <  2:544196/149‑1 (MQ=255)
                                            gtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGa                                                                                          >  4:318982/1‑149 (MQ=255)
                                                         tGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAg                                                                             <  2:463396/149‑1 (MQ=255)
                                                             cATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAg                                                                                                     <  4:105864/121‑1 (MQ=255)
                                                             cATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAg                                                                                                     >  3:105864/1‑121 (MQ=255)
                                                               ttATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAg                                                                       >  2:36016/1‑149 (MQ=255)
                                                                 aTCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCg                                                                     >  1:369952/1‑149 (MQ=255)
                                                                            cgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTaa                                                                                    <  2:223475/123‑1 (MQ=255)
                                                                            cgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTaa                                                                                    >  1:223475/1‑123 (MQ=255)
                                                                                  aGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGcc                                                    >  1:263854/1‑149 (MQ=255)
                                                                                              gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTTACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGAt                                                                                                            >  4:511570/1‑81 (MQ=255)
                                                                                              gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGAt                                                                                                            <  3:511570/81‑1 (MQ=255)
                                                                                              gACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                        <  4:47704/149‑1 (MQ=255)
                                                                                                 cTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATgg                                     <  3:318982/149‑1 (MQ=255)
                                                                                                       cgcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTgg                               <  2:495449/149‑1 (MQ=255)
                                                                                                        gcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTggg                              <  3:289383/149‑1 (MQ=255)
                                                                                                                aaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCaa                      <  2:163806/149‑1 (MQ=255)
                                                                                                                    aCTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTgg                  <  1:36016/149‑1 (MQ=255)
                                                                                                                         aCATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGaa                                                                                 >  3:485019/1‑81 (MQ=255)
                                                                                                                         aCATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGaa                                                                                 <  4:485019/81‑1 (MQ=255)
                                                                                                                                    aaCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTAc  <  2:187195/149‑1 (MQ=255)
                                                                                                                                    aaCCACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaaa                                                                                        <  3:420091/63‑1 (MQ=21)
                                                                                                                                    aaCCACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaaa                                                                                        >  4:420091/1‑63 (MQ=21)
                                                                                                                                                    |                                                                                                                                    
TACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTAC  >  NC_000913/420798‑421078

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