Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,097,938 G→T 7.3% T251K (ACG→AAG) ‡ qor‑II ← quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,097,9380GT 7.3% 171.6 / 3.9 55T251K (ACG→AAG) ‡qor‑IIquinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base T (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GGGAAGCATCGCCAGACCCGAGGCTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACACCTTGGCCTTGGG  >  NC_002947/4097791‑4098079
                                                                                                                                                   |                                                                                                                                             
gggAAGCATCGCCAGACCCGAGGCTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGt                                                                                                                                              >  7:92592/1‑149 (MQ=255)
 ggAAGCATCGCCAGACCCGAGGCTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTg                                                                                                                                             <  7:216250/149‑1 (MQ=255)
              gACCCGAGGCTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATgg                                                                                                                                <  3:53453/149‑1 (MQ=255)
                 ccGAGGCTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAg                                                                                                                             >  7:287152/1‑149 (MQ=255)
                       cTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTg                                                                                                                       <  3:259642/149‑1 (MQ=255)
                        tGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTg                                                                                                                       >  4:170137/1‑148 (MQ=255)
                              aTACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGGGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGa                                                                                                                <  2:318333/149‑1 (MQ=255)
                               tACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAg                                                                                                               >  1:108593/1‑149 (MQ=255)
                                 ccTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGt                                                                                                             <  3:170137/149‑1 (MQ=255)
                                        cAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTgccgcc                                                                                                       >  4:56527/1‑148 (MQ=255)
                                        cAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTgccgcc                                                                                                       >  2:196159/1‑148 (MQ=255)
                                         aaCCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGCCATTGAGGTTGCCACGCTCAATGCCaatggc                                                                                                                         <  6:323404/129‑6 (MQ=38)
                                         aaCCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGCCATTGAGGTTGCCACGCTCAATGCCaatggc                                                                                                                         >  5:323404/1‑124 (MQ=38)
                                           ccGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAAc                                                                                                   >  1:353736/1‑149 (MQ=255)
                                              cTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGCGTCAGGTACAGCGAGCCATAGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGa                                                                                                 >  3:205644/1‑148 (MQ=255)
                                                  cAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATAcc                                                                                             >  1:155241/1‑148 (MQ=255)
                                                                cgcgcTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCAcc                                                                              <  3:91103/149‑1 (MQ=255)
                                                                      gctgGTACTCCTCGACAGTGCGGGTGTGTGCTGGGAGGGATGGCCGGGTCAGGTTCAGCGTGCCTTTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAgg                                                                        <  4:205644/149‑1 (MQ=255)
                                                                       ctgGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAgg                                                                        <  8:215881/148‑1 (MQ=255)
                                                                         ggTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTg                                                                     >  1:315195/1‑149 (MQ=255)
                                                                          gTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATa                                                                                               >  8:324466/1‑122 (MQ=255)
                                                                          gTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATa                                                                                               <  7:324466/122‑1 (MQ=255)
                                                                           gacgcctcGACAGTGCGGGGGTGTGCGGCGAGCGATGGCCGGGGCAGGTACAGCGAGCCATGGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGcc                                                                   <  6:104/145‑1 (MQ=255)
                                                                            aCTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCTGAACGATACCAACAAGCCGCGCACCCTCAGGCTGcc                                                                   >  5:287614/1‑148 (MQ=255)
                                                                                ctcGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGTGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAgg                                                                        >  3:63771/1‑139 (MQ=255)
                                                                                ctcGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGTGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAgg                                                                        <  4:63771/139‑1 (MQ=255)
                                                                                ctcGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGAt                                                                                                <  6:241909/115‑1 (MQ=255)
                                                                                ctcGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGAt                                                                                                >  5:241909/1‑115 (MQ=255)
                                                                                  cGACAGGGCGGGTGTGTGCTGCGAGGGATGGCCGGGTCTGGTACAGCGAGCCATTGGCATTGTGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGTGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAg                                                            <  4:100368/149‑1 (MQ=255)
                                                                                   gACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGa                                                             <  3:56527/147‑1 (MQ=255)
                                                                                          cGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAg                                                    <  1:196159/149‑1 (MQ=255)
                                                                                                 gtgCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACAcc                                             <  2:155241/149‑1 (MQ=255)
                                                                                                       cgagcgaTGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGcc                                        <  2:353736/148‑1 (MQ=255)
                                                                                                        gagcgaTGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCCGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCa                                      >  2:111262/1‑149 (MQ=255)
                                                                                                              tGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCTGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAAt                                <  6:287614/149‑1 (MQ=255)
                                                                                                                  cGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGTACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTa                            <  1:327604/149‑1 (MQ=255)
                                                                                                                  cGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGTACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTa                            >  2:327604/1‑149 (MQ=255)
                                                                                                                  cGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCaacc                                                   >  6:127258/1‑124 (MQ=255)
                                                                                                                  cGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGc                                                   <  5:127258/126‑1 (MQ=255)
                                                                                                                      tCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTacac                         >  6:121938/1‑148 (MQ=255)
                                                                                                                          gTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACgg                    <  5:121938/149‑1 (MQ=255)
                                                                                                                          gTACAGCGAGCCATCGGCATTGAGCGTGGCAACCCCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACgg                    >  8:330281/1‑149 (MQ=255)
                                                                                                                           tACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGc                   <  8:287152/149‑1 (MQ=255)
                                                                                                                            aCAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGc                   >  5:251008/1‑148 (MQ=255)
                                                                                                                              cgCGAGCCATTCGCATTGAGCGGGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCcac                <  6:251008/148‑1 (MQ=255)
                                                                                                                               gCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCcaca               >  3:260819/1‑149 (MQ=255)
                                                                                                                                 gAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAgg                                                                        >  7:225744/1‑90 (MQ=255)
                                                                                                                                 gAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAgg                                                                        <  8:225744/90‑1 (MQ=255)
                                                                                                                                    ccATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACAccttg          <  4:260819/149‑1 (MQ=255)
                                                                                                                                         gacATTGAGGTTGCCACGCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACAccttggcctt     >  4:139941/3‑149 (MQ=255)
                                                                                                                                          gCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACAccttggccttg    >  6:158068/1‑149 (MQ=255)
                                                                                                                                          gCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACAccttggccttg    >  2:152717/1‑149 (MQ=255)
                                                                                                                                           cATTGAGGTTGCCACGCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACAccttggccttgg   <  3:139941/149‑1 (MQ=255)
                                                                                                                                           cATTGAGCGTGGCAACCCCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACAccttggccttgg   <  7:330281/149‑1 (MQ=255)
                                                                                                                                             ttGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACACCTTGGCCTTggg  >  8:197143/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GGGAAGCATCGCCAGACCCGAGGCTGGACGATACCTTCGGCAACCGCTGCCAGCACATCGGCCGCGCGCTGCTGGTACTCCTCGACAGTGCGGGTGTGTGCTGCGAGCGATGGCCGGGTCAGGTACAGCGAGCCATTGGCATTGAGCGTGGCAACCTCAATGGCAGGCGGTGCACCGGAGGTGCCGCCGAACGATACCAACAAGCCGCGCACCCTCAGGCTGCCCAACGAGGCCTCAAGCGACACCTTGCCCACGCAATCGTACACTACGGCCACACCTTGGCCTTGGG  >  NC_002947/4097791‑4098079

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: