Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,714,397 C→T E666E (GAG→GAA bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,3970CT95.8% 72.5 / NA 24E666E (GAG→GAAbisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/14);  minor base A (1/0);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCG  >  NC_000913/3714275‑3714524
                                                                                                                          |                                                                                                                               
ggTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTGCTGGGCATGGTCTGGAGGAATGGTGGCAGGATCACGATTTg                                                                                                                       >  2:159141/1‑133 (MQ=255)
ggTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTg                                                                                                                       <  2:451491/133‑1 (MQ=255)
          ccTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTa                                                                                                             >  1:271395/1‑133 (MQ=255)
            tCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACa                                                                                                           >  2:53289/1‑133 (MQ=255)
            tCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCACGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACa                                                                                                           >  2:20383/1‑133 (MQ=255)
               cTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAAtt                                                                                                        <  2:403340/133‑1 (MQ=255)
                     aCCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCa                                                                                                  <  2:113217/133‑1 (MQ=255)
                        gCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCagag                                                                                               <  2:85907/133‑1 (MQ=255)
                         cTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAgaga                                                                                              >  2:420326/1‑133 (MQ=255)
                          tCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGaa                                                                                             <  2:419103/133‑1 (MQ=255)
                                        gCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTa                                                                                        <  2:103900/124‑1 (MQ=255)
                                        gCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTa                                                                                        >  1:103900/1‑124 (MQ=255)
                                         cccGCGTGAGTTCCACAAGCGACCAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAgctgg                                                                              <  1:221598/133‑1 (MQ=255)
                                                              aaCAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGa                                                         <  1:409359/133‑1 (MQ=255)
                                                              aaCAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGa                                                         >  2:409359/1‑133 (MQ=255)
                                                                  gTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATggg                                                     <  1:6221/133‑1 (MQ=255)
                                                                         cATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAg                                              <  1:420326/133‑1 (MQ=255)
                                                                           tCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGCTGCACCAGATGGGCAGAAAGTa                                            >  1:246171/1‑133 (MQ=255)
                                                                             ttGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTAcc                                          <  1:273564/133‑1 (MQ=255)
                                                                             ttGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTAcc                                          <  2:54556/133‑1 (MQ=255)
                                                                                       cgcTCCTGGGCATCGTCAGGATGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTACCTGCAACTGtt                                >  1:501730/1‑133 (MQ=255)
                                                                                                           aTGAATGGTGACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGc            <  1:20383/133‑1 (MQ=255)
                                                                                                                    gACAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTcc   >  2:295907/1‑133 (MQ=255)
                                                                                                                     aCAGGTTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCg  <  2:36253/133‑1 (MQ=255)
                                                                                                                          |                                                                                                                               
GGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGCCCGCGTGCGTTCCACAACCGAACAGTATCGCCATCTTGTATGCCGCGCTCCTGGGCATCGTCAGGATGAATGGTGACAGGCTCACGATTTGCCACCGCGTACAATTCGCGCAGAGAACTGTAATTCAGCTGGCTGTGCAGGCGGTGCGCCGGATGGGCAGAAAGTACCTGCAACTGTTCTGGTTCGGCATTGCCCTGCCATTCGTCCG  >  NC_000913/3714275‑3714524

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