Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 419,671 Δ1 bp coding (81/1320 nt) brnQ → branched chain amino acid transporter BrnQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913419,6700T.91.7% 81.6 / 0.8 24coding (80/1320 nt)brnQbranched chain amino acid transporter BrnQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (7/15);  minor base G (2/0);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGGCGGCGGTG  >  NC_000913/419539‑419798
                                                                                                                                   |                                                                                                                                
aTTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTcc                                                                                                                      >  1:685966/1‑143 (MQ=255)
    gCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATg                                                                                                                  >  1:505952/1‑143 (MQ=255)
    gCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACACACTGTCCCTCCAATg                                                                                                                  >  2:248065/1‑143 (MQ=255)
          aaaTTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGGAACAGCTTTTCCCTCCAGTGGCCGGc                                                                                                            >  2:485662/1‑143 (MQ=255)
                      acaccaCCATAGCCGCAGGCAGTTCTGTGTATTACCCATCAATTAAGATCGCGCGATATCATCGCTCTTGGCTTTACGCCATTCTCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACaggcaggc                                                                                                <  2:123158/143‑1 (MQ=255)
                        aacaCCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGCTATCATCGCTTTGGGCTTTATGACCCTCGCGTTGTTCGCCGGCGCAGGTAACGAATTTTGCCTCCAATGGTCCGCTTCCAGGCCGGCGa                                                                                              >  2:230373/1‑143 (MQ=255)
                                aTCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCt                                                                                      <  2:773742/143‑1 (MQ=255)
                                   cacaGGCAGTATGATTTATCACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGa                                                                                   <  2:149076/143‑1 (MQ=255)
                                    acaGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACCGGCCGGCGAAACCGTCTggtc                                                                                  >  1:284625/1‑141 (MQ=255)
                                             aTGATTTATGACCCGTCAATTCAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCATGCGAACACGTCTGGACTGCGGCAtt                                                                         <  2:322319/143‑1 (MQ=255)
                                                           aTCAATTAAGTTCGCGCGCTATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCATTGGGCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTa                                                           <  1:179846/143‑1 (MQ=255)
                                                                    gATCGCGCGCTATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTg                                                  <  1:159474/143‑1 (MQ=255)
                                                                           cgctatCACCGCTCTGGTCTTTAAGACCCTTTCTTGGTTCATCGGCGCAGGTAACAGATTTTCCCTCCATTCGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGcc                                           <  2:585635/140‑1 (MQ=255)
                                                                               atcagcGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTa                                       <  1:563816/143‑1 (MQ=255)
                                                                                  atcGCTCTGGGCTTTATGACATTTGAGTTGTTCGGTATTGCAGTTTTCAGTTTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTa                                    >  1:537069/1‑143 (MQ=255)
                                                                                              tttATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTCTTCACGGTGGTGgcg                        >  1:102301/1‑143 (MQ=255)
                                                                                                 aTGCCATTTGCGTTGGTCGTCGGCGCAGGTAACATATTCTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATCACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTg                     <  1:611651/143‑1 (MQ=255)
                                                                                                       tttGCGTTGTTCGTCGGCGCAGGTAACATATATTCCGTTCAAGCGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCaaaa               <  1:606604/143‑1 (MQ=255)
                                                                                                       gttGCGTTGTTCGTCGGCGCAGGTAACATACTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGTCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCaaaa               <  1:337737/142‑1 (MQ=255)
                                                                                                          gCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCCGCTTACAGGCAGGCGAACACGTCTGGGCTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGtt            <  1:57770/143‑1 (MQ=255)
                                                                                                            gttgttCGGCCGCGCAGGTAACAGATGGCCCCTCCGAGGGTCGGCTTACAGGCAGGCGAACACGGCTGGACTGCGGCATTCGGCTTCGTCATTGCTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTgg          <  2:439299/143‑1 (MQ=255)
                                                                                                             ttgttCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTggc         <  1:612898/143‑1 (MQ=255)
                                                                                                                 tcgtcgGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGgcggcg     <  1:159412/143‑1 (MQ=255)
                                                                                                                    tcgGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGGTGGCCTGCCGGTATTAACGGTAGTGGCGCGGGCAAAAGTTGGCGGCGgtg  >  1:606865/1‑143 (MQ=255)
                                                                                                                                   |                                                                                                                                
ATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGGCGGCGGTG  >  NC_000913/419539‑419798

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