Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 419,670 Δ1 bp 81.1% coding (80/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.81.1% 90.9 / 24.5 37coding (80/1320 nt)brnQbranched chain amino acid transporter BrnQ
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base . (17/13);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGGCGGCGGTGTTGACAGCCTCAG  >  NC_000913/419533‑419811
                                                                                                                                         |                                                                                                                                              
ccAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCCTCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATAtttt                                                                                                                                          >  1:421748/1‑143 (MQ=255)
 cAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTCTGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTCCGTTCTTCGTCGGCGCAGGTAACATATTTTc                                                                                                                                         >  1:592076/1‑143 (MQ=255)
  tgcgATGGGCACAAAAATTCAACGTTACATCCCCACATCCACGGGCAGTATGATTTATGACCCATCAATGAAGATCGCGCGATATCATCGCTCTGGGCTTTATGTCATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTcc                                                                                                                                        <  1:519378/140‑1 (MQ=255)
    ggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCt                                                                                                                                      <  2:197782/143‑1 (MQ=255)
       ttGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGAGCAGGTAACATATTTTCCCTCCa                                                                                                                                   <  2:969758/143‑1 (MQ=255)
             caAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGt                                                                                                                              <  2:771397/143‑1 (MQ=255)
             caAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTc                                                                                                                             <  2:288023/143‑1 (MQ=255)
                aaaTTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCCATGGTCGGc                                                                                                                          >  1:494162/1‑143 (MQ=255)
                  aTTTAACGTTACAACACCACATCCCCAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATTACATTTGCGTTGTTCATCGGCCCAGGTAACATATTTTCCCTCAGATCGTCAGCtt                                                                                                                       >  2:790207/1‑143 (MQ=255)
                    ttAGCGTTACAACTTGACATCCACAGGCAGTATGATTTAGGACCCATCAATTAAGCTCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTAc                                                                                                                     <  2:890820/143‑1 (MQ=255)
                        cGTTACAACACCAAATCCGCAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTAcagg                                                                                                                  <  2:43378/143‑1 (MQ=255)
                         gTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTAcaggca                                                                                                                <  2:130986/143‑1 (MQ=255)
                           tACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTAcaggcagg                                                                                                              >  1:458628/1‑143 (MQ=255)
                           tACAACACCACATCCACAGGCAGTATGATATATGACCCATCAATTAAGATCAAAAGTGAGCATCGCCCTGTTCATTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTTAcaggca                                                                                                                >  2:1077940/1‑143 (MQ=255)
                                caccacATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTAcaggcag                                                                                                               <  1:695602/137‑1 (MQ=255)
                                caccacATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTAcaggcag                                                                                                               >  2:695602/1‑137 (MQ=255)
                                caccacATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGCac                                                                                                         >  2:987289/1‑143 (MQ=255)
                                  ccacATGCACAGGCAGTATGATTTATGCCACATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAacac                                                                                                       <  2:592076/143‑1 (MQ=255)
                                   cacaTCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAAGATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGCACACg                                                                                                      >  1:278410/1‑143 (MQ=255)
                                     caTCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCCGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCACGCGAACCCGTc                                                                                                    >  2:319996/1‑143 (MQ=255)
                                     caTCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGACTTTATGACAGTTGCGTTGTTCGGCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTc                                                                                                    >  1:697180/1‑143 (MQ=255)
                                       tCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTg                                                                                                  <  1:679037/143‑1 (MQ=255)
                                                                cATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGGTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGAGCCGCGTACAGGCAGGCGACCACGTCTGGACTCCGGCATGCGGCATcgtgac                                                                          >  2:830650/1‑138 (MQ=255)
                                                                    aaTTAAGATCGTGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTACTCTCCCCTGGTCCGCTCGGTGGCGGGCGCACACGTCTGGACTGCCGCCTCCGCCTTCCTCATTACTg                                                                     >  2:783952/1‑143 (MQ=255)
                                                                     aTTAAGATCGCGCGATATCATCGCTCTGGGCTTGATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCGCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGc                                                                    <  1:971583/143‑1 (MQ=255)
                                                                       tAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGTAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTCCCg                                                                  >  1:142641/1‑143 (MQ=255)
                                                                           aTCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTgg                                                              >  2:114904/1‑143 (MQ=255)
                                                                            tCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGGTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCCTTCGCCTCCCTCATTACTGCCGTTGGc                                                             >  2:873628/1‑143 (MQ=255)
                                                                              gcgcgATATCATCGCTCTGGGCCTTATGACATTCGCGTTGTTCGTCGGCGCAGGTAACATATCTTCCCTCCAATGGCCGGCTTACAGGCCGGCGAACACTCCTGGACTGCGTCATTCGGCTCGCTCGTCACCTCCGTCGGcca                                                           >  2:268871/1‑142 (MQ=255)
                                                                                                      tctGACATTTGCGTTGTTCGTCGGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCt                                   <  2:171504/141‑1 (MQ=255)
                                                                                                        tGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCt                                   <  1:1077940/143‑1 (MQ=255)
                                                                                                                    tgttCGTCGGCGCAGGTAACATATCTTCCCTCCAACGGTCGGCTTACAGGCCGGCGAACACGTCTGGACTGCGGCATCCGCATTCCTCCTTACTGCCGTTGGCCTGCCGGCATTAACGGCAGTGGCGCTGGCAAAAGTTggcg                     <  1:496609/143‑1 (MQ=255)
                                                                                                                      tttgtcgtcTCAGGTATCATCACTTTCGCTCCAAGTGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGTTGGCAAAAGTTggcgg                    <  2:462331/140‑1 (MQ=255)
                                                                                                                        cgtctGCGCAGGTAACATTATTTGCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGgcggcg                  <  2:740373/143‑1 (MQ=255)
                                                                                                                        cgtcgGCGCAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGGcggcgg                 <  2:697180/143‑1 (MQ=255)
                                                                                                                              cgcAGGTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCACAGGTTGGCGGCGGTGTTGa           >  2:33012/1‑143 (MQ=255)
                                                                                                                                   gTAACATATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGAGGCATTCGGCTGCCTCATTACTGCCGTTGGCCTGCCGGTGTTAGCGCTAGGGGAGCTGCGAAGAGTTGGCGGGGGTGTTGACAGcc      >  2:33454/1‑143 (MQ=255)
                                                                                                                                     tagcttattGTCCCTCCAGTGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCAGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGGCGGCGGTGTTGACAGCCt     <  2:878987/139‑1 (MQ=255)
                                                                                                                                        attattTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAGGTTGGCGGCGGTGTTGACAGCCTCAg  >  1:912047/1‑143 (MQ=255)
                                                                                                                                         ttattTTCCCTCCAATGGTCGGCTTTACAGGCAGGCGAACACGTCTGGACCGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCCAAAGGTGGCGGCGGTGTGGACAGCCCCAg  >  2:257599/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACATTATTTTCCCTCCAATGGTCGGCTTACAGGCAGGCGAACACGTCTGGACTGCGGCATTCGGCTTCCTCATTACTGCCGTTGGCCTGCCGGTATTAACGGTAGTGGCGCTGGCAAAAGTTGGCGGCGGTGTTGACAGCCTCAG  >  NC_000913/419533‑419811

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: