Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,092,973 G→A 7.8% G89G (GGC→GGT)  ygiD ← 4,5‑DOPA‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,092,9730GA 7.8% 150.2 / 2.5 51G89G (GGC→GGT) ygiD4,5‑DOPA‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base A (2/2);  total (25/26)
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.70e-01

ACGGGAATGCGCGCATCCGGGTACATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCGGCCGGGGCAGCG  >  NC_002947/2092836‑2093121
                                                                                                                                         |                                                                                                                                                    
acGGGAATGCGCGCATCCGGGTACATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCGAGCCCgggcttgc                                                                                                                                           <  3:362219/149‑7 (MQ=255)
          gcgcATCCGGGTACATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCAcc                                                                                                                                 >  4:71880/1‑149 (MQ=255)
            gcATCCGGGTACATGAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGc                                                                                                                               <  5:237933/149‑1 (MQ=255)
              aTCCGGGTACATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCa                                                                                                                              >  6:372496/1‑148 (MQ=255)
                      aCATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCt                                                                                                                     >  5:253483/1‑149 (MQ=255)
                       cATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCtg                                                                                                                    <  3:71880/149‑1 (MQ=255)
                        aTCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCtgg                                                                                                                   >  4:280400/1‑149 (MQ=255)
                        aTCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCtgg                                                                                                                   >  1:220776/1‑149 (MQ=255)
                        aTCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCACCCGGCGCCGGGTACTGCACCGCATACAGGGCtgg                                                                                                                   >  2:130315/1‑149 (MQ=255)
                        aTCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCACCCGGCGCCGGGTACTGCACCGCATACAGGGCtgg                                                                                                                   >  2:135516/1‑149 (MQ=255)
                         tCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCtggt                                                                                                                  <  5:372496/149‑1 (MQ=255)
                          cAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTg                                                                                                                                     <  8:279184/129‑1 (MQ=255)
                          cAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTg                                                                                                                                     >  7:279184/1‑129 (MQ=255)
                              gACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCa                                                                                                                              <  6:143476/132‑1 (MQ=255)
                              gACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCa                                                                                                                              >  5:143476/1‑132 (MQ=255)
                                        cccAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAg                                                                                                   >  6:260056/1‑149 (MQ=255)
                                        cccAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGCCCAGCAAGCCAGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAg                                                                                                   <  3:193973/149‑1 (MQ=255)
                                         ccAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGcc                                                                                                          >  7:132787/1‑141 (MQ=255)
                                         ccAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGcc                                                                                                          <  8:132787/141‑1 (MQ=255)
                                         ccAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGt                                                                                                  >  6:170426/1‑149 (MQ=255)
                                          cAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCACCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTc                                                                                                 <  1:135516/149‑1 (MQ=255)
                                          cAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCACCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTc                                                                                                 <  1:130315/149‑1 (MQ=255)
                                           aGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCa                                                                                                >  7:160587/1‑149 (MQ=255)
                                           aGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCa                                                                                                >  6:71790/1‑149 (MQ=255)
                                              cgcCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGc                                                                                             <  3:280400/149‑1 (MQ=255)
                                                        gAATGGGCGTTGCATGTCCAGGCGCGCGGGTAACCGCGCAGCCTTCAACCGCTCACTGACCTGGCGAGCAAGCCCGGGCTGGCCCGGCGGCGGGTTCTGCACGGCATACAGGGCTGGTGGGAAGCCGTAGAATTCATGCCAGGTTTCg                                                                                    <  1:344220/148‑1 (MQ=255)
                                                        gAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCg                                                                                    <  5:170426/148‑1 (MQ=255)
                                                            gggCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAg                                                                                                   >  4:329482/1‑129 (MQ=255)
                                                            gggCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAg                                                                                                   <  3:329482/129‑1 (MQ=255)
                                                               cGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTgc                                                                            <  7:102905/149‑1 (MQ=255)
                                                                   gCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTg                                                                        >  1:195578/1‑149 (MQ=255)
                                                                     aTGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGcc                                                                      <  8:160587/149‑1 (MQ=255)
                                                                     aTGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGcc                                                                      >  7:351520/1‑149 (MQ=255)
                                                                     aTGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGcc                                                                      <  6:253483/149‑1 (MQ=255)
                                                                          cAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTg                                                                  <  5:260056/148‑1 (MQ=255)
                                                                          cAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGa                                                                 >  4:341906/1‑149 (MQ=255)
                                                                           aGCCGCGCGTGTAACCCCGCAGCCTTCAACAGCTCACTGACGTGGCCAGCAAGCCCTGGCTCGCCCGGCGCCGGGTACTGCACCTCATACAGGTCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGAc                                                                <  2:28271/149‑1 (MQ=255)
                                                                           aGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGAc                                                                >  2:195919/1‑149 (MQ=255)
                                                                                cgcgGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACcagca                                                           <  5:280759/149‑1 (MQ=255)
                                                                                   gggTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCt                                                        >  6:337730/1‑149 (MQ=255)
                                                                                         ccccGCAGCCTTCAACCGCTCACTGACCTGGCCAGCACGCCCGGGCTCGCCCGGCGCCGGGTACTGCCCCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGc                                                  >  4:227141/1‑149 (MQ=255)
                                                                                                  cTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAgt                                         <  2:195578/149‑1 (MQ=255)
                                                                                                     cAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAgtgt                                       <  5:71790/148‑1 (MQ=255)
                                                                                                            tCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGacac                               >  6:330990/1‑149 (MQ=255)
                                                                                                                       gCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCg                    <  1:195919/149‑1 (MQ=255)
                                                                                                                       gCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCg                    >  4:366375/1‑149 (MQ=255)
                                                                                                                       gCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCg                    <  3:341906/149‑1 (MQ=255)
                                                                                                                           gCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCttt                >  8:29512/1‑149 (MQ=255)
                                                                                                                             aaGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCg              <  3:366375/149‑1 (MQ=255)
                                                                                                                                cccGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCGGcc           <  2:220776/149‑1 (MQ=255)
                                                                                                                                cccGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCGGcc           >  5:163514/1‑149 (MQ=255)
                                                                                                                                         gCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCGGCCGGGGCAgcg  >  4:40712/1‑149 (MQ=255)
                                                                                                                                         gCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCGGCCGGGGCAgcg  >  4:32923/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
ACGGGAATGCGCGCATCCGGGTACATCAGCGACAGCGGCACCCAGGCGCCATGGTCGAATGGGCGTTGCATGTCCAGCCGCGCGGGTAACCCCGCAGCCTTCAACCGCTCACTGACCTGGCCAGCAAGCCCGGGCTCGCCCGGCGCCGGGTACTGCACCGCATACAGGGCTGGTGGGAAGCCGTAGAAGTCATGCCAGGTTTCGGGCTGTGCGCTGCCGGTGACCAGCAGCTCGCGGCTTTCCCAGTGTGCCGACACCACCACGATCGCTTTCGGCCGGGGCAGCG  >  NC_002947/2092836‑2093121

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: ‑