Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,726,620 G→A 9.2% D254D (GAC→GAT)  brnQ ← branched‑chain amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,726,6200GA 9.2% 151.8 / 4.6 54D254D (GAC→GAT) brnQbranched‑chain amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base A (2/3);  total (25/29)
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.03e-01

ACTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGA  >  NC_002947/4726472‑4726765
                                                                                                                                                    |                                                                                                                                                 
aCTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATg                                                                                                                                                   >  7:147100/1‑149 (MQ=255)
    ggCGCAGGCACAGGTCAGGCCAACCGCGGTCACCATGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTc                                                                                                                                                 <  7:131790/147‑1 (MQ=255)
    ggCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGt                                                                                                                                               <  4:197491/149‑1 (MQ=255)
           gCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCTAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCACGACATCGCGGCtgac                                                                                                                                        <  2:200646/149‑3 (MQ=255)
            cacaGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCtggtg                                                                                                                                       <  5:16896/149‑1 (MQ=255)
            cacaGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCtggtg                                                                                                                                       <  5:79684/149‑1 (MQ=255)
            cacaGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCtggtg                                                                                                                                       <  6:4253/149‑1 (MQ=255)
                 gTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCAcc                                                                                                                                  >  5:304340/1‑149 (MQ=255)
                           ccGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGGACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAAc                                                                                                                        <  8:147100/149‑1 (MQ=255)
                                   aCCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGAt                                                                                                                >  2:136581/1‑149 (MQ=255)
                                     caggcaggCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATg                                                                                                               >  1:114714/1‑148 (MQ=255)
                                        gcaggcCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGa                                                                                                           >  8:113927/1‑149 (MQ=255)
                                               aGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGCACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGc                                                                                                    >  1:290811/1‑149 (MQ=255)
                                                gcgcGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGc                                                                                                    >  4:135498/1‑148 (MQ=255)
                                                     aTCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGccc                                                                                              <  6:304340/149‑1 (MQ=255)
                                                      tCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGccc                                                                                              >  1:89279/1‑148 (MQ=255)
                                                         gCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGaca                                                                                          <  7:113927/149‑1 (MQ=255)
                                                                 aGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGAt                                                                                  >  5:133018/1‑149 (MQ=255)
                                                                   gAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGAt                                                                                  >  2:280288/1‑147 (MQ=255)
                                                                      gCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGc                                                                             <  1:280288/149‑1 (MQ=255)
                                                                          cTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGc                                                                         <  2:114714/149‑1 (MQ=255)
                                                                              ccAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATg                                                                                                               <  2:164818/107‑1 (MQ=255)
                                                                              ccAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCgat                                                                      >  4:309043/1‑148 (MQ=255)
                                                                                   gAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGc                                                                <  3:135498/149‑1 (MQ=255)
                                                                                    aGCCGAAGGTATGCTGCACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCg                                                               <  2:290811/149‑1 (MQ=255)
                                                                                        gAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGc                                                           >  7:114079/1‑149 (MQ=255)
                                                                                          aGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCgg                                                         <  1:136581/149‑1 (MQ=255)
                                                                                               tGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAAt                                                    <  3:309043/149‑1 (MQ=255)
                                                                                                  tGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGac                                                                                           <  8:274068/107‑1 (MQ=255)
                                                                                                  tGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGac                                                                                           <  8:269747/107‑1 (MQ=255)
                                                                                                  tGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGac                                                                                           >  7:269747/1‑107 (MQ=255)
                                                                                                  tGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGac                                                                                           >  7:274068/1‑107 (MQ=255)
                                                                                                    tACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCCGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGcc                                               >  3:189917/1‑149 (MQ=255)
                                                                                                             aTGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACt                                      >  4:124766/1‑149 (MQ=255)
                                                                                                             aTGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACt                                      >  1:189114/1‑149 (MQ=255)
                                                                                                                 agcaCCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGa                                   >  1:166174/1‑148 (MQ=255)
                                                                                                                   caCCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCggg                                                        <  2:263150/125‑1 (MQ=255)
                                                                                                                        ccgcACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACcgcgc                           >  7:14888/1‑149 (MQ=255)
                                                                                                                          gcACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCgcga                         <  8:14888/149‑1 (MQ=255)
                                                                                                                               gTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCg                                                               <  4:77504/106‑1 (MQ=255)
                                                                                                                               gTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCg                                                               >  3:77504/1‑106 (MQ=255)
                                                                                                                                  ggtggCATCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATg                  <  4:172781/148‑1 (MQ=255)
                                                                                                                                       caccAGCCACGACATCGCGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGc                                                                <  6:302763/94‑1 (MQ=255)
                                                                                                                                       caccAGCCACGACATCGCGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGc                                                                >  5:302763/4‑97 (MQ=255)
                                                                                                                                        gTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGTCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCgcgag                        <  2:104693/136‑1 (MQ=255)
                                                                                                                                        gTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCgcgag                        >  1:104693/1‑136 (MQ=255)
                                                                                                                                        accAGCCACGACATCGCGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCAtt             <  6:208529/145‑1 (MQ=255)
                                                                                                                                        accAGCCACGACATCGCGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCAtt             >  5:208529/3‑147 (MQ=255)
                                                                                                                                           aGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCTATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACaa        <  2:89279/149‑1 (MQ=255)
                                                                                                                                            gCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAAt       >  2:78160/1‑149 (MQ=255)
                                                                                                                                            gCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAAt       <  8:114079/149‑1 (MQ=255)
                                                                                                                                                 gATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGAc    <  6:101265/147‑1 (MQ=255)
                                                                                                                                                 gATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGAc    >  5:101265/1‑147 (MQ=255)
                                                                                                                                                 gATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGa  >  8:229725/1‑149 (MQ=255)
                                                                                                                                                 gATGTCGTGGCCGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGa  <  4:189917/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACTCGGCGCAGGCACAGGTCAGGCCAACCGCGGTCACCAGGCAGGCCAGCGCGATCAGCACGGCCAGGAAGCCGCTGCCCAGGGAGCCGAAGGTATGCTGTACGTAGGCATGCAGCACCGCCGCACCGTTGGTGGCGTCAGCCGCGATGTCGTGGCTGGTGGCACCCAGGCGGAACAGGCTGATGTAGACCAGGGCCAGCCCGACACCGGCGATCAGGCCGGCGATGATGGCGTAGCGGGTAATCAGCCTCGGCGACTCGACACCGCGCGAGCGAATGGCATTGACAATGACGA  >  NC_002947/4726472‑4726765

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