Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 137,824 A→G T248A (ACC→GCC)  cueO → cuprous oxidase/multicopper oxidase with role in copper homeostasis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913137,8240AG97.1% 120.8 / ‑4.8 35T248A (ACC→GCC) cueOcuprous oxidase/multicopper oxidase with role in copper homeostasis
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (12/22);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCACCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGAT  >  NC_000913/137592‑138047
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               
aaTGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCg                                                                                                                                                                                                                           >  1:396287/1‑239 (MQ=255)
       ccAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGc                                                                                                                                                                                                                   >  2:690227/1‑240 (MQ=255)
          aaaCAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCg                                                                                                                                                                                                                <  1:341410/240‑1 (MQ=255)
                 ggggTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTAt                                                                                                                                                                                                          >  2:704982/1‑239 (MQ=255)
                     tATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGAt                                                                                                                                                                                                     >  2:78413/1‑240 (MQ=255)
                          atgatgTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGcc                                                                                                                                                                                                                  >  2:265836/1‑222 (MQ=255)
                          atgatgTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGcc                                                                                                                                                                                                                  <  1:265839/222‑1 (MQ=255)
                            gatgTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCg                                                                                                                                                                                                                <  2:235448/222‑1 (MQ=255)
                            gatgTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAg                                                                                                                                                                                               <  1:245061/239‑1 (MQ=255)
                                       gATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCg                                                                                                                                                                                                                <  1:85289/211‑1 (MQ=255)
                                          cGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGc                                                                                                                                                                                 >  2:211133/1‑239 (MQ=255)
                                          cGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGc                                                                                                                                                                                 >  1:456084/1‑239 (MQ=255)
                                             tCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTAcc                                                                                                                                                                             <  1:258546/240‑1 (MQ=255)
                                               aGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCt                                                                                                                                                                            >  2:258976/1‑239 (MQ=255)
                                                ggATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGa                                                                                                                                                                          <  1:542382/240‑1 (MQ=255)
                                                       aaaTTTAGCGCGGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGt                                                                                                                                                                    <  1:315583/239‑1 (MQ=255)
                                                                                           gATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCg                                                                                                                                >  1:317528/1‑239 (MQ=255)
                                                                                                 atgaCCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAAGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGa                                                                                                                          <  2:103880/239‑1 (MQ=255)
                                                                                                           ccgTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCtggtgg                                                                                                               <  1:144222/240‑1 (MQ=255)
                                                                                                               gggCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAgg                                                                                                            <  2:396282/239‑1 (MQ=255)
                                                                                                                   tGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTaa                                                                                                        <  1:538584/239‑1 (MQ=255)
                                                                                                                           cGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACaa                                                                                               <  1:614390/240‑1 (MQ=255)
                                                                                                                                   tGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTg                                                                                       <  1:216343/240‑1 (MQ=255)
                                                                                                                                     cTGACCAACGGTGCCATCTACCCGCTACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTTTAATGCCCGTTCGCTCCATGTCGCCACCAGCGGCAATCGCCCGCTGTATGTGATTGCCAGCGTCGGTGGTCTTCTACCTGGACCTGTGa                                                                                                                                                                  >  2:203871/1‑163 (MQ=255)
                                                                                                                                                  cAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGGGa                                                                              >  2:164173/1‑234 (MQ=255)
                                                                                                                                                          ccGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGc                                                                <  1:524207/240‑1 (MQ=255)
                                                                                                                                                                                 ttGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTgcgc                                          <  1:79171/239‑1 (MQ=255)
                                                                                                                                                                                 ttGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGcc                                         >  2:686269/1‑240 (MQ=255)
                                                                                                                                                                                   ggCTGCGCCTGCGTTTGCTCAATGGCTGTCATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCg                                        <  1:318243/239‑1 (MQ=255)
                                                                                                                                                                                      tGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCCGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACaa                                                                                               >  2:209648/1‑181 (MQ=255)
                                                                                                                                                                                      tGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCCGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACaa                                                                                               <  1:209651/181‑1 (MQ=255)
                                                                                                                                                                                             gCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCtggtg                                                                                                                <  2:425135/157‑1 (MQ=255)
                                                                                                                                                                                                                  tGCCCGTTCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCg                                                                      <  2:339370/178‑1 (MQ=255)
                                                                                                                                                                                                                        ttCGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGAt  <  1:690239/240‑1 (MQ=255)
                                                                                                                                                                                                                          cGCTCAATTTCGCCGCCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGAt  <  2:50974/238‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               
AATGCTGCCAAAACAGTGGGGTATCGATGATGTTCCGGTGATCGTTCAGGATAAGAAATTTAGCGCCGACGGGCAGATTGATTATCAACTGGATGTGATGACCGCCGCCGTGGGCTGGTTTGGCGATACGTTGCTGACCAACGGTGCAATCTACCCGCAACACGCTGCCCCGCGTGGTTGGCTGCGCCTGCGTTTGCTCAATGGCTGTAATGCCCGTTCGCTCAATTTCGCCACCAGCGACAATCGCCCGCTGTATGTGATTGCCAGCGACGGTGGTCTGCTACCTGAACCAGTGAAGGTGAGCGAACTGCCGGTGCTGATGGGCGAGCGTTTTGAAGTGCTGGTGGAGGTTAACGATAACAAACCCTTTGACCTGGTGACGCTGCCGGTCAGCCAGATGGGGATGGCGATTGCGCCGTTTGATAAGCCTCATCCGGTAATGCGGATTCAGCCGAT  >  NC_000913/137592‑138047

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