Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,014,259 Δ1 bp intergenic (‑20/+28) yqeB ← / ← yqeC XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold/putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,014,2580C.98.3% 232.3 / ‑2.8 58intergenic (‑19/+29)yqeB/yqeCXdhC‑CoxI family protein with NAD(P)‑binding Rossman fold/putative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (21/36);  minor base T (1/0);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGCGCACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGT  >  NC_000913/3014121‑3014389
                                                                                                                                         |                                                                                                                                   
gcgcgcACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATct                                                                                                                                    >  1:146600/1‑138 (MQ=255)
 cgcgcACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATctt                                                                                                                                   >  2:1176451/1‑137 (MQ=255)
  gcgcACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATcttt                                                                                                                                  <  1:1243262/139‑4 (MQ=255)
    gcACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCa                                                                                                                               <  2:209279/139‑1 (MQ=255)
      aCTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGt                                                                                                                             <  2:911450/139‑1 (MQ=255)
             tttGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCgg                                                                                                                      <  1:257515/139‑1 (MQ=255)
               tgtgCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCGACCTCAATGGATATCTTTCAGTAACCCGGaa                                                                                                                    <  2:146600/139‑1 (MQ=255)
                        tGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTa                                                                                                                            <  2:532462/122‑1 (MQ=255)
                        tGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTa                                                                                                                            >  1:532462/1‑122 (MQ=255)
                        tGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGc                                                                                                           <  2:398307/139‑1 (MQ=255)
                           cGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgt                                                                                                        <  1:1127999/139‑2 (MQ=255)
                           cGGAGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAATCGGAGAATTTTGCTCTTCGTGTTTTGCCGCCTCTGTGAAAATATTCATACCAACCTCAATGGTTTTCTTTCACTAACCCGGAATACCCGGGCCgt                                                                                                        >  2:57919/1‑138 (MQ=255)
                             gggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                      <  1:628150/139‑4 (MQ=255)
                             gggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                      <  1:352988/139‑4 (MQ=255)
                               gAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAAt                                                                                                    <  1:143286/139‑1 (MQ=255)
                                aGTCGAGCCACGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAAtt                                                                                                   >  1:1225935/1‑139 (MQ=255)
                                  tCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTca                                                                                                 >  1:821108/1‑139 (MQ=255)
                                  tCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTca                                                                                                 >  2:1024512/1‑139 (MQ=255)
                                              tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                     <  1:1024512/139‑1 (MQ=255)
                                                 caacaaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGa                                                                                  <  1:945029/139‑1 (MQ=255)
                                                       tttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAg                                                                            >  1:1144449/1‑139 (MQ=255)
                                                        ttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                           <  2:821108/139‑1 (MQ=255)
                                                         tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          <  2:225294/139‑1 (MQ=255)
                                                         tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          <  2:1144449/139‑1 (MQ=255)
                                                         tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          <  1:278918/139‑1 (MQ=255)
                                                          gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCg                                                                                                              >  1:408231/1‑102 (MQ=255)
                                                          gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCg                                                                                                              <  2:408231/102‑1 (MQ=255)
                                                           cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        <  1:388610/139‑1 (MQ=255)
                                                                          aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        <  2:579402/124‑1 (MQ=255)
                                                                          aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        >  1:579402/1‑124 (MQ=255)
                                                                            aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCGTACCAACCTAAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                           >  2:811055/1‑119 (MQ=255)
                                                                            aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCCTACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                           <  1:811055/119‑1 (MQ=255)
                                                                              ttttGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGCGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                     <  1:57919/139‑1 (MQ=255)
                                                                                  gCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGa                                                 <  1:817733/139‑1 (MQ=255)
                                                                                        cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        >  1:485360/1‑110 (MQ=255)
                                                                                        cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        <  2:485360/110‑1 (MQ=255)
                                                                                             tttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                      <  2:832889/139‑1 (MQ=255)
                                                                                             tttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                      <  1:1015801/139‑1 (MQ=255)
                                                                                                        tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCgtta                                                                                                      <  1:499303/64‑4 (MQ=255)
                                                                                                        tgtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                      >  2:499303/1‑64 (MQ=255)
                                                                                                         gtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                           <  1:793553/90‑1 (MQ=255)
                                                                                                         gtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          >  2:151693/1‑91 (MQ=255)
                                                                                                         gtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          <  1:151693/91‑1 (MQ=255)
                                                                                                         gtgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCCCGATTGCAgg                                                                           >  2:793553/1‑90 (MQ=255)
                                                                                                          tgAAAATATTCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                         <  2:1310063/139‑1 (MQ=255)
                                                                                                                  ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          >  2:1224380/1‑82 (MQ=255)
                                                                                                                  ttCATACCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                          <  1:1224380/82‑1 (MQ=255)
                                                                                                                       aCCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                      >  1:1315292/1‑113 (MQ=255)
                                                                                                                       aCCAACCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                      <  2:1315292/113‑1 (MQ=255)
                                                                                                                          aaCCTCGATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGcc                                                         <  2:865194/91‑1 (MQ=255)
                                                                                                                          aaCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGcc                                                         >  1:865194/1‑91 (MQ=255)
                                                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        >  1:174879/1‑75 (MQ=255)
                                                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                        <  2:174879/75‑1 (MQ=255)
                                                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg        >  1:371358/1‑139 (MQ=255)
                                                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATAGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                            <  2:12369/103‑1 (MQ=255)
                                                                                                                           aCCTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATAGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                            >  1:12369/1‑103 (MQ=255)
                                                                                                                            ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                            >  2:1266238/1‑102 (MQ=255)
                                                                                                                            ccTCAATGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                            <  1:1266238/102‑1 (MQ=255)
                                                                                                                                 aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                           <  1:818557/66‑1 (MQ=255)
                                                                                                                                 aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                           >  2:818557/1‑66 (MQ=255)
                                                                                                                                 aTGGATATCTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCCACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgt  >  1:990419/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                   
GCGCGCACTAACATTTGTGCAGAATGGCGGGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGT  >  NC_000913/3014121‑3014389

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